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ANDRZEJ SOLTAN 
(1897—1959) 


The 10th of December 1959 Polish physics suffered a tremendous loss by the 
sudden death of Professor Andrzej (Andrew) Soltan, a leading Polish nuclear phy- 
sicist. This quite unexpected passing away of Professor A. Soltan occured just in 
the period when Polish physics, especially the nuclear physics in this country, ur- 
gently needed his care and guidance. 

Everybody who knew or only approached Professor Soltan was not only under the 
influence of his noble and attractive personality, but also full of admiration for the 
clarity and vividness of his mind, for the amplitude and depth of his knowledge 
in all, even very remote, branches of physics and for his excellent experimental 
talents. 

Andrzej Sottan was born in Warsaw the 25th of November 1897. He completed 
his secondary education in Petrograd (Russia) during World War I. After his return 
to Poland he began to study physics at the Warsaw University. Very early, even 
before being graduated, Andrzej Soltan became assistant at the Institute of Experi- 
mental Physics of the Warsaw University and worked in close collaboration with 
Professor Stefan Pieńkowski who just iniciated the organization of his well-known 
research school in molecular physics. Simultaneously with his numerous didactic 
and administrative duties in the Institute, Andrzej Soltan began under the leadership 
of Prof. Pieńkowski experimental research work on the band spectrum of mercury. 
In particular this work concerned the emission spectrum of mercury hydrate mole- 
cules under electrical discharges. The first results of his investigations Andrzej 

- — Soltan reported in 1924 at the II General Meeting of the Polish Physical Society 
held in Cracow. In 1925 and 1926 two papers appeared in the “Comptes Rendus 
des Séances de la Société Polonaise de Physique” containing the results of these 
and subsequent works on the same subject. Presenting his doctor thesis un band 
spectra in mixtures of hydrogene and mercury vapour, Andrzej Sean un the 
Ph. D. degree in 1926. Shortly afterwards he published together with S. Szezeniowski 
in the same Proceedings of the Polish Physical Society a paper entitled “Caracteristi- 
ques de sensibilité de quelques émulsions photographiques". 

The following year, due to the support of Prof. Pieńkowski, Dr Soltan was 

— granted a fellowship enabling him to spend a year in pone on research work in X-ray 

spectroscopy at the laboratory of Maurice de Bun Together with He} Thibaut 
they worked on the so called «intermediate domain" and succeeded to fill this empty 
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gap in the experimental spectrum of electromagnetic radiation between the farthest 
ultraviolet and the region of the softest then known X-rays. In these investigations 
they applied a method of ruled diffraction gratings and beams of X-rays incident at 
nearly glancing angles. Thus this disquieting gap in the electromagnetic spectrum 
has been experimentally closed. The conspicuous results of these investigations of 
Soltan and Thibaut have been published in the “Comptes Rendus des Séances de 
l'Académie des Sciences de Paris" and in the *Journal de Physique et le Radium". 

After his return to Warsaw Dr Soltan is again very active, mainly in organiza- 
tion and administrative work, and in enlarging and modernizing the Institute of Experi- 
mental Physics of the Warsaw University. In this work he became again the main 
collaborator of Prof. Pienkowski. Due to these extensive efforts the Institute was 
transformed into one of the best organized and equipped physical laboratories in 
Europe. In spite of his extreemly absorbing occupations Dr Soltan still had some 
time for experimental research work. The results of this work were published in two 
papers which appeared in the “C. R. Soc. Polon. Phys." and in the “Acta Physica 
Polonica". The first of these papers contains the description of a new model of 
a capillary mercury lamp constructed by Dr Soltan. This type of mercury lamp 
proved to be very convenient in different optical investigations and is still widely 
used in many spectroscopical laboratories. In the second paper Dr Soltan reports 
the results of his newly commenced experiments carried on with unusual tenacity 
and aiming at discovering the influence of chemical bonds on the absorption coef- 
ficients of X-rays. 


In 1933 Dr Soltan received a fellowship of the Rockefeller Foundation enabling . 


him to spend a year on research work in the Kellog Laboratory, California Institute 
of Technology, Pasadena. There he started his new line of research work — nuclear 
physies, to which he remained faithful through the rest of his life. In this so dramatic 
and important period in the development of nuclear physics, soon after the first 
completely artificial transmutation of elements had been performed, and the neutron 
and the positive electron detected, Dr Soltan entered a team doing experiments 
on nuclear reactions caused by artificially accelerated ions. In collaboration with C. C. 
Lauritsen and H. R. Crane results were obtained of fundamental importance for the 
further progress of nuclear physics. Crane, Lauritsen and Soltan have proved the 
emission of neutrons by different targets of light elements bombarded by high speed 
deuterons, this being the first experiment on the production of neutrons by artifi- 
cially accelerated particles. The results of these important researches were published 
in four papers, one in the “C. R. Acad. Sci. de Paris" and three in the “Physical Review". 

Back in Warsaw, Dr Soltan began immediately the difficult task of creating in 
the Warsaw Institute of Experimental Physics of the University experimental facili- 
ties enabling him to continue the research work started in Pasadena. In spite of the 
tremendous amout of administrative duties at the post of a chief assistant (adjunkt) at 
the quickly growing Institute, he succeeded to build, almost with one’s own hand, 
a cascade ion accelerator of the Greinacher type up to energies of about 800 kV. Using 


| 
| 
i 
x 


Andrzej Soltan (1897—1959) 10 


this accelerator as a source of fast neutrons, he performed various experiments on nu- 
clear reactions produced by these neutrons. In these investigations Dr Soltan was help- 
ed by Prof. Ludwik Wertenstein, director of the Radiological Laboratory of the 
Warsaw Scientific Society, and Dr Jan Cichocki, who at the end 1938 replaced 
Dr Soltan at the post of chief assistant at the Institute of Physics of the Warsaw 
University. The results of these works were published in four papers in 1937, 1938, 
and 1939, one of them in common with L. Wertenstein and two in common with 
J. Cichocki. 

In 1938 on the ground of the paper on radio-isotopes of bromine Dr A. Soltan 
habilitated at the Warsaw University and was promoted to the rank of docent 
(assistant professor) of experimental physics at this university. 

About a year before the outbreak of World War II, Dr Soltan gave up his post 
of a chief assistant at the Institute of Experimental Physics, Warsaw University, 
and became director of the physical research laboratory in the Philips works in War- 
saw. Continuing his research work in the Institute of Experimental Physics, Warsaw 
University, he started simultaneously in his new, reorganized and enlarged, technical 
laboratory to carry into effect a project of building a small cyclotron. Unfortunately 
the fulfilment of this ambitious enterprise was interrupted by the outbreak of World 
War ll. 

During the German occupation in Poland Dr Soltan was intensively engaged 
in the secret teaching of physics in the Warsaw underground University. In 1944 
with the whole Warsaw Philips Works Dr Soltan was evacuated to Vienna, as well 
as his cyclotron, which has been afterwards completely destroyed during the air 
bombardment of Vienna. 

When the Soviet troops entered Vienna, Dr Soltan returned to Poland. Since 
Warsaw was at that time completely devastated, Dr Soltan started his work of rebuild- 
ing higher education in Poland on the post of professor of physics and director of 
the Physical Laboratory in the Institute of Technology in Łódż. In a very short time 
he performed there the great and difficult task of organizing lectures and student 
classes in physics. 

At the same time the rebuilding of the Warsaw University together with its Ins- 
titute of Experimental Physics carried out under the most energetic leadership of Pro- 
fessor S. Pieńkowski has made surprisingly quick progress. À new chair of atomic 
physics has been created in the Warsaw University and Professor Soltan has been 
invited to occupy it. Soon he moved to Warsaw retaining still for some years the 
directorship of the Physical Laboratory at the Lodz Institute of Technology. 

In 1946 on the invitation of the U. S. Governement Professor Soltan together 
with Professor Pienkowski represented Polish science at the atomic bomb tests at 
the Bikini atoll in the Pacific Ocean. During both his crossings of the American 
continent Professor Soltan renewed his pre-war bonds with several American nuclear 
this was of no minor importance for his further work. 


hysicists ; i. 
od t efforts of Professor Soltan were directed towards the creation in the 
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Institute of Experimental Physics of the Warsaw University of contemporary research 
facilities in nuclear physics. The Institute was completely rebuilt and substantially 
enlarged, a special new large hall was built to house a big high voltage accelerator, 
by a special financial support of the Governement Professor Soltan was enabled to 
order in the known Swiss firm E. Haefely, Basel, a one million volt cascade generator 
with an accelerating tube. This generator was mounted in the Warsaw Institute in 
1950. Then Professor Soltan started a very intensive programm of work consisting in 
preparing and providing additional apparatus required for research work and, what 
was certainly his most important job, in teaching and preparing a sufficient of young 
physicists to carry on the research work in nuclear physics. He was successful to 
prepare in Warsaw in the course of several years a group of capable and well trained 
young physicists on which the further development of the research work in nuclear 
physics could be based. 

In 1950 Professor Soltan was elected corresponding member of the Polish 
Academy of Science and Letters (Polska Akademia Umiejętności) in Cracow. In 
1952, he became corresponding member of the iu created Polish Academy of 
Sciences in Warsaw. 

In 1952 the Institute of Physics of the Polish Kon of Sciences began its 
activities and Professor Soltan became the director of its department of nuclear 
physics (called “Laboratory of Radioactive Isotopes”), which was housed together 
with his Laboratory of Atomic Physics in the building of the Institute of Experi- 
mental Physics of the Warsaw University. The creation of the new laboratory gave 
Professor Soltan more financial means and more scientific and auxiliary personnel 
to extend the possibilities of the experimental work in nuclear physics. 

Much more advanced possibilities for research work in nuclear physics in Poland 
arose after the establishment in 1955 of the Institute of Nuclear Research. Prof. 
Soltan’s Laboratory of Radioactive Isotopes was transfered from the Institute of 
Physics of the Polish Academy of Sciences, to the newly founded Institute. Prof. 
Soltan became the first Director General of the Institute of Nuclear Research. The 
organization of this multispecialized Institute left him only very little time to devote 
himself to nuclear physics. Nevertheless, he was very successful in organizing both 
in the building of the Institute of Experimental Physics at the University and later 
in the Centre of Nuclear Research in Swierk near Warsaw model laboratories for 
research work in nuclear physics of low and medium energies. At Swierk this work is 
mainly concentrated on nuclear neutron reactions, basing on the research reactor 
EWA, which started to operate in June 1958, and on nuclear spectroscopy. A d + t 
fast neutron generator is now being put into operation and a 10 MeV linear proton 
accelerator is approaching its completion. 

About a year before his death Professor Soltan resigned from the post of the 
Director General of the Institute of Nuclear Research in order to have more time for 
individual research work in nuclear physics. He retained the post of the Director 
of the Physics Division of the Institute, as well as that of the Director of the Centre . 
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of Nuclear Research at Swierk. At the same time he still remained President of the 
Scientific Council of the Institute of Nuclear Research. Professor Sottan always payed 
vivid attention to every work in physics done in the Institute and was personally le- 
ading many research works of his numerous pupils and collaborators. A most excellent 
and universal specialist in the techniques of experimental work in various branches of 
physics, he was always an unsurpassed and generally admired adviser to anybody 
looking for his help. E.g. his approval of any mechanical shop work, in which he 
was also an expert, was always regarded as a highest praise. 

In the course of the last two years of his life his personal interest was mostly 
directed towards plasma phenomena and thermonuclear reactions. He led in the Ins- 
titute a small but entusiastic group of research workers who did not only make theo- 
retical studies but performed already some small scale experiments, on these subjects. 

Besides many other official duties Professor Soltan was an active member of 
both the State Council for the Application of Nuclear Energy and of the Committee for 
the Peaceful Applications of the Atomic Energy of the Polish Academy of Sciences. 
In the period 1951—1952 he was also acting as President of the Polish Physical 
Society. 

Professor Soltan was a member of the Scientific Council of the Joint Institute of 
Nuclear Research at Dubna (USSR) from the moment of its establishment. In this 
capacity Professor Soltan took an active part in all sessions of the Council, exerting 
a beneficial influence on the development and activity of this international Institute. 

Professor Soltan attended as an official member of the Polish Delegation both 
the Ist in 1955 and the 2nd in 1958 International Conferences for the Peaceful 
Application of the Atomic Energy at Geneva. He attended also the International 
Conference on Nuclear Reactions in Amsterdam in 1956, the Rehovoth Conference 
on Nuclear Structure in 1958 and several others. 

The role which Professor Andrzej Soltan played in the history of Polish science 
is of a very great importance. The premature and quite unexpected loss of this 
eminent, very active and hard working man, exhibiting excellent health up to the 
last moment of his life, created an emptiness in Polish physics which it will be im- 


possibile to fill for many many years. 


Henryk Niewodniczanski 
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A. Promieniotwérczosé by Maria Curie-Sklodowska, PWN, Warszawa, 1953 (revised 
new edition by A. Soltan). 

B. Stany wzbudzone jader atomowych, Postepy Fizyki, t. III, zeszyt 2, 151—167 
(1952). 

C. Lecture Notes, Electrical Faculty, Lödz Institute of Technology: 

; a) Mechanics (107 pages). 

! b) Electricity and Magnetism (105 pages). —— 
D. Ernest O. Lawrence, Postepy Fizyki, t. X, zeszyt 2 (1959). 
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ZUR LOKALISIERUNG DER ENERGIE IN DER ALLGEMEINEN 
RELATIVITATSTHEORIE 


A. UHLMANN 
Theoretisch-Physikalisches Institut der Universitat, Jena 
(Eingegangen am 3. Juni 1959) 


Aufbauend auf Ergebnissen von Möller wird untersucht, in welchem Maße ein 
Energiebegriff, der sich auf eine affıne Tensordichte Ui mit th ae 0 stiitzt, eine Lokali- 


sierung gestattet. 


ik 


Ist M die Mannigfaltigkeit der Raum-Zeit-Punkte und erfüllt der auf M gegebene 
metrische Tensor g;, die Einsteinschen Gleichungen 


Rate, R = — x Ta 
so sollte der Begriff der Energie E nur von der metrischen Fundamentalform 
ds? = g, dx dx* und von einer (möglichen) Gesamtheit voneinander unabhängiger 
Beobachter abhängen. Es sollte deshalb E ein Funktional nur von ds? und einer 
raumartigen Punktmenge 2 aus M sein: 
E = E [ds2, Q]. -(1) 
Es ist dann E der Energieinhalt des in Q lokalisierten (Teiles des) physikalischen 


Systems. i ia j 
Wir können zur Vereinfachung annehmen, daß Q ein Stück einer raumartigen 


(3-dimensionalen) Hyperfläche ist. Man wird dann erwarten, daß E das Integral 
über eine Energiedichte ist, das über (2 erstreckt ist: 


B= f hav | (2) 


Versteht man hierbei unter dV das durch die Beschränkung der Metrik ds? auf Q 
auf dieser Hyperfläche definierte Volumenelement, so ist dV nur von ds? und Q 
abhängig. (1) ist deshalb genau dann erfüllt, wenn h nur von ds?, von (2 und von 
den auf Q liegenden Weltpunkten P abhängt: 
mes plas GP C Q (3) 
(133) 
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Z: 


Sei Q eine Hyperfläche und P, einer ihrer Punkte, Es gibt dann eine Umgebung 
U von P, und eine dort erklärte Funktion f mit 


PeQNU<>f(P)—f(P) =0, PEU (4a) 


und 


De | +0 für cr (4b) 


Genau dann ist bekanntlich 2 raumartig in den eben betrachteten Punkten, wenn 
grad f ein zeitartiges Vektorfeld ist. 

Wir können 2 somit durch Funktionen der Eigenschaft (4) charakterisieren. 
Damit die Funktionen f und f’ die gleiche Hyperfläche 2 charakterisieren, ist not- 
wendig und hinreichend, daß die Pfaffschen Formen df und df" gleichzeitig auf Q- 


verschwinden. Sind x!, x?, x? drei Funktionen, die f zu einem Koordinatensystem 


a ee ne 


{x} = (xl, 22, 43, x} mit «4 = f ergänzen, so heißt das 


dj = A df bi dx" b. dee) Bodas B DEO für PEQ. (5) 
Charakterisieren wir nun Q (stiickweise) durch ein Koordinatensystem {x'} und die 
Bedingung xt — x* (Py) = 0. Wenigstens für genügend kleine Stücke von Hyper- 
flächen können wir dann . i 


h= h [ds?, (X), Pl, | a*(P) — 24(P,) = 0 (6) 
setzen mit irgendeinem Punkt P, € Q. E Br i 
Wegen (3) mab dann gelten E50 
£ ; eae d 
| h [ds?, ey P| = h [ds?, (P ED d 
È f alls ER » 


hiis xs Jis 1133 Be i=! 0 für k= 123 4 l x* (P) = 
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Gegenüber diesen ist die Möllersche Energiedichte [1], [2] À invariant. Möllers 
Energiedichte ist aber nicht invariant gegenüber allen in (7) 
mationen. 

Ist z. B. Q durch x? = 0. gegeben, so ist Q auch durch y4 = 0 gegeben, wenn 
xt = s: y* und s = 0 auf 2 ist. Wir können nun x’ und y” (r = 1, 2, 3) so wählen, 
daß sowohl für die Komponenten g, des metrischen Tensors bezüglich {x°} als auch 
für seine Komponenten g;, bezüglich fy!) die Relation g; = gj, = 0 für i = 1, 2,3 
gilt. Betrachten wir mit den so konstruierten Koordinatensystemen die nach (7) 
zulässige Transformation 


zugelassenen Transfor- 


Wh == A* (xt, x 49 X), k= 12,3, vise sae. (9) 


Gl. (7b) ist erfüllt, weil x? auf Q konstant ist. 
Wegen g; = 0 für ¿= 4 gilt nach Möller [1], [2] 


2 
h = = AV— &a- (10) 


Hierbei ist A der Laplace-Beltrami-Operator, der zu der auf Q induzierten Metrik 
gehört. A hängt also nur von ds? und 2 ab. Beim Koordinatenwechsel (9) wird aber 


gua > S* gua = Baa (9a) 


für Punkte von Q und somit wird 
2 QAM T am 44 
h = — AV—gu >= — As Veh (9b 
x 


Offensichtlich kann s so gewählt werden, daß A # À' ist. Also ist in Möllers For- 


malismus 


Bern (dst@, Pa (biel ee A [ds?; 2] (11) 


4. 


Betrachten wir nun (11) im Hinblick auf den Möllerschen Eindeutigkeitssatz [2]. 
- Hiernach gibt es nur eine Größe TER die erstens eine affıne Tensordichte ist, zweitens 
die Gleichung On — 0 erfüllt und deren zugeordnete Energiedichte h (mindestens) 
unter der Gruppe (8) invariant ist. Diese durch obige Forderungen eindeutig a 
Energiedichte À ist nach [2] gerade durch (10) gegeben. Da die Forderung, daß “O; 
eine affıne Tensordichte ist, im Wesentlichen aus der Forderung nach Invarianz der 
Gl. TE, = 0 entspringt, erlaubt Formel (11) folgenden Schluß: , 

Es gibt keine affine Tensordichte, die nur von der Metrik abhängt und die es 
gestattet, die Energie als Funktional raumartiger Gebilde zu sp 

Die beiden Forderungen: “O ist affine Tensordichte und “Of , = 0 sind also 


mit der Forderung nach strenger Lokalisierbarkeit der Energie nicht vereinbar. 
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Diese SchluBfolgerung als auch die des Abschnittes 6 gelten allerdings nur unter 
der Annahme, daß ^C? höchstens Abteilungen zweiter Ordnung von g; enthält. Es 
erscheint jedoch unwahrscheinlich, daß sich durch Einführung höherer Ableitungen 


die Situation grundlegend ändern könnte. 


5. 


Der eben erwähnte Möllersche Eindeutigkeitssatz gestattet es abzuschätzen, bis 
zu welchem Grade eine Lokalisierung der Energie vorliegen kann, wenn eine affıne 
Tensordichte mit ©? , = 0 zugrunde gelegt wird. Zu diesem Zwecke betrachten 
wir die vorgegebene raumartige Hyperfläche Q als eingebettet in eine einparametri- 
sche Schar &(A) von Hyperflächen. Wir können annehmen, daB 2(A) für Werte von: 
À in der Nähe von A = 0 definiert ist und daß Q = 2/(0) ist. Eine solche Schar können 
wir beschreiben durch Gleichungen der Form f(P) = A wobei grad f ein zeitartiges 
Vektorfeld ist. Insbesondere können wir f als vierte Komponente eines Koordinaten- 


systems {x°} gewählt denken: 
PEZ) «—— x*(P)-—A. 


Ist {yf} ein anderes Koordinatensystem und soll y4 (P) = 4’ dieselbe Schar von 
Hyperflächen (jedoch möglicherweise mit einem anderen Parameter) beschreiben, 
so muf die Gl. (7b) nunmehr für alle Punkte der pen gelten. D.h.an 
Stelle von (7b) tritt die schärfere Forderung 


9y* A 
— —0 für k=1,2,3 und alle P. (13) 


Es muB also sein: 
yf = A(x4)s yt = A (xt, 92,03), k=1,2,3, (14) 


wobei wir noch g;4 = gj = 0 für i + 4 vorausgesetzt haben. 
Aus (14) folgt 


Möllers A geht deshalb wegen (10) in 


2 DA 9A 2 
= au (24) ou an 


x 


über, da 34/2x* auf jeder Hyperfläche der Schar X(À) konstant ist. Mit der Kon- 


stanten 
RAN sil 
Qut Porn % 
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transformiert sich also die Möllersche Energiedichte unter (14) gemäß 
h>ath (15a) 
und somit 


fr dV > a fr dV. (15b) 
Q Q 


Aus mathematischen Griinden erscheint es unwahrscheinlich, koordinatenunab- 
-hangige Bedingungen einzuführen, die bewirken, daß nur diejenigen Transforma- 
tionen (14) mit y* = x* zulässig sind. 

(15a) zeigt aber, daß wenigstens der Quotient A(P,)/h(P,) für zwei Punkte 
von Q koordinatenunabhängige Bedeutung besitzt. Dies kann man so aussprechen: 
Sei h Möllers Energiedichte, dann ist 


h(P) _ 
5 (P; = @ [ds?, Q, X (4), Py, Pa] 
mit P,, P, € Q und Q = E (0) (16) 


Hierbei ist X (A) eine Schar raumartiger Hyperflächen, die Q enthält. 


6. 


Obwohl unwahrscheinlich, wollen wir doch annehmen, wir könnten eine zu- 
sätzliche Bedingung « einführen, die die zugelassenen Koordinatentransformationen 
nicht nur auf (14) sondern sogar auf (8) beschrünkt. Unter solchen günstigen Um- 
ständen hätten wir an Stelle von (16) 


h(P) = @' [ds?, Q, E (A), z, P] 
mt P€ @, 2=-20 (17) 


zu schreiben. J - 
Wir bemerken, daß nach dem Vorhergehenden und nach Móllers Eindeutig- 


keitssatz durch (17) die maximale Möglichkeit für die Lokalisierbarkeit des Energie- 
begriffes unter den gemachten Voraussetzungen über “GË gegeben ist. Jedoch ist (17) 
und noch prägnanter die daraus folgende Abhängigkeit 

E = E[ds2, 2, Z (å), z], Q = 2 (0) (18) 
sehr fragwürdig. Tatsächlich müßte in die Messung der Energie durch auf Q pw 
findliche Beobachter die „Geschichte der Beobachter“ (und nicht die des Systems !) 
eingehen. D. h. die Energie würde nicht nur vom Bewegungszustand der 


bhängen. Dies ist aber im Rahmen der klassischen 


Beobachter bei der Messung a Š 7 
Physik schwerverstindlich. Will man also, diese Härten der Interpretation 


\ 
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umgehen, dann muB die Vivace daB Co eine affine Jen a mit “OF p= 0 


ist, fallengelassen werden. Eine kleine Chance besteht vielleicht auch in der Ein- 
führung höherer Ableitungen. 
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WAVE EQUATIONS INVARIANT UNDER DISCONTINUOUS 
GROUPS AND THE PROBLEM OF NUCLEAR FORCES 


BY P. G. O. FREUND* 


Chair of Physics of the Facultatea de Matematica si Fizica Timisoara 
(Received June 8, 1959) 


It is shown that physical space has only the role of a representation-space of 
a certain group. In a consequently discontinuous theory, space-time, therefore, may 
be left topological but the fundamental group is to be taken discontinuous. A 2-di- 
mensional theory of nuclear forces invariant only under a discontinuous rotation group 
is developed. Hence an angular periodicity of the nuclear forces results. 3 and 4-di- 
mensional generalizations lead to the non-euclidianity of space and space-time and 
thus may serve as physical supports to Einstein’s gravitation theory. It is shown that 
theories invariant under discontinuous groups don’t contain conservation laws. 


§ 1. About the deduction of wave equations from groupal invariance conditions 


The wave equations of quantum physics are completely determined by groupal 
invariance conditions. Thus e. g. the Klein-Gordon eq. is the most general homo- 
geneous, linear partial differential eq. of second order with constant coefficients! 
invariant under the Lorentz group (L. G.). Yet, all considered groups (in the’ case 
above mentioned the L. G.) are continuous. But in the case of a consequent quantum- 
theory such topological (continuity) pretentions are not justified, as they contradict 
the fundamental conception of quantum physics according to which the structure of 
universe is discontinuous. Though concerning atomic phenomena the topological 
approximation concords with experiment, its extrapolation to nuclear phenomena 
where discontinuous structure could be observed is not at all justified. W. Heisenberg 
(1936, 1938, 1939, 1939) tried to give an account of disconutinuous universe structure, 
by introducing an elementary length / ~ 10715 em. Yet a lattice model of universe 
with the periodicity constant l rises paradoxes (Schild 1948). In order to analyse 
whether such a method is indicated or not we shall reason as follows. In the wave eqs. 
there appear the spacial and temporal coordinates is ct. On p hand we should 
interprete them as notations for the physical quantities: space, time. On the other 
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1 Such eqs. we shall call in the following simply: linear while in the case of variabie coefficients we 
shall speak of semilinear es. 
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hand they appear as coordinates in a representation-space of a linear representation 
of the group under which the wave eq. is invariant. 

At the construction of the formalism attached to the respective wave eq. latter 
interpretation is necessary and sufficient. Insisting yet upon a unification of the two 
interpretations we are led to the idea that physical space-time is but a representation 
space of a group. In order to make clear what is to be understood by this, let us 
consider an example. In the transformation 


x =xcosy—ysing 


y = xsin p + y cos 9 (1.1) 


x, y have no more meaning than that to be the components of a vector in a space in 
which there is defined a representation of the group of plane rotations (P. R. G.) 
Concerning the representation (1.1) it should be observed that this representation is 
irreducible in the field of real numbers as there doesn't exist any real one-dimensional 
subspace of the plane invariant under P. R. G. Thus if x, y appear in a differential eq. 
to be constructed of invariance conditions under (1.1) it is sufficient to remain at 
this interpretation for x, y. 

If the group is continuous its representation-space will be topological too. But 
a representation-space of a group is representation-space of any of its subgroups 
(under the induced representation). Thus the same topological space may be the 
representation-space also of a discontinuous subgroup of the initial continuous group. ` 
Thus it is not necessary to pass to a discontinuous physical space not even in a con- 
sequently discontinüous theory (i. e. a theory invariant under a discontinuous group) 
as the role of physical space is limited to that of a representation-space of the group 
which, as we have just seen, may be topological even for a discontinuous group. 

Thus a consequently discontinuous theory can start from a differential eq. (and 
not e. g. from an eq. with finite differences) as this means to claim the continuity of 
X, y, z 1. e. the continuity of the representation-space which is, without any physical 
restrictions as we have seen, admissible. Concerning the discontinuous group, we 
shall choose it as a subgroup of the continuous group of the approximative theory 
and the quantization of the group will lead us to a universal constant (not necessary 
with physical dimensions e. g. an angle). The discontinuos group evidently is not 
allowed to possess infinitesimal transformations. 

A theory covariant under a discontinuous group will be called in the following 
Q-covariant under that group. E. g. a theory covariant under the discontinuous L. G. 
(discontinuous P. R. G.) will be called Q. L. covariant (Q. P. R. — covariant). 


$ 2. Q-covariant formulation of the mesonic theory of nuclear forces 


As it is well known the scalar mesonic theory of Yukawa is L. G.-covariant. 
The consideration of a static (time independent) problem in this theory must therefore 
start from an eq. invariant under the static subgroup of the L. G. i. e. the space rota- 
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tion group (S. R. G.) Such an equation was proposed by H. Yukawa: 
A, G — u? G = — g ó(r) (2.1) 


where: 
22 22 22 


Ag eu 9x? | Oy? | 922 


r=\x +y? + 22 (2.2) 


From this eq the well known potential with short ray of action of Yukawa: 


(Qu 
G —-- 


dues 


(2.3) 
is obtained. 


For not complicating too much our further considerations we shall limit them 


to the equivalent plane problem i. e. the problem invariant under P. R. G. when 
instead of eq. (2.1) we have: 


4, G — u? G = — g d(r) (2.4) 
where: 
22 9? 
ze T yz 


r= y x? -—— y (2.5) 


Also in this case a potential with short ray of action is obtained. 
& 
G= £ h (lal D) (2.6) 


(for the proof see Appendix I). 
Before passing to the Q. P. R.-invariant eq. let us observe that also the eq.’ 
As c- go got) (2.7) 


invariant under the continuous S. R. G., leads to potentials with short ray of action. 
In fact the singular solution of this eq. is (see Appendix II): 


N rE] tex, 


j-1\* 
v n 2.8 
x X (n deep imt Kad 
k=0 


Eq. (2.7) differs of eq. (2.1) by being semilinear while (2.1) is linear. 
XE 
2 | being a fundamental length of the theory (like the Compton wavelength u”! in (2.1)) 
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Looking for the Q. P. R.-invariant eq. we start from the plane equivalent of eq. 
Qu 


]-1 


A G IC) (2.9) 


which as (2.7), also admits solutions with short ray of action. 


We shall quantize the P. R. G. as given by (1.1) to: 


x’ = x cos ngo — y sin ng 


y = x sin np, + y cos ngo (2.10) 
Po = 2a/v 
v = universal constant 


where x, y are coordinates in the same topological representation-space as x, y were 
in (1.1). This is admissible as we have seen in § 1. 

For passing to a Q. P. R.-invariant eq. it is sufficient that for the coefficient of G 
in the first member of eq. (2.9) instead of — /~1/r, which is a P. R. G.-invariant function, 
there should be chosen a function which should be Q. P. R. invariant without being 
also invariant under the continuous P. R. G. (evidently any function invariant under 
the continuous P. R. G. is also Q. P. R.-invariant, but not also conversely so that 
such a choice is possible). In order to choose this function we shall reason as follows. 
The function — /-!/r is up to the constant — I”! equal to the reciprocal of the function 


z = r= yx? + y? representing in the space x, y, z a cone, the generatrix of which 
makes an angle of 45? with the plane xOy (fig. 1). This surface is evidently invariant 
under the P. R. G. of the xOy plane. A surface Q. P. R.-invariant, yet not also 
invariant under the continuous P. R. G. of the xOy plane is the regular v-faced 
pyramid? inscribed in the cone of Fig. 1, » being the constant of (2.10). 

The eq. of this pyramid is: 


sin (m +1) = — sin m= 


iA peg Ais Tae 


a OT 
sin — 
v 
(2.11) 
cos m 2 — cos (m + dae 
d EE Rue ums 
sin — 
y 
(m = 0, 1,2, ... ... ,»— 1) 


E 3 It should be remarked that this pyramid is invariant also under a lot of symetries with respect to 
planes passing through the axis of the pyramid. These superflue symetries could be eliminated choosing 
instead of the pyramid a surface like that in fig. 2. 
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Now the Q. P. R.-invariant eq. will be: 


[-1 
ANG ee OT) (2.12) 


with z given by (2.11). 
| Evidently the singular solution of this eq. is a periodical function of y with the 
period gy (of course it also possesses a short ray of action). Thus the potential of 


Fig. 1 Fig. 2 


nuclear forces and thus nuclear forces themselves are periodical functions of Q. 
This fact could be verified experimentally; a determination of the constant v could 
be realised. Of course pọ is very small and consequently » very large. 


3. Observations 


1. The angle 9, is neglectible with atomie phenomena, nevertheless it ought to 
be observed with nuclear phenomena. 

2. Our method of deduction of the eqs. (2.9), (2.12) includes a rather arbitrary 
step by the choice of the function z (2.11). Yet the arbitrariness of this step is not 
greater than that is done when linearity is claimed from eq. (2.1). A better choice of 
the function z must start of experimental data. 

3. We treated a simplified 2.dimensional model. Really we ought to have started 
from a 3-dimensional corresponding to the S. R. G. (eq. 2.7). In this case as O: S. RG 
would have appeared the group of a regular polyhedron. But it is well known that 
there exist but 5 groups bacause of the fact that the sum of the angles of a triangle 
must be 180°. These five polyhedra corespond to very high values of pọ what is in- 
admissible. Yet this is not an impediment to the generalisation of our method because ` 
of the following. Regular polyhedra with an arbitrary high number of faces can be 

constructed in a non-euclidean space where the sum of the angles of a triangle is 
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high enough (of course, > 180° i. e. an elliptical space). This effect seems to be very 
likely, as by this Einstein’s gravitation theory could be interpreted as the theory of 
global curvature of space due to the statistical behaviour of many individual curva- 
tures due to Q-covariance of nuclear phenomena. Of course generalization to 4-dimen- 
sional space-time and other transformation properties of the wave function (e. g. 
spinorial) comport but mathematical difficulties. 

4. As conservation laws in theories covariant under continuous groups are due 
to this covariance by means of the Noether-theorem, where the existence of infini- 
tesimal transformations is essential, in a Q-covariant theory conservation laws are 
excluded. This means that in nuclear phenomena conservation laws are not anymore 
valuable. Such a result would also offer a physical justification to Hoyle’s (1948, 1949) 
macroscopical conceptions. 

5. It should be observed that the intrinsic structure of physical space is not 
essential in the theory, its significance being limited to that of a representation space. 
Thus in physics not space is fundamental, but the group i. e. the wave equations, 
i. e. the particles. This idea is related to some of the philosophical ideas recently ex- 
pressed by W. Heisenberg (1958). 

The author wishes to thank Prof. J. Rayski for a critical remark and Prof. 
B. A. Rosenfeld for an interesting correspondence concerning the geometrical pro- 
blem discussed in the third observation. 


Appendix I. 


Proof of relation (2.6) 
The singular solution of eq. (2.4) is:! 


dus sd hes. " 
rds kia (1.1) 


Choosing the unity vector k, in the direction of 7: 


kr = kr cos | (1.2) 
and 
d? k = dk, dk, = k dk dy (1.3) 
By (L1), (L2), (1.3) we obtain: | 
oo 2n 
& kdk ikr cos : 
ca-£ |; E 95 9 dy (1.4) 
0 0 


But (see Watson (1945) § 2.2: (5)): 


2n 
f er cos do, — 2z J, (kr) (L.5) 
0 


. 1 The solution of the homogeneous eq. is put as usually equal zero. 
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By (1.4), (1.5) we have: 


G() = £ I og (L6) 


(O28 J Kb gu 
0 


There exists the relation (see Watson (1945) 8 13.6: (2)) 


eo 


k+! J, (kr) Pw y | 
la Lapu 2 (A + 1) Ky—a (Iul r) (1.7) 


0 


Putting v = A = 0 into (1.7): 


k Jo (k 
IE e dk = ko (|u| r) (1.8) 
0 
By (L6) and (1.8) finally: 
G(r) = 5 ko (ll n) (1.9) 


(1.9) is just (2.6). 


Appendix II. 


Proof of relation (2.8). 

In order to find the singular solution of eq. (2.7) (the solution of the homogeneous 
eq. is put as usually — 0) we shall appeal to the following method usual with semi- 
linear eqs. Let be D a semilinear, homogeneous differential operator and let us consider 


the eq.: I 
| ve (11.1) 
Be g, the eingenfunctions of the eq.: 

Dg, = Yn En (11.2) 


coresponding to the eingenvalues y, If the g, form a complete orthonormalized 


system then 


O(r) = >; an En (11.3) 


where: 


a = gy (0) 


where the asterisk indicates complex conjugation. 


146 P. G. Oliver Freund 


By (II.1) and (II.3) we obtain: 


CG) = -daia oora sS Dat (> a, £,) (15 - 


Considering (11.2), we can write (IL5): 


G0) — —gD $a = Lg y a DD 


pvp P CPI IPTE 


Are i 
3 
Finally: 
G= z ut (11.6) 
In the case of eq. (2.7): 
| E la 3 
D = A, — 3 ; (0.7) 


In this case (see Sommerfeld (1951) ch. II, 81: (17), (4a), (2): — 


; noir Er 


Smet condition with suos to à and Gi adt di 0f 3-575 al 
i ae ae a Tied AP SOWA GR. BEER; ow ir gie: £ys120g ii 


abis: ; vg sel bus Joie IPTE ipli pns - anh ER NS us 


(a) 8, (r) = No? Lr ORS os 2 Pr (cos 0) ene E 3 
S (b) I| <! m | (13 — 
4 | LL 4 a 


p- 
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There is (see Sommerfeld (1951), ch. II 82: (25), (25a), (20)): 


NE ooo S (n—1)! 
L, (0) = do In (0) = (n —1)! 1! [> EHRT —k— 1)! (k + 1)! k! a) 


as all terms of the sum for k >n vanish. 


Thus: 


n—1 


oS 2 o* 
Ind = (md? 5 (Dr (n—k— 1) 4- DERI ue) 


k=0 


It should be mentioned that some authors write (II.12) with a supplementary factor —1 
(Schiff (1955) : (16.22)). 
By (11.10), (II.11), (II.12) we have: 


n— 


1 
Lud? aaa o* 
En (r) Vans 2nle EA 1) G-ECDREX DIE (11.13) 
Hence: 
0h 
ga (0) = & (0) = às nh (11.14) 
mt 


(11.3), (11.4), (11.13), (II. 14). yields: 


co n—1 
EIER = : oF 
ölr) = D: ay nine => ee ties (11.15) 
n=1 = 
But (see Sommerfeld (1951) ch. IL 1 : (10)) 


s EN (11.16) 
di 2n 


? 


By (11.6), (IL. 8c), (11.14), (11-16) we obtain finally: 


oo 


1 LE (n—1)! 
GO) = z s= >. ° ee eS 
n=1 
111% 
x NO n IL17 
ue a n—k—11G Dk | qun 
k=0 


o see that the series in (II.17) converges for any value of 


That is just (2.8) It’s easy t 
r= 0. The exponential factor indicates the short ray of action of the potential. 


148 | P. G. Oliver Freund 


REFERENCES 


Heisenberg W., Zs. f. Physik 101, 533 (1936). 

Heisenberg W., Zs. f. Physik 110, 251 (1938). 

Heisenberg W., Ann. d. Physik 32, 20 (1939). 

Heisenberg W., Zs. f. Physik 113, 61 (1939). 

Heisenberg W., Naturwiss. 45, 227 (1958). : 
Hoyle F, Monthly Not. Roy. Astron. Soc. 108, 372 (1948). — 
Hoyle F., Monthly Not. Roy. Astron. Soc. 109, 365 (1949). 

Schiff L. I, Quantum Mechanics 24 ed. Mac Graw Hill N. Y. 1955. 
Schild A., Phys. Rev. 73, 414 (1948). 

Sommerfeld A., Atombau und Spektrallinien II. Band Vieweg & Sohn Braunschweig 1951. 
Watson G. N., Theory of Bessel functions 24 ed. N. Y. 1945. 


re er are AAAS 


| o pru | 
T lij F: —— v p . = M J £ tea La h: ^ 


E a 

O uo elie IE 

è ' LN p 
tw 20013 eae: "2 x í 

E he "mu ah 

ys = P Tamia E 


=i > 


Vol. XIX (1960) ACTA PHYSICA POLONICA Fasc? 2 


A MOLECULAR THEORY OF LIGHT SCATTERING IN GASES 
AND LIQUIDS 


By S. KIELICH 


Institute of Physics, Polish Academy of Sciences, Poznan 


(Received June 15, 1959) 


The general principles of a statistical-molecular theory of light scattering in an 
isotropic medium consisting of polar-anisotropic molecules are given. 
A general, fundamental equation for the scattered intensity I is derived, containing 
the molecular factors Fj, and F,,;, which account for isotropic and anisotropic light 
scattering as brought about by the molecules of the medium. The molecular factors Fj, 
and F,,;, are discussed in detail for the case of gases and liquids possessing molecules 
whose linear dimensions are small with respect to the light wave-length. In such gases 
and liquids, F;, depends on the mean polarizability of the molecule and on the radial 
intermolecular correlations, whereas F nis depends on the symmetry and anisotropy of 
polarizability of the molecules and on the orientational intermolecular correlations. 
From the fundamental equation giving the scattered intensity, general expressions for 
the optical anisotropy 4? and degree of depolarisation D of the scattered light, and for 
` Rayleigh's ratio S and the extinction coefficient h are derived. For ideal gases, these 
expressions ,reduce to the known formulas of Rayleigh-Born-Cabannes. For liquids 
the molecules of which exhibit the axial symmetry, the expressions obtained for D, S 
and À differ from those of Cabanues-King-Rocard in the anisotropic term by the 
angular correlation factor Romp which also appears in the formulas for the Cotton-Mouton 
and Kerr constants. For compressed gases the molecular factors Fj, and Fynis are 
expanded in inverse powers of the volume of the system. The coefficients of the expan- 
sion, ae Vim and Be penis, termed, respectively, the first and second virial coeffi- 
cients of isotropic and anisotropic light scattering, are computed for spherical molecules 
of variable polarizability and for axially symmetrical anisotropic molecules possessing 
a permanent dipole moment. Moreover, following Buckingham’s method, the effect of 
the internal molecular field and the hyperpolarizability of the molecules on light scattering 
in liquids is accounted for. Finally, general relations between the quantities D. S and 
h and formulas relating these quantities and the anisotropic term in Kerr's 
constant Kanis are derived. The relations thus obtained contain no molecular parame- 
ters and may serve for checking the theory by experimental results. 


1. Introduction 


The theory of classical scattering, in which the wave-length of the light is pre- 
served, is due to Lord Rayleigh (1899), who considered a gas of spherical, non-absor- 
bing molecules of linear dimensions small with respect to the wave-length. 

Smoluchowski (1908) was the first, to prove that light scattering in optically 
homogeneous media arises from the spontaneous, thermal Spree of their gen 
Einstein (1910), who developed Smoluchowski’s ideas, proposed a phenomenologica 

i (149) 
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i. e. thermodynamical theory of scattering in liquids whose molecules are spherical”. 
Gans (1923) generalized the Smoluchowski-Einstein theory for anisotropic molecules. 

The ,,molecular” or electronic theory of the scattering of light in gases and 
liquids possessing optically anisotropic molecules is due to Born (1917, 1933), Ca- 
bannes (1921, 1929), King (1923), Rocard (1928) and others. 

Investigations by the foregoing authors proved that light scattering by a medium 
consisting, quite generally, of anisotropic molecules is given by the superposition of 
the scattering on fluctuations of its density and of the scattering occurring on fluctu- 
ations of the optical anisotropy of the molecules. Moreover, Rocard (1928) proved 
that scattering, especially in the neighbourhood of the critical point, is affected by the 
fluctuations of the intermolecular field. The foregoing theories assumed molecular 
interaction between the molecules, spherical or anisotropic, to depend solely on their 
mutual distances. Clearly this was a weak point of these theories, since both in real 
gases and in liquids consisting of anisotropic or polar-anisotropic molecules the 
interaction between the latter is strongly dependent on their mutual orientation. 
Hence, in addition to radial correlations, a theory of light scattermg by condensed 
media should also account for orientational intermolecular correlations. 

Müller (1936) was the first to consider orientational-molecular interaction, using 
the method of Fowler-Debye, in his theory of scattering in non-polar liquids. Anselm 
(1947), Benoit and Stockmayer (1956), and Prins and Prins (1956) proposed theories 
of the degree of depolarisation of light scattered by real gases and liquids, wherein 
the orientational-molecular interaction is described by Zernike and Prins’ (1927) 
correlation function, as generalized to account for orientation (cf. also Volkensteyn, 
1951, and Weill, 1958). Buckingham and Stephen (1957) proposed a general theory 
of the degree of depolarisation of light scattered by a dense medium; their theory 
accounts i. a. for scattering on assemblies of spherical molecules and for the effect 
of hyperpolarizability of the polar molecules. 

In a former paper (Kielich, 1958a), the author proposed a generalization of 
Born’s ,,gas” theory (1933) to liquids; therein, in addition to light scattering on 
fluctuations of the density and of the anisotropy of polarizability of the molecules, 
account is taken of scattering arising from the orientational interaction of molecules 
of arbitrary symmetry. In the present paper a general, statistical-molecular theory of 
light scattering in an isotropic medium consisting of anisotropic-polar molecules of 
arbitrary symmetry is proposed. 


2. Fundamental Statistical-Molecular Theory of Light Scattering in an Isotropic 
Medium 


. We shall be considering an isotropic medium of volume V containing a large 
number N of identical, optically anisotropic, polar molecules. The origin 0 of the 
fixed system of reference (X,) is located at the centre of the volume. Let molecular 


*) A molecular theory of light scattering by optically isotropic molecules of a fluid is given 
by Fixman (1955). j 
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systems of reference (x(?), p = 1, 2, ... N be rigidly attached to each molecule of the 
medium. The position and orientation of the molecules in volume V are given by their 
configurational variables v = + (r, œ); here, r = (ri, v,, ... Py} are those determining 
the position, whereas œ = (0, Wy, ... 0 n} determine the orientation of the molecule. 

The probability for a mutual configuration of the molecules in an element dt of 
the configurational space is 


dP (z) = f(t) dv, 


with f(t) being the configurational statistical distribution function. By classical statis- 
tical mechanics, we assume f(t) to be given by Gibbs’ canonical distribution function 


f(t) = Ce Wen (2.1) 


wherein C is the normalisation constant, k — Boltzmann's constant, T — the Kelvin 
temperature of the system and U(r) — its total potential energy. 

Electric field of scattered light. Consider a parallel monochromatic light beam 
of electric field 


E(t) E e à)" 


incident on the volume V of the medium. E denotes the electric field amplitude, 
y — the frequency, and t — the time. 

Thus, the electric field of the incident light beam at the centre of the p-th molecule 
of the medium situated at the distance t, from the origin 0 of the system (X) is of 


the form 
E(t)” EA a 1 (2.2) 
wherein 
2n 27v 
k= "t So = ia So; 


with c and A denoting the velocity and wave-length of the incident light, respectively ; 
So is the unit vector in the direction of propagation. 
The electric field given by (2.2) induces the time dependent dipole moment 


m (n)? = m’ eat (2nvt — k: rp) x (2.3) 


within the p-th molecule of the volume V. Here, m? is the amplitude of the dipole 
moment of the p-th molecule. By classical electrodynamics, the dipole induced in the 
molecule radiates electromagnetic waves in all directions. In the present paper, only 
Rayleigh scattering, i. e. scattered light of the same wavelength À as that of the in- 


cident light will be considered. 


At considerable distances from t 
light wave scattered by the p-th molecule of the medium is of the form (cf. Born, 1933) 


he scattering system, the electric field of the 


Ag? 


i(2nvt — k. r. ze 
E (jo — — ggg (R, x [R, x mn ra ae 
p 
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wherein R, = Ry — r, is the vector directed from the centre of the p-th molecule 
to the point at which the radiation is observed. The latter point is at the distance Ry 
from the origin 0 of the system (X,). The distance at which the point of observation 
is situated is great with respect to the linear dimensions l of the scattering volume 
and with respect to r,; thus, the following expansion may be written: 


" 
R, = |R; — r| = Ro — s? : T, + [rj p er = ot, - (2.5) 


wherein S? is the unit vector in the direction of the observation vector R, = s! Ro, 
thus determining the direction of propagation of the scattered wave. 

If Ry > l> r, then R, may be replaced by Ry everywhere in eq. (2.4) with the 
exception of the exponential factor, wherein, with respect to the expansion (2.5), the 
substitution is R, = Ry — S° ` r,; this yields 


; An?  -i2n (or =) : ; 
E'()9 = — ye a” ls? x (s? x mh? (P) ei@k—k') rp)}, (2.6) 
š 
with 
poi ene 
A c 


denoting the wave vector of the scattered light. 
By eq. (2.6), the total electric field resulting from the light wave scattered by all 
the molecules in volume V at the point of observation is given by 


R Age? f ; 
a 2 Fa) = — aR X (s' x Meam), (2.7) 
wherein 
N ; e 
M = » MP) RR): rp (2.8) 
p-1 


denotes the amplitude of the total dipole moment induced in the medium of volume V 
by the electric field of the incident light wave of amplitude E. 

Intensity of scattered light. We shall now consider the total intensity of the 
light scattered by the volume V, in order to determine the component thereof trans- 
` mitted by a Nicol prism analyzer at the point of observation. Let the direction of 
the vibration transmitted by the nicol be that determined by the unit vector n per- 
pendicular to the vector of observation Ry = R, S’. Since now $°: n = 0, the com- 
ponent of the electric field E" in the direction of the unit vector m is given by: 

An? 


E n= PR, M nei, (2.9) 
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The intensity of the light scattered by the volume V of an isotropic medium 
and transmitted by a nicol is defined as follows: 


I, = CE n (EB? -n)*5g ; (2.10) 
herein, the symbol < g stands for the statistical mean value in the presence 
of the electric field E of the incident light beam, and the asterisk denotes a complex 
conjugate. Let  — € (z, E) denote an arbitrary state function of the system in the 
presence of the electric field E; then, with Gibbs' distribution (2.1), the statistical 
mean value of Ø is determined at thermodynamical equilibrium of the system as 


follows: 
U(t, E) 


che ft, Eje * ddr... dry: (2.11) 


here, U (z, E) is the total potential energy of the system under consideration in the 
presence of the electric field E of the incident light wave. 

Substituting (2.9) in (2.10), the fundamental equation for I, may be written 
as follows in tensor notation: 


ee "(MyM 2.12 
n R 5 c noa» (2. ) 


the summation indices o, v assuming the values 1, 2, 3. 

At not too great intensities of the incident light, the dipole moment M induced 
within the volume V may be assumed to be a linear function of the electric field 
E. The dipole moment may then be written as follows: ; 


M, (v, E) = Agr E, = As e, E, (2.13) 
wherein 
2 M, 
Ar = (ee (2.14) 


is the tensor of the optical polarizability of the medium of volume V, and e is the 
unit vector in the direction of the electric field, E = e E. p 
By (2.13), eq. (2.12) yields for a spherical scattering sample of volume V in vacuum 


4 2 2 
ala 22) aem em 
0 ‘ 


um) 
n? +2 
of refractive index n; the symb 
mean value in the absence of an electric field (E — 0): 
UG) 


p= C Í if yt f d()e * dw dr, ... dt: | (2.16) 


Since the unit vectors m and € satisfy the relations 
Ey 


: .` 
wherein I, = [ EE’ is the incident light intensity within the spherical sample 


ol < > without lower indice denotes the statistical 


é, n = cos Os sree, le ls 


` 
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averageing over all possible directions of m and e with respect to the fixed system 
of reference (X,) yields 


Nee, = = ((8 cos? Q, — 1) (0,,0,, TÓ, 9, ,) + 2(2 — cos? Q,,)0,,0,.), (2.17) 


with Q,, denoting the angle subtended by the unit vectors e and m, and 6,, being 
Kronecker’s symbol: 


pom I foro 7 
Wie We es ee M 
The optical polarizability tensor 4,, being Hermitian, 4,, = A,, and with 
respect to (2.17), eq. (2.15) yields the following general, fundamental equation for 
the intensity of the light scattered by an isotropic medium: 


I, = Q {5 cos? Q2,,F., + (3 + cos? Qon) Fanio (2.18) 
wherein 
Re (2.19) 
d us x 3 Ç (36, Ov é 3 Ögr 650) Aas A, = (2.20) 
with 
(m on M [n2 Pru 
Vespa e Q2 


The quantities F;, and F nis will be termed the molecular factors of isotropic and 
anisotropic light scattering, respectively. With (2.8) and (2.14), Fis and F nis may by 
rewritten in general form involving molecular parameters 


am?) Qm», š 
Fis (s) = < Oye 2E» e puts. e» £ (2.22) 
pa d e 


Im’) | Im’ 


l v —-is:r 
he. (s) == ° (80, ss = On: Ove) Y JED FEE) e 2 (2.23) 
pa Š 


wherein the quantity 
S:T,, — (k — k): Ty (2.24) 


determines the difference in phase between the light scattered by the p-th and q-th 
molecules of the medium, with r,, = r,— r, being the vector connecting their centres, 
and 


s=k—k’ —<(s,— 8). (2.25) 


The polarizability tensor of a molecule immersed in the medium. As a first 
approximation, it may be assumed that within a moderately condensed medium the di- 


aee cq > 
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pole moment is induced in the molecule only by the electric field E of the incident light 
wave. If the medium is a condensed one, however, an internal field F® due to the 
charge distribution of all the other N—1 molecules and superimposing its effect on 
that of the field E(? appears and acts on the p-th molecule. The dipole moment in- 
duced in a molecule immersed within the condensed medium is given by (see Bucking- 


ham and Stephen 1957) 


Sp Sp) (p) (p) (p) p 1 
ma? = as (EP + FP) + UBER (EP + FP) (EP + FP) + a yap (BQ + 


Je (E DS EEUU FI) (2.26) 


where a? is the polarizability tensor, and Be up are the hyperpolarizability tensors 
describing the field dependence of the polarizability. These tensors are symmetric 
in all suffixes and have been discussed by Buckingham and Pople (1955 a). 

The field Ff? may be computed by Kirkwood's method (1936). Namely, if 
the medium consists of anisotropic-polar molecules, and considering only dipolar 
interactions, we have 


N 
122] = — I TO (ng 4 mg) > (2.27) 
= 

wherein ue is the permanent electric dipole moment of the isolated molecule, and the 
tensor of dipolar interactions is given as follows: 


1 Da 
T4 = — (pa Sap — 3 Toga Toa); (2.28) 
"pq ^ 

the latter assumes non-zero values for p Z q and is zero for p — q. | 
From eqs. (2.26) and (2.27), the following formula is obtained for the field po 

at the centre of the p-th molecule of the medium: 
puo TS Fo Tel > T oe Hs ar >= m u Tim at Ee i. F 
q qr 

_ >) Tp agp TH? a) TS af) EP + (2.29) 


ars 


(s) T(GD 4G pO. /! 
DTG ag rp TP ah Tip eg BP — 
qrst 
with F® denoting the intermolecular field.existing at the centre of the p-th molecule 
gs y L Li .` LI ` 
immersed in the medium when the electric field E of the incident light wave is ab- 


sent (E = 0). The field FG is given quite generally as follows: 
N 
FQ = — TS m (2.30) 
q=1 ; 


wherein 
Fm = WP raf FE + (23 
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is the total dipole moment of the molecule immersed in the medium at E = 0, or 
explicitly : 
(pa) (Q) (pa) e Ca) Tar) „(N 
ne > T us + > Top py lys Mi 


< y BA a B m THO ju pu (2.23) 


qrs 
By (2.26), the polarizability tensor of the p-th molecule immersed in the con- 
densed medium is given as follows: 


dm? Ol soy) () ple) p) eg 
— = Meere as Boas FY dec Yoaßy I $5 ar . Cay i: vro di 
IEP J E-o : OE; 
or, with respect to the expansion (2.29), explicitly: 
Im P) : 
[255]. "e X 0 rne neo 
X = q 
ik - (pa) (Q) ) 4G) PCs) „(9 pik: rps | 
+ Dal TEP og TID al erm — yat TY? ag Ti ap TO ag erred 
qr qrs 
= > we De a Tt a 208 a T a) e ik-rpt 7. (2.34) 


qrst 


The foregoing expansion shows the polarizability tensor of a molecule immersed 
in a condensed medium to be constant in the first approximation only; generally, 
however, this tensor is a function of the intermolecular field. 

Potential energy of molecular interaction within the medium. From the defini- 
tion of eq. (2.16), in order to compute effectively the molecular factors F, (s) and F, (5) 
given by (2.22) and (2.23), respectively, the potential energy U (v) of the mutual 
interaction existing between the N molecules immersed in the medium of volume V 
should be known. The energy U (z) is dependent, in general, on the structure of the 
molecules and on the thermodynamical state of the medium, and, in principle, con- 
sists of two terms 


U (7) = U® (r) + UO (rw), (2.35) 


U! (r) denoting the potential energy inherent in radial interaction of the molecules 
only, and U? (r, œw) denoting that due to radial-orientational interaction. 
For U(r), the Lennard-Jones potential (1924) is usually applied (cf. koples 1954), 


of the form í 
LM YS N 
U(r) = 4e* (2) — (=) i (2.36) 


wherein the first term represents the repulsive energies and the second — the attrac- 
tive. The quantity e* represents the negative value of the lowest potential energy, 
r is the intermolecular distance and r, — the value of r at which the attractive and 
repulsive energies exactly balance. 
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It is customary to include in U(r, w) the energy of the Keesom dipole-dipole 
interaction, as well as those of the Debye-Falkenhagen permanent dipole-induced 
dipole interaction, the London anisotropic-dispersive interaction and, finally, the 
multipole molecular interactions. For anisotropic-dipolar molecules, by (2.32), 
Barker’s expression (1953) in tensor notation is obtained for U(r, œ): 


Uo (r, ©) = 1 2 u? 22 um Res 1 > p> D e p 2 no zs 
par 
tt Doe Tap? cj Ty) am Ta pa — i (2.37) 
herein, the consecutive terms determine interaction within assemblies of two (pq), 
three (pqr), four (pqrs), ... of the N molecules present within the volume V of the 
medium. 

The fundamental equation (2.18) containing the factors Fj, (s) and F,, (s) as 
determined by (2.22) and (2.23), together with the polarizability tensor of the mole- 
cule as given by (2.33) and the potential energy of molecular interaction (2.35) jointly 
determine the general form of the statistical-molecular theory of light scattering by 
an isotropic medium of an arbitrary nature. 


3. Light Scattering in Gases and Liquids. 


Factors F,, and Fg, for molecules exhibiting constant polarizability. If the 
medium is but a moderately condensed one, the effect of the intermolecular field on the 
polarizability tensor of the molecule may be neglected in the expansion of eq. (2.34), 
and the expressions (2.22) and (2.23) may be written in the first approximation as 


follows: 
F;G) = (0, 0,, ] a9) aQ" e*t "aay, (3.1) 
m 
p (s) = i qo an om = Ose 650) D a) a" e is: va) à (3.2) 
pa 


For further discussion it is convenient to refer the tensor a? (and, similarly, xD) 
to the molecular system (x); -1,2,3 attached to the p-th (or,. respectively, q-th) 


i ` 
molecule of the medium. Thus, we have the following transformation formulas: 


PAG ADS 

4p =o of a, a = of of aff ; (8.3) 
2 : ; ul (p) 

in the case under consideration, the fixed (X );= 1,2,3 and molecular (x;’); —1,2,3 

the transformation coefficients o, ... eX? have the 


systems being rectangular ones, | 
meaning of the cosines of the angles subtended by the axes of these systems. The 


cosines satisfy the general relationship 


ó, for p= q; "` i 
| IL COR Od 86. 3.4 
Ôi, wP o = o5 Ook Mm for p# 9 2 
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w$? denoting the cosine of the angle subtended by the i-th and k-th axes of the systems 
(x) and (x(?) attached to the p-th and q-th molecules, respectively. 

The tensor a?) determines the constant polarizability of the isolated molecule, 
1. e. the N, independent of intermolecular fields. In optically inactive sub- 
stances and throughout spectral ranges widely remote from those in which absorp- 
tion appears, the tensor ate) is symmetrical and real. 

If all the scattering molecules present within the volume V are of one kind 
aff? = a = a, then, by (3.3) and (3.4), eqs. (3.1) and (3.2) may be rewritten as 
follows: 


(au A eee V oe 
Fis(s) = (asas 2 s spa (3.5) 
1 sinsr 
Lai (s) = — (aan ) Bog off? — dy du) 3.6 
| DER. 2 " (3.6) 
wherein 
27 At .. 9 | 
Fe ig (3.7) 


9 denotes the angle between the direction of propagation Sọ of the incident wave and 
that of the scattered wave 8’, and is termed the angle of scattering. 


Isotropic light scattering. With the radial correlation function g (r) introduced. 


into the theory of X-ray scattering by crim and Prins (1927), e (3. >) id 


Fis (se = Z; snl ge der sin sr ul m jak] » š 


here g o) Az r? dr is the probability for the presence of a molecule | within as 
volume element 4a r? dr at a distance r from the. molecule whose position | is fixed. 


If the linear dimensions of the scatterin nz molecules are small as compared tothe — 
Zu. ote D ‘i SE Mp OP 


light wave-length . A, we aye. ZA E b EE a ; 
sin sr $r 4 


Ss "apna han) ear m ee Kurt 1, f 


TNT. Se ee RR cg e 


~ Ge 
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denoting the mean polarizability of the isolated molecule. It is seen that F; depends 
on the mean polarizability of the molecules « and on a quantity determining the 
radial intermolecular correlations (see, Zernike and Prins 1927): 


N 
yR = l+ is | [ev = x r? dr = x kTBr, (3.12) 
0 
wherein 
Te URN ar ji 
D A IDA DAT E 


is the isothermal compressibility coefficient of the medium. 

Smoluchowski (1908) and Einstein (1910) showed that the quantity yg accounts 
for the light scattering resulting from spontaneous density fluctuations in the medium. 
As both @ and vg are isotropic quantities, Fis determines isotropic scattering only. 


Anisotropic light scattering. If the scattering molecules are small in terms of the 
light wave-length, then, with respect to (3.9), eq. (3.6) yields 


pn i Xij Cy Ç > (369 ox? dies m y) = (3.14) 
pa 


It is seen that the factor Fynis 
mined by the polarizability tensor %;; and on the angular intermolecular correlations 


determined in general by the statistical distribution of the directions of thé molecular 


depends on the anisotropy of the molecules as deter- 


axes. 
If the axes of the molecular systems coincide with the principal axes of the respec- 


tive molecules, eq. (3.14) may be rewritten in the form 
Per; a; < (3 cos? Oo — 1) (3.15) 
ij pa 


with 9g denoting the angle between the i-th and j-th principal axes of the p-th and 


` q-th molecules, respectively, and a; — the polarizability tensor referred to the prin- 


cipal axes of the molecule. j | 
If the molecules possess the axial symmetry, the axis of symmetry being the 


3-axis, then & = % Z %» and (3.15) yields 
Fas = 3 (zs — %)®< UG cos? goose (3.16) 
pa : 


OCD denoting the angle between the axes of symmetry of the p-th and q-th molecules, 


and «4,0, — the polarizabilities of the molecule in the two principal directions 


parallel and perpendicular to its axis of symmetry. — 
If, in particular, the molecules exhibit the spherical symmetry, & = 9 and, 
by eq (gno Fo = 0. This means that F nis determines the anisotropic light scat- 
2 3 ? ` anis í 
tering only. The latter arises chiefly from the anisotropy of the molecules and from 


. ¥ . 
the angular intermolecular correlations. 
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Light scattering in ideal gases. In ideal gases, the molecules are mutually inde- 
pendent and may assume all possible configurations within the volume V. Since now 
exp |- | = 1, the statistical mean value < > of eq. (2.16) reduces to 
the isotropical average. As shown by Smoluchowski (1908), the density fluctuations 
in an ideal gas yield yp = 1, so that eq. (3.10) assumes the form 

gas P3, = 9 a? N. (3.17) 


as determined by eq. (3.14) may now be expanded as follows: 


The factor Faris 
Fais = $ @;; G N (3603 0; — 050, + 2i (Bog? of = 0,0,)); (3.18) 


gas” anis 
qp 
with respect to 
o? of? = $ ôy Ôm for p Z q, (3.19) 
we finally obtain 
ee x= i (3 a Qj — Xii "T N. (3.20) 


In ideal gases, the anisotropic scattering depends only on the anisotropy of the 
isolated molecules, whereas isotropic scattering depends on their mean polarizability. 


Factor accounting for angular intermolecular correlations. Eq. (3.16) may be 
written as follows: 
anis — (a — aj)? NRom; (3.21) 
the quantity x 
Ro SE x x (3 cos? Oa) — 5 — (8.22) 
q 
appearing therein is termed angular intermolecular correlation factor, if the mole- 
cules possess the axial symmetry. A factor Roy of the same form also appears in the 
theory of molecular orientational effects (cf. Buckingham 1955, and Piekara and 
Kielich 1957, 1958). For an ideal gas U(t) = 0, and eq. (3.22) yields: 
Rey = 1+ d», B-cos? O°? — 1) = 1. (3.23) 
qp i 
` The correlation factor (3.22) accounts for the angular interaction of the molecules 
entering molecular assemblies which coalesce and disappear continually within the 
medium; such assemblies present a momentary axis of maximum polarizability aris- 
ing statistically from the effect of ordering of the molecular axes. From the foregoing 
formulas it is seen that in an assembly of mutually ,,coupled“ molecules light scat- ` 
tering differs from that resulting from free molecules, for which Roy = 1. 
For instance, for the case of molecules entering aggregates of two, Rom as first © 
computed by Piekara (1939, 1950) is given by 


| ii 
Rom —2—37, (3.24) 
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wherein L = L(y) is the well-known Langevin function and y = 


m the energy 


of dipole coupling in kT units. This refers to nitrobenzene, for which Rojas calcu- 
lated according to eq. (3.14) for the Cotton-Mouton effect (hence the indices CM) 
fits well with the experimental data. It is of interest that Rom for y#0 is always 
larger than unity. 

Anselm (1947), and Benoit and Stockmayer (1956) introduced the correlation 
function g (r, o) into the theory of the degree of depolarisation of scattered light; 
this function has the properties of the radial distribution function g(r) of Zernike 
and Prins (1927) as generalized for orientational correlations of the molecules. With 
the function g(r, œ), the correlation factor (3.22) takes the form 


Rom =1+22 [Je cos? @ — 1) E w) — M r?drdo, (3.25) 


r œ 


wherein @ is the angle subtended by the axis of symmetry of a given molecule and 
that of an arbitrary one at distance r. j 


4. Optical Anisotropy and Degree of Depolarisation of Scattered Light 
Optical anisotropy. The ratio of F,,, and Fis which determine, respectively, 


anisotropic and isotropic light scattering, i. e. 


Fani | 4.1 
A B x 
ee (4.1) 


yields a measure of the optical anisotropy of the scattering medium. 
Substituting herein the expressions of eqs. (3.10) and (3.14), we have 


a; au CQ Got wf? — di; Sus)? 


A? zm ducali alils: Mi E ee VT p (4.2) 
2a; & N L aF in f E (r) i $i 4o ar) 
0 


this is an expression yielding a measure of the optical anisotropy of a liquid consisting 
of molecules of arbitrary symmetry. 

Hence, for ideal gases, we have aA 
of an isolated molecule 


2 — 62, wherein ó, is the optical anisotropy 


3 = 3a; Qij — Ri Gi : (4.3) 
: 2a; 0j 


or, if the polarizability tensor be referred to the principal axes of the molecule, 


2l (z — ag*- (0 = ta)? + (05 — t? ! (4.4) 
Oa = 2 (a, + % + ts)” 


162 S. Kielich 


For molecules having the axial symmetry the foregoing expression reduces to 


2 2 
2 Hy — 01 % — &ı 
== — = 5 4.5 
de (2 3^ ze) | 3a $ 
clearly, for spherically symmetrical molecules ô, = 0. 


In condensed gases and liquids whose molecules have the axial symmetry, eq. 
(4.1) yields, with respect to (3.10) and (3.21): 


A? = ó Rom ; (4.6) 
YR 


this means that the optical anisotropy of condensed gases and liquids depends on 
that of the isolated molecule, 6,, and on the ratio of the orientational (Rom) and 
radial (y,) correlations. | 

Introducing Roy from (3.25) into eq. (4.6), we obtain the Benoit-Stockmayer 
(1956) formula: 


r- afl f foro nre ead e 


the latter has since been generalized by Benoit and Weill (1956) een by Weill (1958) 
for dilute solutions of liquids. 


EN ratio of scattered light. in a medium Bun anet AAC 


EEN ieee T 


scattering (F nis Æ 0), the light scattered undergoes depolarisation. The degree of ` 


depolarisation is usually defined as the ratio of the lowest (Znin) and highest (L nax) 
possible values of the scattered intensity (see, Born 1933): ` 


E — i o ERN 
Trin. and Jax ‚being ES from the fi ere atio 
Ifthe RM beam incident. on the ae im dium 
light scattered is observed with. a ed m we ha 
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When the nicol is rotated, the scattered intensity transmitted attains its maxi 
value with the former transmitting vibrations parallel to the plane of the ME 
light peels E only, i. e. at @, = 0. The minimum scattered intensity is observed when 
the nicol is set so as to transmit vibrations perpendicular to the plane of E, i. e. at 


> -> 
Fig. 1. Space arrangement of the unit vectors s and s in the direction of the incident and scattered light 


> > 
beams, and of unit electric vectors e and n in either light beam, respectively. 


„ = 90°. Thus, by the definition of eq. (4.8) and with respect to (4.11), a general 
expression for the degree of depolarisation of light scattered by an isotropic medium 


is obtained as follows: 
5 cos? @ Fig (s) + (6 + cos? 9) Fanis (5) (4.12) 


Di) = 5 Fi (s) + 7 Fuss (s) 


When observation of the light scattered is carried out at an angle of 90° with 
respect to the direction of incidence and if the scattering medium satisfies the condi- 
tion of eq. (3.9), the following relation is obtained for d = 90^: 


eX. 6 Fanis 
D()mD-rp TE (4.13) 


The degree of depolarisation D thus defined is related to 


3 Fanis 
D = E 4.14 
5 Fis + 4 Fanis ) 


plane polarized incident light, by the formula 


the degree of depolarisation for 


2D’ 
i = ——. 4.15 
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The formula (4.13) (or (4.14)) together with expressions (2.22) and (2.23) is 
identical with the one derived by Buckingham and Stephen (1957) in their theory of 
the degree of depolarisation of light scattered by a dense medium. 

Degree of depolarisation in some special cases. Substituting Fis and Finis» as 
given by (3.10) and (3.14), in the general expression (4.13), an expression for the degree 
of depolarisation of light in condensed gases and liquids possessing molecules of 
arbitrary symmetry is obtained (Kielich 1958a, Piekara and Kielich 1958): 


65; an (Bw? wif — ó; Ôr) > : 
De IN — s L0. Ns 
lO aj; z; Nyg + 7 a5 oer C21 (3 of? cff? — à; ôr) Y 
pq 


For substances whose molecules exhibit the axial symmetry, (4.16) reduces to 


692 Ro 


ON cH (4.17) 
9yR--7 ôz Rom 


where ó, = (x4 — x.) | 3a, and Roy is given by (3.22). If the angular correlation 
factor Roy is given by (3.25), the eq. (4.17) is identical with the result obtained by 
Anselm (1947), and Benoit and Stockmayer (1956). 


For an ideal gas, (4.16) yields the Gans-Cabannes formula: 


er; 6 (3a; wiz — c; A) 6 ba 


wo Das; T ms ay are) SHE (4.18) 


For a substance whose molecules possess the spherical symmetry, 6, = 0, and 
the degree of depolarisation assumes the value zero, which means that in such a me- 
dium light undergoes no depolarization. 


5. Rayleigh's Ratio and Coefficient of Extinction 


Rayleigh’s ratio. With an incident beam of natural light, and if observation is 
carried out without a nicol, eq. (4.11) should, additionally, be averaged over all 


values of g, from 0 to 2x, thus 
^ 2x 


: 1 
LO) t> il I (Gn, 9) do, ; (5.1) 
0 
hence, with respect to (4.11), 
I (9) = 3 Q {5 (1 + cos? 9) Fis (s) + (13 + cos? 0) Fanis (s) }- (5.2) 
The scattering constant, also known as Rayleigh’s ratio, is defined as follows: 
I(9) Ro 
5(8) = 1€) Bo ; : (5.3) 
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o 


by (5.2), the following general expression for Rayleigh’s ratio is obtained 


I2» (in? +2 V? 
se = ap |) | 3 ) (60.6 costó) Fuld) + (18 + e050) Fa): (5.4) 


If observation is at an angle of 90° with respect to the direction of incidence 
(9 = 90°) and if eq. (3.9) is satisfied, the foregoing expression reduces to 


Ae la 1 (2m\* [n3 + 2\? 
sees | > ) (5 Fis + 13 Fais). (5.5) 


Extinction coefficient. In passing through the medium, the intensity of the 
light beam diminishes as a result of scattering (the effect of absorption is neglected 
here). On emerging from a layer of thickness J, the intensity is given by the 
well-know Lambert law 


IL=Ihe™, (5.6) 


wherein Ah is a coefficient dependent on the properties of the medium and on the light 
wave-length A, and determines the weakening of a beam of intensity J, resulting from 
scattering. By definition (see, Cabannes 1929), 


2n x 
h = f [ S (0) sind dd dg, (5.7) 
00 


i. e. by (5.4), a general expression for the coefficient of extinction (turbidity) is hence 


derived: 


X [2% Sin? 42 "t 2 
ebd | 3 js {5 (1 + cos 9) Fis (s) + (13 + 
0 


+ cos? #) Fanis (s)} sin ddd, (5.8) 
wherein F, (s) and F4, (s) are given by eqs. (2.22) and (2.23). 


In particular, for small molecules, i. e. if the condition of eq. (3.9) is satisfied, 
the foregoing equation yields 


axial 9m AÉ d % 
S em MS M (Pic 2 Fai): 5.9 
e £a 3 (Fis + ) 62 


S and h for special cases. Substituting F; and Fanis from eqs. (3.10) and (3.14) 
in eqs. (5.5) and (5.9), the following expression is obtained for real gases and liquids 


consisting of molecules of arbitrary symmetry: 


4 2 NY. 
si (2:3) N yaa |t [1 | (ep) | + 


0 


fs ign < (3 o£? wf? — 9; >. (5.10) 
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8x [2n\* (n? +2 V N i N 4 
h = ds B = N auo su| 1+ 4 | [s 3. "a| - B= 
E 


1 | 3 
bs + 2, (3 o£? wif? — db, tu) >t. (5.11) = 


If the tensor of optical polarizability be referred to the principal axes of the j 
molecule, these equations yield: | 


1 [22|* [+2 PN < | 
s= i (*| | 3 ) ee et on ae TON “gh (3 cost Of — 1), >), 


¿j= 1 1 
(5.12) 
3 f 
8z [2=)% [n2 + 2)? N « 1 "MN 
h= 97 EJ | 3 y 2. & Xj CES Wi 2 (3 cos? OF -5|. 
(5.13) 9 
. With the Lorentz-Lorenz formula 
nà—l 4a N Š E 
pig ee (5.14) 


and for molecules having the axial symmetry, eqs. (5.12) and (5.13) reduce to 


2(52 _N2 
E (ret 9 em 
Me 3 (n2 — 1)2 | DER 
o lawn) m 
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For an ideal gas consisting of spherical molecules the Born-Cabannes formulas 
reduce to the well-known formulas given by Rayleigh: 


, 22 (n— DV 3223 (n — 1)? V 
gas) jx» DIY N , gash T 324 N o (5.20) 


6. Second Virial Coefficients for Light Scattering in Compressed Gases 


General form of virial coefficients. Virial coefficients have been computed by 
Harris and Alder (1953) and by Buckingham and Pople (1955) in the theory of electric 
polarisation, and by Buckingham (1955) in that of Kerr’s effect. The virial coefficients 
for light scattering in compressed gases will now be computed by the method of the 
latter authors. For this, the molecular factors Fi, (s) and Fanis (S) given by eqs (2.22) 
and (2.23) are expanded in powers of 1/V: 


T Bs (> 
Bo a er (6.1) 


Bs. Cer 


Fanis (s) = AS" + 3: 


ARE (6.2) 


V y 
the coefficients Ü pee pss pnis. isd cais + are termed the first, . second, 
third, ... virial coefficient of isotropic (F,,) and anisotropic (Fa) light scattering, 
respectively. 


The first virial coefficients are defined as follows: 


is — lim Fi, (5), Zu = lim Fu (5); (6.3) 


V> V-oo 


by (2.22) and (2.23) we have - 
Ai 9a? IN, Amit 902 ON, (6.4) 


wherein & is the mean polarizability of the isolated molecule, and 6, denotes its ani- 
sotropy as given in general form by eq. (4.3). Comparison with eqs. (3.17) and (3.20) 
shows that the first virial coefficients A5 and Ag“ account for light scattering by 


& ” 
free molecules (“gas molecules"). 


By definition, the second virial coefficients are given by 


Bis = lim ([F,,() — As] V), (6.5) 
V= 
; B5“ = lim Fa) — AS] V}; (6.6) 


V> 
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by (2.22), (2.23) and (6.4), they may be written in the following general form: 


Im P? Im,® ij —is'r 2 EN 
ar ett 9s N >Vt, (67) 


Au. 
BE — lim. 14 85,0, 
> Pq 


V oo 


mi?» dU Cen 
BP = lim C; 5 (895. Öre — dor NA SEO E i e it "pa — Oa? 62 Nr]. 
(6.8) 


As the second virial coefficients contain only contributions accounting for the inter- 
action of molecules in pairs, eqs. (6.7) and (6.8) may be rewritten as follows 


A 2 *(1) m.» 2(2)\* ah 
s SO | | lan pum) fe Frame (9) 
panis 2x N? Im (am IY 

EU Xs | (e in dar Pi 9 pO si JEL 


Ui; 
— 18a? 62 ES r? dr do, (6.10) 


wherein m», m, 2 are the dipole moments of molecules 1 and 2, Uj, is the potential 
energy of interaction between a pair of molecules, and Q = f dw is the integral over 
all angular coordinates. u. 

Spherically symmetric molecules. For simplification, assume the scattering 


molecules to have spherical symmetry. Then, by (6.4), the first virial coefficients are 
AS = 952 NX. ASEO, (6.11) 
i. e., the first virial coefficient of anisotropic scattering vanishes. 

To compute the second virial coefficients, the expansion (2.34), which accounts 
for the effect of the intermolecular field on the polarizability of a molecule immersed 
in a condensed medium, will be applied. Since for spherically symmetrical molecules 
the tensor &,, reduces to the isotropic tensor a; = œ Ó;, and with respect to ex- 
pansion (2.34) and condition u): eqs. (6.9) and (6. 10) yield 

T at ate 
Bs = 36 2a? N2 ee ar. tube RS r3 dr, (6312) 
: a? a3 at _ U(r) 
Bs = 727 a? N? — $247 = + «Je kT 72 dr. (6.13) 
r r r 
r. 


Pople (1954), assuming U(r) in the form of the Lennard-Jones potential (2.36), 
proved that 


Un r omn j ; 
free e AT 2 dr = m y H, (y), (6.14) 
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wherein 
27—n © 6m "n 3 
H, = 6 Om A] ym 
o = y De | x ) = 6.15) 
A g* % l 
with y = 2 le The functions H, (y) have been discussed and their values tabu- 


lated by Buckingham and Pople (1955b). 
With eq. (6.14), we obtain finally from (6.12) and (6.13) 


XE 24m æt N? a i 
Boc 007 NV RE = —=— = 2 IM (ES HS yas ACA VIE 2 (6.16) 
ro y To To 
anis 6r a* N? & : 
I a {Ha (9) +2“ Ho) + 7% Ha 0) + Jem 
TQ y TQ TQ | 


wherein yg is defined in analogy with (3.12). 

Thus, in compressed gases consisting of spherical molecules, Au: Ba) 
and 427s — 0, whilst B$" > 0. This means that free spherical molecules give rise 
to isotropic light scattering in a medium. However, in compressed gases, intermole- 
cular forces cause the spherical molecules to form assemblies, which represent aniso- 
tropic scattering centres for the incident light beam. It is seem from eq. (6.17) that even 
the presence of assemblies consisting of no more that two spherical molecules is suf- 
ficient for the medium to exhibit additional, anisotropic scattering. Hence, on passing 
through such a medium, light may undergu depolarization. 

Axial symmetrical molecules with dipole moment. If the anisotropic molecules 
| possess the axial symmetry, then, neglecting all terms in the higher powers of the tensor 
æ; in expansion (2.34), we may write, by (6.9) and (6.10) 


2 A Un 
ps = d ff kT 72 dr dw, (6.18) 


252 172 US 
pu — Den 1 ih (3 cos? O02 — 1)e AT r? dr do, (6.19) 
wherein 942 is the angle subtended by the axes of symmetry of molecules 1 and 2. 

Combining eqs. (2.36) and (2.37), the following expression is obtained for the 
potential energy of interaction of two dipole molecules (the permanent dipole-induced 


dipole interaction is neglected): 


Uy, = 4e* die e =) = = (3 cos OW cos OP — cos 6?) (6.20) 
r 


r r 


with cos 912 = cos 0? cos 9 + sin GC? sin QD cos (pg — Pı), Wherein the angles 
Q02, @%, ... g, are explained in Fig. 2. 
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Substituting (6.20) in (6.18) and. (6.19) we obtain, by (6.14), 


3 
za? rox? N? wayt 


is 3 š 
Bs = 9a? NV (yn = 1) “Ten Be eel IM, (y) +0 Ba (y) + a $ (6.21) 


D 


4% 
; 


Xe ER > at 
Fig. 2. Angles determining the mutual orientation of two dipoles (D) and Z@) in the fixed system of 
reference, (mayo thie et 
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7. Effect of the Molecular Field and Hyperpolarizability of the Molecules on Light 


Scattering in Liquids 


In a condensed medium, light scattering depends also on the molecular field 
and on the hyperpolarizability of the molecules (see Buckingham and Stephen 1957). 
In the present case, when computing the molecular factors F; and F,,;,, the tensor 
of the polarizability of the molecules, as given by the expansion (2.33), should be 


applied. The expressions of (2.22) and (2.23) may be written as follows: 


Fis) = bor de V] (a a + VER) eos s m, (7.1) 
24 


Fais (5) = 5 ; (85. Öze — dar Öse) >: (a ase + Verve) 008 B: TS (7.2) 


wherein the tensor 


US P a?) pa mn ze CER pa mue ae pa pa Ee F9 as 


vey OTY ory M voô 
(p ) ) 
Lt hale yD FO FD + am, PO IY + o 


(p) (a) P) gpa 
+ (a? oe al ps Hy Toe > (7.3) 


oa Gp Eo Ope + Sar 5 JEO + 5 FG) 2 E@ 


determines the effect of the internal molecular field and hyperpolarizability of the 
molecules on light scattering in a condensed medium. 


Light scattering on assemblies of spherical molecules. In order to obtain insight 
into the part played by the internal molecular field in the scattering of light, it will be 
assumed for the sake of simplicity that the condensed medium consists of assemblies 
of spherical molecules. Thus, m; = Bj; = 0 and aus Yu reduce to the isotropic 
tensors given by 

Xi; = M Ygkd — 3 $ y(Ó;; Oe + 0; Ó + Ôi m (7.4) 


and the foregoing expressions yield, quite generally : 


: (b) (a) Eo 2 Fe 
= > | 2 Fa 3 9Fg ; 
Fis (s) = AF |“ 9+3 3h epe Li T 9E 


+ Say QR + Fp FA) + A Bd 1) 


2 po 9 po 9 Fo 2 Fe S US 16 
Fuis() = > OS t Ë EP OEY EP TA Su Las 22 s 
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If 2 > r, the foregoing relations yield a general expression for the optical aniso- 
tropy 4? of a condensed medium consisting of spherical molecules: 


Qe IF? IFP IFP ary > 
IEP IED IE? OEY 
At = : ee, (7.7) 


OPPO NORD Or! UR 
2 9--3—*- 13 4 
Sl Sn E 


Hence, it is clearly seen that the optical anisotropy of such a medium arises 
from fluctuations of the molecular field. If the internal molecular field F is here 
assumed to be given by the expansion (2.29), which, in the case of spherical molecules, 
takes the form 


FP — F® — a Y) THEO + a? Y) TOP TS) ES — ..,, (7.8) 


then eq. (7.7) reduces to 


eX TET) 
At = o No UN 7.9) 
6Nyg + 202 >) (TY? T FIRE ( 


pars 


This yields, by approximation, 


At un (X TO TOY = 4n a° V 


-4 
er =: ET Br r£ (r) dr: (7.10) 


As Fus Z 0, the light scattered on assemblies of spherical molecules undergeos 
depolarisation (D = 0). Buckingham and Stephen (1957) showed that, in the case 
under consideration, the degree of depolarisation D may be expressed by the difference 
between the molecular refraction of the fluid (Ry) and that of the perfect gas (Rio. 

Polar liquids. Buckingham and Stephen, in the theory proposed by them (1957), 
proved the value of the degree of depolarization D in polar liquids to be strongly 
dependent on the hyperpolarizability of the molecules. 

If an isolated molecule possesses an axis of symmetry, and if this is the 3-axis, 


then the following components of the teusurs 5, @; and ß,;, have non-zero values: 


Hag Un Os Oey S= Be F Os, = 03, Bug = Boos = Bis Z Bs = Bss- (7.11) 
As a first approximation, expressions (7.1) and (7.2) may be reduced to the form 


Fal) = (Y (es +20) + (Bos + FE, (7.19) 
= ` STpq 


Feais(s) I E <> {læs — 5) + (Bas — Bis) Fos}? (8 cos? G0 — 1) HS : 
(7.13) 
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"mmt and Stephen (1957) showed that, in the case of Onsager's model 


po26-1o0-—10D 5 
t Soin ni ee 


(7.14) 


wherein e is the static dielectric constant of the medium. Thus, for small scattering 
molecules, eqs. (7.12) and (7.13) reduce to 


2 (e — 1) (n? — iB |° 
En = 9a? hi 4 ( 3 ct n?) ) E Nyr; (7.15) 
an n? — 2 
Fans = 9a? h + Uu m N 62 NRcm> (7.16) 
1 s 
P= (at Px, a= mim. (17) 


with f denoting the mean hyperpolarizability of the isolated dipolar molecule, and 6, 
yielding a measure of the anisotropy of its hyperpolarizability. 
Substituting (7.15) and (7.16) with Roy = 1 in eq. (4.13), the Buckingham-Ste- 


phen formula for the degree of depolarization in polar liquids is obtained: 


4 (e — 1) (n? — 1) uB 58] ia 
E fı di uus (25 + n?) a? bq da 


5 nero . (7.18) 
4 (e — 1) (2 — 1) HV ye fr 1) (n? —1) rel 


She 
sp 


BE ne 3(2e+n) aq 


Similarly, by (5.5), the following formula for Rayleigh’s ratio is obtained: 


(ns Zi = 4 (e —1) (n? —1 e YR + 


Site 214 N 3Qe+n*) a? 


— 4(e—1) (n? — 1) HB 56] 2 d 
AE mee |} Sol 


These are formulas which, if the quantities D and S are known from experimental 
data, yield the value of the hyperpolarizability of the isolated molecule. 


8. General Relations 


From the theory, as discussed in detail in §§ 3 to 7, it is seen that, in order to 
D, S and h for a given condensed substance, the quantities 


compute the quantities , 
; and angular intermolecular correlations should be known, 


accounting for the radial 
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in addition to those describing the properties of the isolated molecules. Since, for 
a given substance, it may well be the case that not all the molecular parameters are 
known simultaneously, it is convenient to compute the quantities D, S and A from 
relations containing only such quantities as are accessible to experimental determi- 


nation. 


Relations between D, S and h. Eqs. (4.1) and (4.12) yield 


5 (D (8) — cos? 91 
6—7 D(®) + cos? 0 ° 


F, anis (s) 11/9 Em. 
OL A? (9) = (8.1) 
a relation whose left hand side coincides with the theoretical definition of the optical 
anisotropy of the medium, and whose right hand side makes it possible to compute 
A? direetly from experimental data on the degree of depolarisation of the light scat- 
tered. 

By (8.1), eq. (5.4) determining S (9) may be rewritten in the form 


_ 1 (2n [n2 + 2)? (6 + 6D (0) sin? |, 
56 = a (2) | 3 ee, ca 


1 [2x [n2 + 2)? {6 + 6D (8)) sin? d 
(2) fee) Penis (s) , (8.3) 


249) xr ie 3 | 'S(D(8) — cos? 8j 
wherein F;, (s) and F,,, (s) are given by (2.22) and (2.23), and 9 is the angle of scat- 
tering. 
If observation is perpendicular to the incident beam ( = 90°), eqs. (8.1) — 
(8.3) yield, for scattering molecules of small dimensions, 


ee Š 
Fis 4 6 — 7D (8.4) 
1 [2mM [n? + 2? 6 + 6D i 
s= gp |) | 3 ) $59 n. eo 


reme (8.6) 


Under identical conditions, by eq. (5.9) and with respect to (8.4), the extinction 
coefficient is obtained as follows: 


_ 8a 2n\* [n2 + 2 2643D 

ne) | 3 o pen 
_ 8a [22|* (n? + 21* 6-- 3D | 
m) | 3 N ra ao 
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Finally, the foregoing expressions yield the general relations: 


i 
D(0) = D (a+ = cos? o), (8.9) 
2 fot Messi 
S (29) = S (1 + 1T D eosta) , (8.10) 
_ 8x 2+D 
ho ants T (8.11) 


which hold for a medium of arbitrary density consisting of scattering molecules of 
small dimensions. Analogous relations were derived by Cabannes (1929) for light 
scattering on isolated molecules. 

If the intermolecular field existing within the condensed medium be neglected, 
the molecular factor F, of isotropic scattering may be expressed by experimentally 
known quantities. Namely, with respect to the Lorentz-Lorenz equation (5.14), eq. 


(8.10) yields 


F. 81 p ==) 


2 
= ! 12 
167? aa Maren S 


Substituting the foregoing expressions in eqs. (8.5) and (8.7), we obtain the 
Cabannes-King-Rocard formulas : 


__ x2 (n? — 1)? 6 + 6D 


S 224 Sp t Br. 
(8.13) 
2 _ 12 
iak 823 (n 1)? 6+3D kT Br, 


344 6—7D 
which contain no molecular parameters and may be checked experimentally. 
Anisotropic light scattering and Kerr effect. The molecular factor F nis account- 
ing for anisotropic scattering may be expressed by the anisotropic term in Kerr’s 
constant. By a general relation, 
1215 n? (n? — 1) kTV K.. 
Fanis = 57, (s — 1) (e + 2) (n? + 2)? i 
wherein & denotes the electric permittivity of the medium. The anisotropic term in 
Kerr’s constant is, in general, given by the following expression (cf. Kielich 1958b): 


(P) Q (q 
z (n? + 2)? (e + 2)* < TERRY Arte S Tie 8.15 
Ka — gone Ga (30s, Oxe oT o) > 9E, Cp ( ) 


(8.14) 


with m® denoting the dipole moment induced in the molecule by the variable electric 
qum Dose uM. and m® — the total dipole moment of the q-th molecule 


v a ` . . 
in the presence of the static electric field Ey inducing optical birefringenve in the 


y 


` medium. j 
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For molecules possessing constant polarizability, eq. (8.15) yields 


_ a(n? + 2)? (e + 2)? aij A (pq) (pa) _ s. s 
Konis = 077916 52 LTE > (30; ik Ol Oy Oxi) 7 (8.16) 


wherein a,; is the electric polarizability tensor of the isolated molecule- In particular, 
if the oleas have the axial symmetry, eq. (8.16) reduces to 


n? —1)(e— 1) (n? -2) (e +2)V 


> 
Kanis = 120z n? kTN 


HS RR (8.17) 


with Roy denoting the correlation factor given by (3.22). 
For the case of a condensed medium consisting of spherical molecules, the ge- 
neral expression (8.15) yields 


z (n? + 2? (e + 2? x a IFP IFO IFP Fp N 
Qu P 


1215 n? kTV a 


Kanis = OE, OE, OF, SE 


or, as a first approximation: 


z (n? + 2)? (e + 2)? a? a? de 
Kan = T ERATES CTP TS >. on 
pars 


The quantities D, S, h, and Kerr's constant K,,;,. Substituting eq. (8.14) in (8.6) 
and (8.8), we obtain the general expressions 


7223 n? (n? — 1) kT 


S= 6-1 63 wy 0 DD Ke rv 
19224 n? (n? — 1) kT (OD RE ` (8.21) 


suits zT) (e + 2) (n? + 2) 


which contain experimentally accessible quantities only, and hold for condensed 
media consisting of molecules of arbitrary symmetry. 
Combining expression (8.4) and eqs. (8.12) and (8.14), we have 


A = A a AN oe, T Kanis ë 

6—7D (n®?—1)(—l(®+2 (+2) Br’ (8.22) 
hence, Gans’ formula (1923) is obtained: 

Kage et Dr I (8.23) 


ann 6— 1D’ 


which yields a.relation between Kerr’s constant and the degree of depolarisation of 
the light scattered. 
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Substituting F,, and F,,;, as given by (8.12) and (8.14) in eqs. (5.5) and (5.9), 
we have 


2(n2 — 1)? 3127 n? Kan; 
S — T (n zi Zn anis 
224 Fiat (n3 — 1) (e — 1) (n? + 2) (e + 2) J' SER 
j ale 2407 n? Kanis 


—1)2 

une near €" 

Conclusions. The quantities D, S and h accounting for light scattering in an iso- 
tropic medium are given (88 4—7) in terms of various molecular parameters inac- 
cessible to experiment; on the other hand, the same quantities are related mutually 
and with Kerr's constant, thus yielding the basis for comparing the theoretically 
predicted results and those obtained experimentally. Thus, a study of light scattering 
may serve not only to determine the shape and properties of the isolated molecules, 
but also to obtain valuable data on the nature of the molecular interaction existing 
in a condensed medium. l l 

The author wishes to express his sincere thanks to Professor Dr A. Piekara for 
his valuable advice and for his helpful interest throughout the present investigation. 
The author is greatly indebted to Dr A. D. Buckingham for his stimulating discus- 
sions of some of the problems under consideration. 
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Beitrag zur Ausarbeitung einer neuen 


Vol. XIX (1960) ACTA PHYSICA POLONICA Fase? 


ZUR FESTSTELLUNGSMOGLICHKEIT EINER PLASTISCHEN 
DEFORMATION MIT HILFE EINES GEIGER-SPITZENZAHLERS 


von J. MADER unD B. Sujak 


Institut für Experimentallphysik der Universität, Wroclaw. 
Physikalisches Institut der Polnischen A. d. W., Wroclaw. 


(Eingegangen am 18. Juni 1959) 


Es wird der Einfluss einer inneren Deformation auf die photostimulierte Koexo- 
elektronenemission von Aluminium-Proben untersucht. Die innere, durch Einpressung 
von Stahlstäben erreichte Deformation lässt sich wegen der langsameren Abklingung 
der Emissionsintensität in der Druckpunktumgebung feststellen. Es werden zwei Deu- 
tungsmöglichkeiten diskutiert: Eine durch Herstellung von Emissionszentren, die aus 
dem Innern der deformierten Probe nach der Oberfläche diffundieren, die zweite durch 
eine mit der Deformation verbundenen Kohäsionsänderung der Oberflächenschicht. 
Im weiteren werden Messergebnisse über die Feststellungsmöglichkeit von Druckeingriffen 
in Polimer-Proben bei Ausnützung des Einflusses der bei Druckeingriffen entstehenden 
Elektrisierungszuständen auf die Arbeitsbedingungen des Zählers, angeführt. 


Kramer (1950) hat auf das Ansprechen eines offenen Geiger-Spitzenzählers auf 
sogenannte „Exoelektronen“ (Bericht über die Diskussionstagung, 1957) hinge- 
wiesen, welche von deformierten Oberflächen einiger Metalle während oder nach der 
Deformation emittiert werden. Einige Jahre später wurde von einem der Verfasser 
(Sujak 1957a, b) gezeigt, das der Geiger-Zähler auch zu Untersuchungen von defor- 
mierten Polimer-Oberflächen sich eigne indem der Einfluss der im Deformations- 
bereich entstehenden Elektrisierung auf die Arbeitsbedingungen des Zählers aus- 
genutzt wird. 

In der vorgelegten Arbeit werden weitere, das erwähnte Problem betreffende 
Ergebnisse angeführt. Das erste Ergebnis betrifft den Einfluss einer sogenannten 
inneren Deformation von Aluminium auf die photostimulierte Koexoelektronenemis- 
sion, das zweite eine Erzeugung von elektrisierten Gebieten beim Deformieren von 
Plexiglas. In beiden Fallen erlaubt die durch den Geiger Zähler, in der Zeiteinheit 
registrierte Impulszahl eine Feststellung wie auch relatives Vergleichen von den vor- 
her plastisch deformierten Bereichen. Die angefiihrten Ergebnisse sollen einen weiteren 

Untersuchungsmethode von plastischen Defor- 


mationen in Metallen und in aus organischen Stoffen hergestellten, für das weisse 


Licht undurchsichtigen’ Elementen bilden. 
(179) 
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Einfluss einer inneren Deformation auf die photostimulierte Koexoelektronen- 
emission aus der Oberfläche des technischen Aluminiums 


Die photostimulierte Exoelektronenemission aus Metalloberflächen nach voran- 
gehender, im weitesten Sinne des Wortes gemeinten Bearbeitung (mechanische Ab- 
schabung, plastische Deformation, Beleuchtung mit elektromagnetischer und korpus- 
kulärer Strahlung, Temperung usw.) erregt Interesse wegen hervortretenden Erklä- 
rungsmöglichkeiten einiger, die Festkörper-Oberfläche betreffender Probleme, wie 
auch wegen Aussichten auf technische Anwendung. Das angehäufte experimentelle 
Material erlaubt einstweilen keine eindeutige Klärung des Emissionsmechanismus, 
obwohl er in groben Umrissen verständlich zu sein scheint. 

Die Exoelektronenemission wird heutzutage mit den durch Elektronen besetzten 
Kristallgitter-Defekten einer vorangehend auf der Metalloberfläche erzeugten Oxyd- 
schicht gebunden. Wenn die erwähnten Kristallgitter-Defekte der Oxydschicht, die 
wie Haftstellen wirken, während einer Bearbeitung entstehen, können sie dann — wenn 
sie sich gleichzeitig mit Elektronen füllen — später als Elektronenquellen dienen (Sujak 
1959). Eine zweite Auffassung sieht den Mechanismus der photostimulierten Exo- 
elektronenemission aus mechanisch bearbeiteten Metalloberflächen als Begleiter- 
scheinung einer chemischen Reaktion (z. B. Oxydation) zwischen der frisch abge- 
deckten Metalloberfläche und der mit dieser kontaktierenden Gasatmosphäre (Sujak 
1959). In welchem Grade beide oben erwähnte Mechanismen miteinander verbunden 
sind, ist vom heutigen Standpunkt aus schwer zu beurteilen. Im Allgemeinen wird 
jedoch angenommen, dass die Oberflächendeformation, die entweder zur Erzeugung 
von Kristallgitter-Defekten in der Oxydschicht und deren Elektronenfüllung, oder 
zur Freilegung frischer Metalloberflächen durch Zerreissen der Oxydschicht führt, 
die Hauptrolle spielt. Im ersten wie im zweiten Falle kann also als Resultat einer De- 
formation dünner Oberflächenschichten die photostimulierte Exoelektronenemission 
einsetzen, wenn jedoch eine für das gegebene Material entsprechende elektromagne- 
tische Strahlung (Licht) zur Stimulierung benutzt wird. Die Oberflächendeformation 
von Aluminium wurde bisher durch Kratzen, Schaben oder Dehnen hervorgerufen 
(Grunberg und Wright 1957). 

In unserem Falle wurden die Aluminium-Proben sowohl oberflächlich wie auch 
innerlich deformiert. Die Oberflächendeformation beruhte auf einer mit Stahl durch- 
geführten Abschabung der Deckschicht, die innere auf einem Eindrücken eines 
Stahlstabes in die Probe. Wie sich zeigte, hat die Deformation des Innern grossen Ein- 
fluss auf die Intensität der photostimulierten Exoelektronenemission. In folge ver- 
schiedener Abklingungsgeschwindigkeiten der Emission konnten die einer plastischen 
Deformation des Innern unterzogenen Gebiete lokalisiert werden. Es wurde eine 
Reihe von Versuchen an Proben aus technischem Aluminium durchgeführt, die —aus 
einem Flachstück geschnitten — das Ausmass 60 x 30 x 8 mm besassen. Die Proben 
wurden bei ständig kontrollierter photostimulierten Koexoelektronenemission der 
Reihe nach einer oberflächlichen, einer inneren und wieder einer oberflächlichen 
Deformation unterzogen. | 
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a) Die Oberflächendeformation 


Nach Durchführung der Oberflächendeformation, die auf einer Abschabung mit 
Feile und Stahldrahtbürste einer ungefähr 0,5 mm starken Schicht beruhte, wurde 
die Probe auf einen verschiebbaren Tisch unter den offenen Geiger-Spitzenzähler ge- 
bracht (Fig. 1). Der Abstand Probenoberfläche-Zähler wurde beim Verschieben der 
Probe konstant eingehalten und betrug etwa 4 mm. Eine 6V-35 W-Wolfram-Lampe 


Fig. 2 


Fig. 1. Schematische Darstellung der Versuchsanordnung: H — Hochspannungsspeisegerät, V — Impuls- 
verstärker, I — Impulsintegrator, A — Schreibender Ampermeter, O — Oszilloskop, Z — Geiger- 
Spitzenzähler, L — Lichtquelle, P — Probe, T — Tisch, R — Papier-Registrierband 


"Fig. 2. Der Verlauf der photostimulierten Koexoelektronenemission von einer abgeschabten Aluminium- 


Probe. Kurve a unmittelbar nach der Bearbeitung gemessen, Kurven b und c entsprechénd 7 und 
14 Minuten spáier. 


beleuchtete durch eine Glaslinse und eine gleichzeitig als Zählerspitzenbefestigung 
dienende Plexiglaslinse die untersuchte Oberfläche. Die durch die photostimulierte 
Exoelektronenemission hervorgerufenen Zählerimpulse wurden nach Verstärkung 
durch ein Integrator mittels eines schreibenden Ampermeters registriert. Die Regis- 
trierrolle des Ampermeters war mit dem verschiebbaren Tisch gekoppelt. Wie aus 
Fig. 2 zu ersehen ist, ist die der untersuchten Probe entlang gemessene Intensität der 
photostimulierten Koexoelektronenemission nach den oben angeführten Eingriffen 
fast konstant. Die Emissionsintensität ist in willkürlichen Einheiten ausgedrückt. 
Die zeitliche Abklingung der Emission verläuft auch fast gleichmässig für alle Ober- 
flüchenbereiche der Probe. Kurve a entspricht einer Messung unmittelbar nach der 
Bearbeitung, Kurven b und c — 7 und 14 Minuten später. Während der Messung 
lig abgetasteten Flächenelemente beleuchtet. Die Emis- 


wurden immer nur die jewei 
sionsintensitat ohne gleichzeitiger Photostimulation konnte nicht ausgemessen wer- 


den, da sie in den Grenzen der Dunkelstosszahl des Zählers lag. 


b) Die innere Deformation 


wurde in eine der Seitenflächen der Probe 
die parallel zu dieser laufenden Seitenfläche 
" wurde die Probe vorher auf eine plane 


In der zweiten Bearbeitungsphase 
ein Stahlstab eingedrückt (Fig. 3). Um 
vor unregelmässigem Druck zu schützen, 
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Holzunterlage gelegt. Nach Durchführung der Stabeinpressung (es wurden Stahl- 
stäbe mit verschiedenen Profilen angewendet) wurde die Probe mit Feile und Bürste 
geebnet. Die gleich darauffolgende Messung des Intensitätsverlaufes der photostimu- 
lierten Koexoelektronenemission wurde entlang der in Fig. 3 eingezeichneten Linie 
a durchgeführt, die ungefähr 8 mm vom Boden der Druckkerbe entfernt war. Kurven 
der photostimulierten Koexoelektronenemission zeigten bei festgelegten geometri- 
schen Bedingungen verschiedene, von der nach Stabeinpressung verflossenen Zeit 


Fig. 3. Aluminium-Probe mit Druckkerbe. E — Einpressstelle, a — Messlinie 


abhängende Verläufe (Fig. 4). Kurve a entspricht einer Messung gleich nach der 
Ebnung mit Feile und Bürste (Oberflächendeformation). Nach 24-stündiger Lagerung 
ohne Lichtzutritt wurde Kurve 5 aufgenommen. Wie man sieht, prägt sich ein 
deutliches Maximum (W in Fig. 4) aus. Die zeitliche Gestalltung eines solchen Maxi- 
mums zeigt Fig. 5. 


Fig. 4 Fig. 5 


Fig. 4. Der Verlauf der photostimulierten Koexoelektronenemission von einer oberflächlich und innerlich 
bearbeiteten Aluminium-Probe. P — Probe, E — Einpresssteile, Kurve a unmittelbar nach der Bear- 
beitung gemessen, Kurve b 24 Stunden später. 

Fig. 5. Maximum-Ausprägung im Druckbereich einer oberflächlich und innerlich bearbeiteten Alumi- 
nium-Probe. P — Probe, E — Einpressstelle, Kurve a unmittelbar nach der Bearbeitung, Kurven b 
bis f entsprechend 6, 15, 21, 27 und 60 Minuten spáter gemessen. 


An einer frisch in obiger Weise bearbeiteten Probe wurde in verschiedenen 
Zeitpunkten der Intensitütsverlauf entlang der Linie a (Fig. 3) gemessen. Kurve 
a entspricht einer Messung unmittelbar nach der Bearbeitung, die Kurven b bis f 
entsprechend 6, 15, 21, 27 und 60 Minuten später. Wie schon vorher erwähnt, wurde 


See Ss eer 
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die Probe nur während der Messung beleuchtet. Wie man sieht, tritt mit der Zeit 
eine deutliche Selektivität der Emissionsintensitäts-Verteilung in der Umgebung der 
Einpressstelle (innere Deformation) ein. Die Emissionsintensität klingt also langsamer 
in unmittelbarer Umgebung der Einpressstelle als in übrigen Oberflächengebieten 
der Probe ab. Wenn als Messlinien parallel zu a (Fig. 3) laufende Linien gewählt 
werden, kann der Wirkungsbereich der inneren Deformation in Abhängigkeit vom 
Profil der eingepressten Stäbe ermittelt werden. Die Aussetzung der innerlich und 
oberflächlich deformierten Probe auf längere Lichteinwirkung (z. B. Tageslicht) be- 
wirkt ein fast vollständiges Herabsinken des Emissionsvermögens. Nach wiederholter 
Abschabung der dünnen (ungefähr 0,3 mm) Oberflächenschicht einer solchen Probe 
wurde keine Entwicklung einer selektiven Intensitätsverteilung der Emission mehr in 
der Umgebung der Einpressstelle festgestellt. 

Die durch Knicken und darauffolgendes Gerademachen sowie durch seitliche 
Einpressung von Stahlspitzen erreichte innere Deformationen führten zu ähnlichen 
Ergebnissen. Der sichtliche Einfluss der inneren Deformation auf die Exoelektronen- 
emission von oberflächlich bearbeitetem Aluminium kann auf zweierlei Weise gedeutet 
werden. Erstens als Herstellung im deformierten (Metall!) Volumen gewisser, nach 
der Oberfläche diffundierenden Zentren, die danach selbst als Exoelektronenquellen 
auftreten oder durch lokale Energieabgabe auf die Exoelektronenemission aus anderen 
Zentrenarten einwirken. Für diesen Mechanismus würde die Tatsache sprechen, das 
sich das Emissionsmaximum in der Umgebung der inneren Deformation unabhängig 
(in erster Näherung) von der Stärke dernach Deformationsvollführung abgeschabten 
Oberflächenschicht ausprägt, sowie auch das Ausbleiben des Maximums nach völliger 
„Bleichung“ durch Lichteinstrahlung und darauffolgender Oberflächenbearbeitung. . 

Der zweite mögliche Mechanismus einer Ausprägung des angeführten Maximums 
würde auf einer durch die innere Deformation hervorgerufenen Kohäsionsänderung 
der Oberflächenschicht beruhen. Diese Kohäsionsänderung könnte sich in einer mikro- 
porösen Oberflüchenstruktur der Deformationspunktumgebung, also in Herstellung 
von efektiv grösserer „frischer Oberfläche“ auswirken (siehe auch Gobrecht 1954). 
Der erschwerte Transport von Sauerstoffmolekülen in einer und Photoelektronen 
(Sauerstoffionen) in entgegengesetzten Richtung würde die langsamere Oxydation 
der mikroporösen Oberfläche im Vergleich mit der übrigen Probenoberfläche beein- 
Aussen. Dieses würde das langsamere Abklingen der Prädisposition zur photosti- 
mulierten Koexoelektronenemission aus der Umgebung der Einpressstelle begrün- 
den, da die übrige Probenoberfläche nur einer seichten Oberflächenbearbeitung 
unterzogen wurde. Für diesen Mechanismus sprechen die Tatsachen, dass die 
Exoelektronenemission nicht im Hochvakuum auftritt (Lohff und Raether 1955, 
Lohff 1956), wie auch das verschwinden der Emissionsfähigkeit von oberflächlich 


bearbeiteten Aluminium-Proben nach Einwirkung einer Hochfrequenzentladung 
Jin Luftatmosphäre (Sujak und Mader 1958). Die Entscheidung, welcher yon diesen 
Mechanismen die Hauptrolle spielt, wie auch der Mittwirkungsgrad beider ist fortge- 


hend Gegenstand weiterer Untersuchungen. 
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Einfluss der Deformation von Polimeren auf die Dunkelstosszahl des Zählers 


Die Feststellung einer plastischen Deformation in Aluminium mit dem Geiger- 
Zähler beruht auf Intensitätsmessungen der Koexoelektronenemission. Bei Polimeren 
wurde jedoch bisher keine während oder nach einer plastischen Deformation ein- 
setzende Exoelektronenemission festgestellt. Es kann sein, dass eine photostimulierte 
Exoelektronenemission aus einer deformierten Polimer-Probe erst nach entsprechend 
energetischer Ultravioletteinstrahlung einsetzen würde, die Einstrahlung von Wellen- 
längen unter 2000 Ä während der Zählermessungen ruft jedoch eine Photoelektronen- 
emission aus der Zählerkathode hervor, was die Feststellung einer eventuellen Exo- 
elektronenemission aus der Probenoberfläche gänzlich unmöglich macht. Um also 
die Deformation einer Polimer-Probe mit dem Spitzenzähler festzustellen und auszu- 
messen muss eine andere Erscheinung als die Exoelektronenemission ausgenutzt 
werden. Während der Deformation durch Eindrücken entstehen Kontakt-Elektri- 
sierungen (Triboelektrizität) sowie auch elektrische Polarisierungserscheinungen, die 
dann mit dem Zähler abgetastet werden können (vergleiche auch z. B. Bazan 1959). 

= Die auf solche Weise elektrisierten Gebiete einer deformierten Polimer-Probe 
können dank der Eindringung ihrer elektrostatischen Felder in das Innere des Zählers 


i 


Fig. 6. Elektrisierungsverteilung einer Plexiglasplatte nach der Druckdeformation. a) Messprobe mit 

eingezeichneten Messlinien I, II und III. Die Pfeile geben die Druckpunkte an. b) Kurven Ib, IIb und 

IIIb entsprechen der Dunkelstosszahl-Änderung des mit einer dielektrischen Kappe geschlossenen 

Zählers, entlang der Messlinien I, II und III gemessen. Kurve IV — Kontrollmessung vor der Druck- 
i : ; deformation. f L 
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qualitativ festgestellt und verglichen werden. Dieses Problem war schon Gegenstand 
früherer Veröffentlichungen (Sujak 1957a, b; Lewowski und Sujak 1959). Dort wurde 
das Arbeitsprinzip wie auch eine Reihe von wesentlichen Versuchsergebnissen über 
die relative Messung von Verteilungen der Elektritizätszustände auf Oberflächen 
deformierter Isolatoren angegeben. Der Geiger-Zähler musste, um das Ionenerzeu- 
gende Gebiet (Zählvolumen des Zählers) von der untersuchten Oberfläche zu trennen 
mit einer dielektrischen Kappe abgeschlossen werden. Die Versuchsanordnung war 
der in Fig. 1 dargestellten ähnlich. Als Beispiel wird die Elektrisierungsverteilung 
einer Plexiglasplatte (62 x 18 x 6 mm), die vorangehend einer schematisch in 
Fig. 6a gezeigten Druckdeformation unterzogen wurde, angeführt. Die deformierte 
Plexiglasplatte wurde unter dem Geiger-Spitzenzühler (Abstand Probenoberfláche- 
Zähler etwa 2 mm) entlang der Linien I, II und HíI (Fig. 6a) verschoben. Der Zähler 
reagierte auf die lokale Elektrisierungszustände mit einer Änderung seiner Dunkel- 
stosszahl wie es Kurven Ib, IIb und IIIb in Fig. 6b zeigen. Die Kurven wurden — 
wie im ersten Teil— durch ein schreibendes Ampermeter auf einem Papier-Regis- 
trierband aufgetragen. 

Die untersuchungen von Polimer-Deformationen mit dem Geiger-Spitzenzähler 
scheinen besonders Vorteilhaft im Falle einer optischen Undurchsichtigkeit des Mate- 
rials der untersuchten Probe zu sein. 
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The purpose of this paper is to ascertain, whether a spin wave resonance effect in 
antiferromagnetics analogous to that predicted by C. Kittel (Phys. Rev., 110, 1295 (1958)) 
for ferromagnetics is possible. By means of semiclassical formalism, the author shows 
the possibility of exciting odd modes of spin waves in antiferromagnetics and calculates 
the resonance values of applied magnetic field for these modes. Each mode can be excited 
by two different values of applied field. The investigated phenomenon might be probably 
observable in films 3—10 times thicker than in the case of ferromagnetics. 


Introduction 


Some months ago Kittel (1958) has predicted the possibility of exciting exchange 
modes of standing spin-waves in a thin ferromagnetic plate by a uniform radio- 
frequency field. He has proposed the term *'spin-wave resonance” for this effect. The 
basic assumption which leads to the exciting of spin-wave modes is the spin 
pinning boundary condition, proposed by Kittel. Here the spins are assumed to be 
pinned down at the surface of the specimen, owing to the magnetic anisotropy energy 
which must be much greater at the surface than in the interior. Consequently, the 
spin waves can be considered practically as the sinusoidal standing waves with nodes 
at the boundaries. It follows, that only odd modes of spin waves are allowed. 

The effect of ferromagnetic spin-wave resonances has been observed by Seavey 
and Tannenwald (1958) for permalloy film of 5600 A thickness. The results agree in 
principle with Kittel’s predictions. ee 

In the present paper we attempt to investigate the analogous ote in antifer- 
romagnetics. We use the approach similar to that of Kittel, although with some modi- 
fications. We apply indeed the semiclassical theory of spin waves, nevertheless we 
consider the discretely localised spins. Moreover, instead of linearly polarized radio- 
frequency field we take into account only one of his two circularly polarized compo- 
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nents, namely the one which rotates in the sense of spin precession. Such approach is 
suggested by the paper of Bloch and Siegert (1940), who have shown, that the other 
circularly polarized component causes only an inobservably small displacement of 
resonance frequency. We get two advantages from such an approach. First, the semi- 
classical equation of spin motion is strictly satisfied by wave-precessional solution 
(otherwise, the equation for S* is satisfied only approximately). And second, the itera- 
tion made by Kittel is not necessary in that case. This iteration of simultaneous differ- 
ential equations for S* and S leads to two series of spin-wave resonances for each œ 
(o = + e, in the notation of Kittel), one of which must be ficticious. 


Derivation of the resonance formula 

Let us consider an antiferromagnetic plate with parallel plane faces. Taking 
je jy. J, as unit vectors of carthesian coordinate system, we choose J, and J, parallel 
to the plate surface. As mentioned above, from the radiofrequency field 

2Hge™ j, = Are™ (j, T ij) + Hae (j, z ij) (1) 

we retain the first term only (therefore our Hp corresponds to Kittel’s H,/2). Further- 
more we take both the anisotropy field H , and the applied field H in the direction of 
z-axis. Y 

In anlogy to Kittel`s quoted paper we consider only a chain of N spins with spacing 
a aligned along the z- axis. Then, the equation of motion of m-th spin is: 


S. — Ss. x (2Jjh) (S. + Sur) d > - 
+ S_ x y [GP + (H A") je + Hg e™ (ie + 3 Qo 


J being the exchange intergral between nearest neighbours, the coefficient y < 0 — = 
the magnetomechanical ratio and A — — the Planck’s constant divided by Zr. We take — 
the solution of (2) in the form: 4 caes ja 
a PLE A 921 
Wa aP AER co NI ey NS EE T Boi fox 

- va x oni 
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= where 
b, = Ry [|o + y (H + Hj) + n) +R n cos ak, 
C, = Rj q cos ak, — RjH [o + y (H — H j) + ri] 


= 45J/h 
We see that the coefficients of sines must be correspondingly equal: 
k ,=% 
which leads to 
Q; Q; 
b, = —P—- er 
p R Mi Re ua (5) 


where 
Q, = (yH — oy + YH} — n [2y H4 + n cos ak, (2 — cos ak,)] 
MJ = w + y H — yH,—1 eos ak, 
M," = — o — yH — yH, — cos ak, 
We shall calculate the coefficients b, multiplying both sides of (4) by sin (gz m/N) and 
summing from m = 1 to m = N. 
Using the relations (see Appendix): 


reat AE, ' 7 N 
Y sin (pam/N) sin (gam/N) = > p, (6) 
m=1 - 
a 7 doctum 
Ye |, aq PETER @ 
^^ m-1l 
we obtain 
spe q 0 for p even 
| by = [asp Bo. for p od | 


er ters cr 


190 Hreofta 


where 


Hg = ny > 0 


For small p we can take: cos? ak ~ 1 — a?k? 
and therefore 


H, = H, + VH,4(H, + 2H5) + Cp? (10) 
where 


C-—HgmnuN* and H, = —o@|y > 0 


ae 


For each p there occur two resonance values of applied field. For p — 0 we get the 
well-known formula for the main resonance. 
Let us introduce the abbreviation: H4 (H4 + 2Hg) = B. Then, the interval 
between two successive odd-mode peaks is 
4C(p + 1) 
VB + C(p + 2)? + VB + Cp? 
s AC (p 4- 1) 
2yB + Cp? 
The effects of anisotropy fields are usually small: H4, —10-? H, Taking N ~ 10? and 
p «20 as in experiments of Seavey and Tannenwald, we have approximately 
B + Cp? ~ 2H , Hg. Consequently we get: 
2H: n? 
N?V2H,4 Hg 


AH, = | Hp+2 — H,| = 


Comparing this with the value of corresponding interval for the ferromagnetic case: 
, m? 
A Her = 2Hg W2 (p + 1) 


where Hj = 4JS/|y|h (J < 0), we see that for the same |J|, S, a and N values the fer- 
romagnetic resonance peaks lie V H,/2H, times nearer than in our case. Hence, the 
antiferromagnetic observations can be made using 3—10 times thicker films than 
in experiments for ferromagnetics (provided that the absolute value of exchange 
integral has the same order of magnitude in both cases). 

It seems, that actually the most suitable antiferromagnetic substance for such 
experiments should be chromium. 
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Appendix 


One can verify the equality (9) as follows: 


m pm mam 1 j = 
s= s ae 29 NE q) am _ _ 


Each of four sums is a finite geometrical series, so that 


Ants exp [ix (p — q)] — 1 exp [— iz (p — g] — 1 
1 — exp [— iz (p — g||N] ` 1— exp [iz (p — DIN] 
exp [ix (p + g)] — 1 exp [— iz (p + gj] — 1 


~ T—exp[— iz (p + gJN] 1 — exp [iz (p + QIN] 
For p — q both the first and the second term tend to N and the other two terms tend 
to zero. For each p Z q however we have 


gi" (p— g) — e 0-29 — ei^ Ota) — e 0*0 — Pi 


For the upper sign we have s = 0; for the lower one the considered sum s also vanishes 
because of the identity 
l 1 


1 — e* 1 — e“ 


II 


Therefore, the relation (9) is verified. 
We can treat the sum (10) approximately because always N>p>1 for the 
spin wave resonances phenomena. We have then 


in - sin aa ee (N + D sin? 7 


N sin — 
2 2N 
1 — — = 
) sin (prrm/N) pn pa 
m=1 sın 2N 2N 


For p even, sin? (p/2) vanishes, for p odd however, sin? (p/2) = 1. Therefore, the 
relation (10) is also verified in a good aproximation. ‘ 
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NOTE 


has read in Physical Review Letters the summary of 


After preparing this text to print the author s the 
an article by R. Orbach and P. Pincus, to be published in Physical Review and treating the same pro- 
blem. It is to be hoped, that the results of Mrs Orbach and Pincus will agree in principle with our results. 
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The purpose of the present paper is the theoretical investigation of the possibility 
of exciting spin wave modes in ferrimagnetics by an uniform rf field. Ferrimagnetics 
with parallel spins as well as those with antiparallel spins have been considered. The reso- 
nance values of applied magnetic field for allowed odd modes of spin waves and fora given 
value of the frequency of rf magnetic field are obtained. For even modes the resonances 
do not appear provided that eddy currents can be neglected and that spin waves have 
nodes at the boundaries. In the special case of ferromagnetics one gets the formula of 
Kittel (Phys. Rev., 110, 1295 (1958)) and in the case of antiferromagnetics the resonance 
formula goes over into that derived by Cofta (Acta phys. Polon. 97, 187 (1960)). The 
author estimates also the interval between successive odd mode peaks. 


Introduction 


In the past year Kittel (1958) has investigated theoretically the possibility of ex- 
citing exchange modes of spin waves in ferromagnetics by an uniform radio-frequency 
magnetic field. He showed that such effects may occu rin thin specimens. In this 
manner a new resonance phenomenon, called “spin-wave resonance”, has been discov- 
ered by means of theoretical research. Seavey and Tannenwald (1958) have observed 
the spin-wave resonances for a permalloy film of 5600 A thickness. These experi- 
ments have confirmed Kittel’s results; as predicted, pronounced peaks occur for odd 


modes only. 


The direct observation of spin-wave resonances 


the spin waves method, even in the simple semiclassical form as used by Kittel, is 


heory of ferromagnetism for an unexpectedly large range of 


a good approach to the t | 
temperatures. Moreover, from such’ measurements one can obtain the value of ex- 


at room temperature shows that 


change integral. 
On account of these reason 
discussion of spin-wave resonance to a 


s it should be very interesting to extend the theoretical 
larger class of magnetic materials. An attempt 
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of such a discussion for the case of antiferromagnetics has been made recently by Or- 
bach and Pincus (1959) and independently by Cofta (1960). 

In a simplified manner anlogous to that of Kittel’s quoted work, we shall consider 
in the present paper the general case concerning all substances for which a significant 
exchange coupling between atomic spins exists, inducing an arrangement of spin 
orientations. We shall call all such substances “ferrimagnetics” and include here these 
with parallel spins as well as these with antiparallel ones. When all spins have the same 
value and are arranged parallelly — we have the case of ferromagnetism; however 
for antiparallel orientation of such spins —- we have the case of antiferromagnetism. 
For both these special cases the resonance formula derived in our paper goes over 
into ferromagnetic or antiferromagnetic formula respectively, both derived in the 
papers! mentioned above (Kittel 1958, Cofta 1960). 

It is perhaps worth while to remark that permalloy, which has been used by Seavey 
and Tannenwald (1958), is strictly speaking a ferrimagnetic. 


Calculations 


We shall use here the semiclassical method similar to that of Kittel (1958). Some 
modifications introduced have been discussed in our previous article (Cofta 1960). 

We shall consider a ferrimagnetic chain of spins, which is directed perpendi- 
cularly to the plane faces of a plate. Let us assume, that the chain contains N spins 
of two kinds, aligned .alternately in the equal distances a. In a carthesian coordinates 
system we can write the equation of motion for m-th spin S,, of the considered chain 
as follows 


AS, = S, x {hy G,, + 2] (S,,-1 + Sn+1} (1) 
where J is the exchange integral for nearest neighbours interactions and h — the 
Planck's'constant divided by 2%, the negative constant y being the magnetomechanical 
ratio. Let us assume the z-axis perpendicular to the plate face and denote by j,, j, 
and j, respectively the unit vectors of x, y and z-axis. Then. the vector G,, includes 
all three magnetic fields which we take here into consideration, namely: 1? the ap- 
_ plied constant field H, directed along the z-axis; 2° the effective anisotropy field H gm ` 
which acts on the internal spins also in the z-direction and may be different for different 
sublattices; 3° the rf field 2H ge" j, alternating with a given frequency of order of 
10* Me/sec or more. We can decompose this field into two circularly polarized compo- 
nents and neglect (Bloch and Siegert 1940) that component which rotates in opposite 
sense to the spin precession (this point is also discussed by Cofta 1960). Then we have 

G,, EN ur a GL j: w Hg ej 
where j = j, + ij,. 


1 Writing the present article the author did not know yet the results of Orbach and Pincus (1959). 


It is to be hoped, that they will agree with the formula for antiferromagnetic special case, derived in the 
present paper. . ê ' 
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When we impose Kittel’s (1958) spin pinning boundary condition then the solu- 
tion of (1) represents a standing wave of spin precession. We write it in the form: 


Sm = Rp ej + Sy, j, (2) 
Substituting (2) into (1) and using the identity 
we get the fundamental relation for R,,: 
Ra 0 [LU le Sees) = 
== $177 55, (ine setae tee AT aE (3) 


where n = 4J/h. The R,,, Sh and (H ‚)„ are different for different sublattices; therefore 
we put for the first one: 


SA ORY =D) 4, sincak ms (Hm Hs. 
P 


and for the second one: 


Se =B, R,= 2 b,snak m, (H), = H; 
P 


where the wave number k, = z plaN on account of Kittel’s boundary condition. 


Moreover, we have to put H* = — H. Hence, we get for each sublattice a relation of 
the same form: 
D o, sin ak, m = y Hg (4) 
P 
but we have 
Ae, — s= 0 + y (Hi — H) + n] — bm 4° f, (5) 
for the first sublattice and 
(0 Bre, = b, [-0 + y (Fg — H) na — aos P7, (6) 


for the second one, f, being the abbreviation for cos akp 
From (5) and (6) we get: 


Cy = ap SE = be Hf (7) 


where 
Q,- [o y 0t — B) E p Ba o y GR MtnAl- AB 
L,— A* [7 o 4- y (Bj — H) + nA" + Pf) 
M, — B, [o + y (Hi — H) + n B* + nA* f.) 


One can obtain the coefficients c, multiplying both sides of (4) by sin qamjN and 


making the summation from m = 1 to m = N (as to the details see the article of 
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Cofta 1960). We get: 
"m 0 for p even 
^ ^ My Hg[np for p odd 


and therefore from (7) we obtain for p odd: 


_ Ha ‚I 

oo TP Q; 
(8) 

en stp Qp 


both these coefficients vanishing for p even. 


Discussion 


The result (8) shows that when the spin pinning boundary condition is fulfilled 
and it is possible to neglect eddy currents, then the odd mode peaks only can be 
expected. As can be seen from (8), resonance occurs for Q, = 0. Taking œ constant 
we get quadratic equation for H, from which we obtain the resonance values of 
applied field for allowed modes: 


H, = H, ++ + HD +5 (niy) (47 + B) + z (In (du Br) + 
+ y (H; — 2 + 4 A* B* n? cos? Jy (9) 
where H, = — w/y > 0. For low values of p we have to write: 
H, = H, x x: Va — Bp (10) 
where 
x = $ (Hi + H3) + š (nly) (4* + B°) 
a = (n/p)? (4* + B9? + 2(n/y) (A? — B) (Hz — HD + (H5 — HD: (11) 
B = 44* B* (ny)? z°JN2 
Now we shall consider the special forms of derived formula for various cases. 
l. The case of ferromagnetism 


In this case we have to put 4* = B* = S and Hj = H; = — H 4. Then, intro- 
ducing the abbreviation 


Hg = Inlyl (12) 
we get from the appropriate branch of (10): 
H, = H, — H4 — 4 Hs n p|N* Mri 


which agrees with Kittel's (1958) result!. 


1 The double sign + ® before the applied radio-frequency appearing in Kittel’ s formula is caused 
by the iterating procedure used by this author. 
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2. The case of antiferromagnetism 
Here we have to put A* = — B* = S and Hj = — H; = — H,. 
Then 
H, = H, + [H,(H, + 2Hg) + Hg a? p2|N2] (14) 
which is identical with the result of Cofta (1960). 
3. The case of ferrimagnetism with parallel spins. 
In this case we put: Æ” = S, B*— S, and H = — M, Hj — — H, 
Then we get for (10): 
x = — š (H, + Hy) + š (nly) (S, + S) 
a = (nly)? (S, + Sa)? + 2(n|y) (Hı — He) (Sı — $9) + (H, — H)? (15) 
B = + Any)? S. S, z°|N° 
4. The case of ferrimagentism with antiparallel spins (antiferri- 
magnetism). 


We put: A% = Si B* = —S,, = — Hy H3 — H,, getting for (10): 


x = $ (H, — H) + $ (ly) (S, — 59 
a = (nly)? (S, — Sa)? + Anl) (H, + Ha) (S, + S) + (H, + Ha)? (16) 
B = — Amy)? S, S, 2P|N | 

For both cases of ferrimagnetism (case 3. and 4.) we have «> B p? for small p values 


and N — 103; consequently we can expand the square root in (10). Then; we get 
the resonance formula (10) in the form 


H,~ H, + x + 1 a + 3 pa% p? (17) 


similar to that for the case of ferromagnetism. 

As can be seen, in cases 2., 3. and 4. we have two resonance series, because for 
each p two resonance values occur, separated mutually by (x — Bp?)*. In cases 3. 
and 4. the separation of two peaks corresponding to the same value of p may be 


written approximately : 
| IH; — Hy | = a — $ Ba" p° (18) 
Let us now estimate the interval i 
AA SH ae ae 
between two successive odd mode peaks of the same series. Using the abbreviation (11) 
we can write for low values of p: 
AH — 3 lla — B (p + 3% — [a — Br” 
= |p| (a — Bp?) (p + D 


Let us put: Sı = S, Sa = 45, where usually 4 < 1 < 1 Considering not too thin 
plates (N ~ 10°) and neglecting the generally small effects of anisotropy fields, we 
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have approximately for the lowest modes (p < 10): 
a — pp? Hal) 
for the case of parallel spins and 
a — Bp? ~ HB (1 — 2)? + AH, Hy (1+) 


for the case of antiparallel spins, where Hg is defined by (12). Hence for parallel spins 


the considered interval is 


NPE PEE (19) 
Ne TA Oa 


which goes over into 


T2 
2Hg sa (p + 1) 


for the case of ferromagnetism (A = 1). For antiparallel spins however we get 


AHg À 


at [1 — 2)? + 4(1 + 4) Hal Hg]? 


D (20) 


which goes over into 
2Hg T? 
Von]; Ni TY 


for the case of antiferromagnetism (A = 1). 
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The article on antiferromagnetic spin-wave resonance by Orbach & Pincus has appeared in Phys. 
Rev. 113, 1213 (1959). The results obtained by these authors are quite the same as those shown in the 
formula (14). The author is greatly indebted to Dr. R. Orbach and Dr. P. Pincus for kindly sending the 
informations on their publication prior to coming of No. 5 (March 1, 1959) of the Physical Review to 
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The absorption spectra and decay curves of phosphorescence of trypaflavine in gelatine 
were investigated as a function of the pH value of the aqueous trypaflavine solution in 
which the gelatine has been coloured. It appeared that they vary with pH of the above 
solution. The observed decay curves and their decomposition into simple exponential 
functions can be explaned assuming the trypaflavine molecule to exhibit one metastable 
level only and by taking into account the observed preexcitation effects. 


Introduction 


The aim of the present paper is to check on ground of investigation of the decay 
curves of phosphorescence of trypaflavine in gelatine, whether this dye exhibits more 
than one metastable level. | 

On ground of investigations on the decay curves of phosphorescence of trypafla- 
vine adsorbed on silicagel (Pringsheim and Vogel 1936) and on observations on 
oxygen quenching, the trypaflavine molecule was supposed to possess two metastable 
levels (Pringsheim and Franck 1943a). Similar suggestions were put forward from 
considerations of the observed variations of the polarization rate in various spectral 
ranges of the phosphorescence and slow fluorescence bands (Pilipovitch and Sveshni- 
kov 1958a), of deviations of phosphorescence decay from the exponential law, and of 
the dependence of the slow fluorescence decay constant on the wavelength within 
the emitted band (Pilipovitch and Sveshnikov 1958b). 

Molecules, for which investigations carried out hitherto have lead to ascribe 
them more than one metastable level, are assembled in Table IL However, the 
interpretation of the above facts is by no means unequivocal. Thus, e. g. Bre- 


1 The structure of the dye molecules given in Table I is to be found in Pringsheim and Vogel's 
«Luminescence of Liquids and Solids” (1943b), and in Foerster’s “Fluoreszenz organischer Verbin- 
dungen” (1951). I 
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TABLE I 


Some dyes supposed to possess more than one metastable level 


Substratum 


Jl Silicagel Pringsheim, P., Vogels, H., 1936 
Pringsheim, P., Frank, J., 1943 
Sugar Pilipovitch, V. A., Sveshnikov, B. I. 
1958 a, 1958 b 
2 Fluoresceine Boric .acid Pilipovitch V. A., Sveshnikov B. I., 1958b 
Sugar Pilipovitch V. A., Sveshnikov B. I., 1958 a 
3 acridine orange Sugar Pilipovitch V. A., Sveshnikov B. I., 1958 b 
4 esculin Sugar Pilipovitch V. A., Sveshnikov B. I., 1958 b 
5 auramine Sugar Pilipovitch V. A., Sveshnikov B. I., 1958 b 
6 carbosol Sugar Pilipovitch V. A., Sveshnikov B. I., 1958 a 
U eosine Sugar Pilipovitch V. A., Sveshnikov B. I., 1958 a 


del (1955) explains the deviations from the exponential law observed in the decay 
of phosphorescence of fluoresceine in boric acid by a hypothesis due to To- 
maschek (1922). The electronic terms of the dye molecule and the probabilities of 
electron transitions connected with phosphorescence emission can vary as a result 
of interactions between the dye molecule and its surroundings. The interactions may 
by different in the case of various molecules. Hence phosphorescence becomes a com- 
plex process and its decay will no longer be described by a purely exponential law. 
These deviations can also be explained, i. a., by a simplified model of luminescent 
centres as proposed by Jablonski (1957a). 

Quite generally, the interpretation of the observed luminescence phenomena 
depends largely on whether differences in the nearest surroundings of the various 
luminescent molecules are taken into account or not. This is a point that is especially 
critical in the various theories of photoluminescence. The last standpoint is represented 
by theories put forward by Vavilov (1942, 1943, 1944, 1950, Feofilov and Vavilov 
1931, Shishlovsky and Vavilov 1934), whilst the other is that of Jablonski’s theory of 
luminescent centres (1936, 1954, 1955, 1957a), by which the observed phenomena 
should be generally considered as the superposition of phenomena occuring within 
the various groups of centres. In very special cases only do boih theoretical ap- 
proaches yield the same form for photoluminescence decay. 

The above differences in theoretical approach are practically of little consequence 
for the problem of the dependence of photoluminescence yield on the concentration: 
of the quenching molecules (Förster 1949, Galanin 1955?, Jablonski 1957 b). 


* Rozman (1958) proved the theories of Förster (1949) and of Galanin (1955) to be identical. 


ae 


for an hour, and subsequently rinsed 
the decay curves were investigated. The try 


0,1 n HCl solution, and of 0,1 n KOH 
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So far, no suggestions have been made on the possibility of solving the problem of 
multitude of metastable levels on ground of investigations of the dependence on 
concentration of the luminescence yield. 

On the other hand, Kislak (1958, 1959), who investigated the slow fluorescence 
decay of above organophosphors, has found it to be strictly exponential. This amounts 
to assuming that a single metastable level is responsible for the slow fluorescence 
observed by this autor. 

It seems to be a problematic question whether there are or not two metastable 
levels in trypaflavine as well as in fluoresceine. The aim of this paper is to elucidate 
this question. 


Experimental 


The decay curves of trypaflavine in gelatine were investigated with a “linear” 
phosphoroscope. A phosphor having the shape of an “driving belt" was set in motion. 
A part of its path was rectilinear. In this segment of the path, the phosphor was excited 
and the decay of phosphorescence measured. An appropriate system of slits served 
to separate the exciting and scattered light from that of the luminescence measured. 
The phosphorescence light was photographed through a polaroid. By orienting the 
latter with respect to the direction of the electric vector of the plane polarized 
exciting light, it was possible to obtain successively photographs of both the J, 
(perpendicular) and Jy (parallel) components of phosphorescence. A detailed de- 
scription of the “linear” phosphoroscope is given in an earlier paper (Frackowiak 
1957a). 


The phosphorescence decay curves were measured in the following conditions: 


1) velocity of displacement of phosphor: 1155 + 60 cm/sec, 

2) duration of excitation: 2,16 x 10-4 sec, 

3) observation of the decay curves started 4,76 x 10-* sec, after cessation of 

excitation and lasted until 10°? sec, 

4) duration of exposition was 30 sec, 

5) measurements were carried out at room temperature (20°C). 

excited with plane polarized white light. Both components 

of the decay curves were investigated throughout the entire band of luminescence. 
The phosphors were prepared as follows: Unexposed fixed and dried “Ultra- 


rapid” films of equal length were kept in appropriate aqueous trypaflavine solutions 
in flowing water for 10 minutes. After drying, 
paflavine concentrations in the phosphors 
which were measured with a Zeiss spectro- 


The phosphor was 


were determined from the absorption bands, 
photometer. 

. The above mentionned tryp 
and 277,7 mg of trypaflavine were 


aflavine solutions were prepared as follows: 27,7 mg 
dissolved in equal volumes of distilled water, of 
solution. Thus, solutions of known trypaflavine 
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concentrations and different pH were obtained. The latter was measured with a 
PHM — 3i pH-meter. The absorption spectra of the aqueous trypaflavine solutions 
were measured with a Zeiss spectrophotometer. 

All data relating to the solutions and trypaflavine concentration in the phosphors 
are given in Table II. 


' TABLE II 


Quantities characteristic of aqueous trypaflavine solutions, and concentration in gelatine 


nai Re VONT — Solutions in gelatine 
No. of Trypafllavine con- Trypaflavine 
solution Solvent centration pH Y concentration mol/l 
x 107? mol/l 
P 

1 0,1 n HCl 0,106 1,0 7,38 

2 1,066 0,5 2,44 0,0071 

3 distilled water 0,106 4,6 1,0 0,045 

4 1,066 19. 0,93 0,444 

5 0,1 n KOH 0,106 1275 2,66 0,057 

6 1,066 10,2 0,87 0,436 


Experimental Results 


Fig. 1 shows the extinction coefficients as a function of wavelength, for trypa- 
flavine in the solutions used for impregnating gelatine with the dye. All extinction 
curves in Fig. 1 have the same area, i. e. the total “intensity of absorption" is equal. 
By comparing these areas with tbat of a band for an arbitrarily chosen concentration, 
the relative value of the factor y in the generalized expression for the oscillator strength 
(Proter and Windsor 1958), accounting for interaction between the absorbing molecule. 
and the medium, can be determined: 


oo 


mc? mc? 
(= 25 [voor s [o nh (1) 
0 - 


0 


where f — denotes the oscilator strength, 
m — the mass of the electron, 
c — the velocity of light, 
e — elementary charge, 
N — the number of absorbing molecules per 1 cm? of the solution, and 
& — the extiction coefficient. 


For the lower concentration of the trypaflavine in distilled water y=1l(f=f,) was 
assumed. The values of y for different solutions are given in Table II. 
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a 600 . mi 


Fig. 1. Absorption spectra of ayueous trypaflavine solutions. 


1 
2 
3 
4 } see Table II 
5 
6 


for the absorption band of trypaflavine in gelatine. These 
have been *normalized" as in the case of the aqueous solutions. From these absorp- 


tion bands, the concentrations C of trypaflavine were determined, as given in Table II. 


It should be stressed that the C values do not necessarily coincide with the true concen- 


trations of trypaflavine in gelatine. As the swelling process takes place in media with 
different pH, gelatine can exhibit different properties (cf. Basinski, “Physical Che- 
mistry of Colloids” 1957). These variations can lead to variations in the interaction 
between trypaflavine and gelatine, and, hence, to variations in the value of the factor y 
in eq. (1). Hence, determining the concentration from the absorption bands, with the 


approximation of eq. (1), yields but a value C* equal to 


Fig. 2 shows curves 


cr=yt, (2) 


, 
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300 400 500 2 600 mu 


Fig. 2. Absorption spectra of trypaflavine in gelatine. 


—..—..—..—.. Phosphor No. 


2 
BOT Ee A 3 
oe 55 4 f see Table II 
d 5 

6 


” 


with C denoting the true concentration of the dye, and y — the factor accounting for 
interaction between the dye and the substrate. 

Figs. 3, 4, 5, 6, 7 and 8 show decay curves for the J, and T, components of the 
phosphorescence of trypaflavine in gelatine, together with their decomposition in 
exponential functions. 

The decay curves for phosphors impregnated in solutions Nos. 3, and 6 (Figs. 5 
6 and 8) are given by expressions: 


> 


Ij (t) = 4 (1 + 2r) exp (— pı t) + A, (1 + 27’) exp (— gat), (3) 
I, () = A, (1 — r) exp (— 9, t) + A, (1 — r”) exp (— 9s t), 


wherein A,, Ag, r, r', 94 and p, are certain constants, and £ — the time. 
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in J 


Fig. 3. Decay curves for trypaflavine in gelatine. Phosphor No. 1 
—e—— parallel components 


perpendicular components 


° 
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Fig. 4. Decay curves for trypaflavine in gelatine. Phosphor No. 2 
e—— parallel components 
perpendicular components 


o: 
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Fig. 5. Decay curves for trypaflavine in gelatine. Phosphor No. 3 
e—— parallel components 
——o—— perpendicular components 
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6 = 
Fig. 7. Decay curves for trypaflavine in gelatine. Phosphor No. 5 
ə parallel! components 
E o perpendicular components 
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\ Fig. 8. Decay curves for trypaflavine in gelatine. Phosphor No. 6 
5 € parallel components 


perpendicular components 


° 
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Hence the time dependence of the total phosphorescence intensity emitted in all 
directions is: 


M(t) = I, (t) + 2I, (t) = 3[4, exp (— pı t) + Ap exp (— pa t)]; (4) 


the values of the constants are assembled in Table III. 


TABLE III 


Constants of phosphorescence decay, emission anisotropies, and initial phosphorescence intensities 
for trypaflavine in gelatine 


No. of Decay constants Emission Values of J); and T| 
: Component "epe : š 
solution x 10? sec anisotropies at £= 0 


I) 34,0 — 75 
1 
1 D = = 2 
Ij 7,8 — 26 
uu" d — 12 
I 40,0 — 245 
2 n 50,0 215 
Ij 3,8 — 42 
P 7,5 24 
I 22,5 0,232 305 
3 iv 22,5 160 
f 2,25 0,139 95 
n 2,25 64 
Ij 33,5 0,269 1080 
n 33,5 513 
4 rA 6,0 0,072 121 
r 6,0 98 
Ij 24,0 E 267 
5 n 21,5 165 
qj 1,4 = 90 
P. 2,0 
ii , 64 
Ij 36,0 0,273 1050 
SOLA, n 36,0 | : 493 
If 6,8 0,120 124 
ri . 68 88 


The decay curves for the phosphors impregnated in solutions Nos: 1, 2 and 5 
(Figs. 3, 4 and 7) cannot be represented by means of eqs. (3) and (4). For these phos- 
phors, Table III brings only the respective decay constants g ans the values of A. 
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Discussion 


Trypaflavine in gelatine is a two-phase substratum organophosphor, resembling 
a colloidal solution (Fröhlich and Mischung 1944). This dye, according to the hypothesis 
of Pringsheim et al. (Pringsheim and Vogels 1936, Pringsheim and Franck 1943, 
Pilipovitch and Sveshnikov 1958a, 1958b), should exhibit two metastable levels. 
Investigation of the decay curves for trypaflavine in gelatine impregnated in solu- 
tions of type 4 and 5 (Table II) throughout the range of low concentrations (10-75 to 
10-8 mol/1) yielded a shape of phosphorescence decay similar to to that observed at 
high concentrations 10-? to 107! mol/l) as shown in Figs. 5 and 6. The fact that for 
low as well as for high concentrations of trypaflavine in gelatine the decay curves can 
be decomposed into two exponential functions only points to independence of this 
decomposition from any interaction whatever between the excited luminescent mole- 
cule and other, non-excited ones, in its neighbourhood (Frackowiak 1957a, 1957b). 
According to Féerster (1946), trypaflavine in aqueous solutions forms associates. 
To determine whether such associates can be made to account for the decompo- 
sitions of the decay curves as observed in the present investigation, a comparison 
was made of the absorption spectra of aqueous and gelatine solutions of trypa- 
flavine. As seen in Fig. l, for aqueous solutions of trypaflavine, the changes 
characteristic of associates in the absorption bands do not appear (Kraviets, Pyeskina 
and Zitkova 1950, Levshin and Baranova 1956, 1958, Fóerster and Koenig 1951). 
Thus, it may be assumed that in the solutions used for impregnating gelatine with 
trypaflavine practically only monomers were present. 

A similar conclusion is reached on investigaung the changes in the absorption of 
trypaflavine in gelatine as shown in Fig. 2, with exception of curve No. 6. Even if 
trypaflavine dimers were to account for the "splitting" of this band, they would still 
have nothing to do with the observed composition of the decay curves, since the shape 
of the curves for the phosphor under consideration (Fig. 8) resembles that obtained 
with phosphors (Figs. 5 and 6) whose absorption bands exhibit no corresponding 
changes. N 

It is seen from Table II that the concentrations C* of trypaflavine in gelatine 
impregnated with solutions of different pH values are diferent, notwithstanding 
the fact that impregnation had been carried out in identical conditions. This 


- circumstance, and also the changes in the absorbtion bands of Fig. 2 point : 
to changes in the properties of gelatine resulting from immersion in solutions 
of different pH values. These are changes that can lead to variations in the 
shape of the phosphorescence decay curves. It will be shown below that the 
decay of the phosphorescence of trypaflavine in gelatine follows one general law 


enabling to describe all observed shapes of the decay curves from the properties of 


the gelatine. 


All the decay curves of Figs. 3, 4, 5, 6, 7 and 8 (and their decompositions) can be 


210 M. Frackowiak and H. Walery$ 


accounted for assuming trypaflavine decay in gelatine to be described by the following 


expressions: 


mO = 35,0 + 24) exp (Tu) 6) 


I, (t) = 2: B, (1 — r) exp (— £t) 


4 
I (t) = h (p) + 21, (t) = 3», B; exp (— T; t) (6) 
i=1 
with B, r 


„ and J’; (i = 1, 2, 3, 4) denoting certain constants. 

A number of simple hypotheses can be proposed which make it possible to des- 
cribe phosphorescence decay by (5) and (6). 

1. The trypaflavine molecules have two metastable levels and exhibit phospho- 
rescence within both phases of the gelatine. 

As well known the medium (substratum) does affect considerably the luminescent 
properties of dyes (Pringsheim ‘‘Fluorescence and Photsphorescence" 1951). If trypa- 
flavine is phosphorescent in both phases of the gelatine, then the phosphorescence 
observed is clearly a superposition of two ,,different” phosphors. Assuming both 
metastable levels to take part in the process of phosphorescence or in slow fluorescence, 
we are lead to expressions (5) and (6). Iwaki (1953a, 1953b) found that the constants of 
phosphorescence decay of trypaflavine and acridine orange in gelatine decrease as 
the exciting light intensity rises. Preliminary investigation of the pre-excited state of 
trypaflavine in gelatine by the present author showed the intensity of phosphorescence 
to rise with the time of pre-excitation. Variations of the intensity of the exciting light 
may be considered to produce effects similar to those of variations of the time of 
pre-excitation at constant intensity of the exciting light. Thus, e. g. for acridine 
yellow in gelatine, as pre-excitation rises, a fall in the phosphorescence intensity is 
observed simultaneously with an increase in its decay constants. The symbol P-'* 
has been employed for denoting such a pre-excited state of the phosphor (Frac- 
kowiak and Heldt 1959). It may be assumed that similar effects occur in trypaflavine 
in gelatine. In this case, however, the pre-excited state is Pt’, meaning that a rise 
in the intensity of phosphorescence is accompanied by a fall in the constants of 
phosphorescence decay. 

This hypothesis, however, does not take into account the effects of pre-excitation 
mentioned. To account for these, it would need functions of greater complexity than 
(5) and (6) to describe the decay of phosphorescence. As already mentioned, eqs. (5) 
and (6) are quite sufficient for description of all decay curves observed (their decompo- 
sitions into simple exponential functions). For these reasons the foregoing hypothesis 
can be rejected. 

2. The trypaflavine molecules have two metastable levels, and are phosphorescent 
within one of the gelatine phases only, and exhibit pre-excitation effects. 

On the assumption that acridine orange, like trypaflavine, has two metastable levels 
(Table I), we can expect the decay process of both dyes in gelatine to be similar. In- 
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vestigation on the subject has in fact proved (Frackowiak 1957b) the decay curves of 
acridine orange and yellow in gelatine to wholly resemble one another. Hence, we can 
compare the decay of phosphorescence in trypaflavine with that in acridine yellow, 
which exhibits pre-excitation effects (Frackowiak and Heldt 1959, Walerys 1959, 
Frackowiak, Grzywacz and Heldt 1959). The fact that the experimental decay curves 
of the phosphorescence of acridine yellow in gelatine (Frackowiak 1957a, 1957b) 
agree well with the theoretical curves for homogeneous rigid solutions (Jablonski 
1957a) indicates that acridine yellow does phosphoresce in one of the two phases of 
the gelatine. It is easy to state that the phosphorescence decay of trypaflavine in gel- 
atine differs fundamentally from that of acridine yellow. Hence also this hypothesis 
can be rejected. 

3. The trypaflavine molecules have one metastable level, are phosphorescent 
within both phases of the gelatine, and exhibit pre-excitation effects. This hypothesis 
yields an explanation of all decay curves as observed for trypaflavine in gelatine. 

Let indices i — 1, 2 in eqs. (5) and (6) denote non-deformed and deformed 
molecules within the first phase of the gelatine, respectively, and i = 3, 4 the respective 
molecules within the second phase. Moreover, let J; denote the respective decay 
constants, B; — the phosphorescence intensities at the time of cessation of pre- 
excitation, and r — the emission anisotropies. 

Investigations on the pre-excited state of acridine yellow in gelatine (Frackowiak 
and Heldt 1959, Walery$ 1959) have shown that the decay constants of the non- 
deformed and the deformed centres are different. Trypaflavine centres may be assu- 
med to have similar properties. Since these centres are situated within two distinct 
media, it may be assumed that I, Z T,Z D, I,, and moreover T} >I, and 
T, I4, because the pre-excited state is expected to be Pt’. 

Let 


I DC ml‘, 2n) 1"p2n,; "mü-— „len 
EA EN e oett (2) 
For (I, — 13) t <1 and (P, — D) t <1, eqs. (5) become 
mk I, +T 
IQ = B1. + 2r) (1 + mk) exp [- . o ai ) + 


nlr, +T. 
+ Bs (1 + 2r ar nl) exp-(— a j. 


m k I, + T. 
I (t) — B, (1 — 7i) (1 + mk) exp [- RAN ) “Ë 


n lT4 +l 
+ B, (1 — r3) (1 + n D) exp (a). (8) 


If the emision anisotropy of the deformed and non-deformed molecules differ but 
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insignificantly, r, = ry and rg e r4, then (cf. eq. (7)) m = m = 1, n = v = 1, and 


eqs. (8) assume the form 


Ij () = B, (1 + 2n) (1 + E) exp (a ) er 
£ B 190 D exp er), 
1,0 9-7) 0-8 (- FEED.) 4 
+ BQ — 1) (1+ Dex (— TD. (9) 
With 
A=Blth, nen =, 
A= B+), ner a (10) 


eqs. (9) reduce to eqs. (3). As already mentioned eqs. (3) account for phosphorescence 
decay in phosphors impregnated by means of solutions Nos. 3, 4 and 5 (Figs. 5, 6 
and 8). It may be concluded that in these phosphors the emission anisotropies of the 
deformed and non-deformed molecules differ insignificantly in both phases of the 
gelatine. 

The ratio of the phosphorescence intensities of the nondeformed and the deformed 
molecules is proportional to the ratio of their numbers, and thus inversely proportional 
to the probability of the formation of deformed molecules, A. As already mentioned, 
the substratum can affect various physical properties of the luminescent molecules, 
and thus, also the probability A. Changes in the properties of the substratum generally 
produce variations of the physical quantities mentioned. Hence, for the two distinct 
phases of the gelatine, A must be different: A, and A,. Since k ~ 1/4, and 1 ~ 1/45, 
and by eqs. (10), p, = gı (k), and g, = 9, (k), we have 9, = m, (À) and g, = 9 (Ay). 
From Tables III and II it is seen clearly that equal changes of the pH value of the 
solutions used for impregnating the phosphors give rise to different changes of g, 
and @,. Hence it is to be concluded that, when gelatine is immersed in media of 
different pH, the two phases change their properties in a different manner. These 
changes, in turn, affect the emission anisotropies r and r’ (cf. Tables III and II). 
The considerable differences in emission anisotropy r and r’ (Table III) can be explained 
by the assumption that, in addition to phosphorescence, trypaflavine at room temper- 
ature emits slow fluorescence which is relatively intense in one of the gelatine phases 
and weak in the other one. In this case rand r” should be considered to represent the 
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emission anisotropies of the total luminescence of the molecules present in the respec- 
tive media. | 

The decay curves of trypaflavine in gelatine impregnated with solutions Nos. 1, 
2, and 5 (Table II) can be described by eq. (8) only. In these cases, decomposition 
of the decay curves yields exponential functions of various decay constants. From 
eqs. (7) this seen to relate to situations when the emission anisotripies of the non- 
-deformed and deformed molecules are different. 

It should seem, however, that the quantities m, n and m’, n' (see eqs. (7)) are 
not related to the properties of the gelatine. From Tables III and II, the decay curves 
are seen to be described by eqs. (9) or (3) in all cases when the gelatine had been impreg- 
nated in solution for which y <1. This means that the emission anisotropies of the 
non-deformed and deformed molecules differ bit insignificantly if interaction between 
the trypaflavine molecules and the medium in aqueous solutions is relatively weak. 
This is in fact the case for solutions, for which 4 < pH < 10. The decay curves of phos- 
phors impregnated in solutions of pH outside the foregoing range can be accounted 
for only by eqs. (8). This is the case for phosphors exhibiting different emission 
anisotropies of the non-deformed and deformed molecules, r, Æ T2, T3 Æ Th 
(m Z m Z 1 and n Z n Z 1). These had been impregnated in solutions for which 
y > l, i. e. in solutions in which interaction between the trypaflavine molecules and 
the medium was very strong. This interaction could produce some insignificant 
structural changes in trypaflavine, resulting in insignificant changes in the luminescent 
properties of the non-deformed molecules but in highly important changes of the 
respective properties of the deformed molecules. It is known from various investiga- 
tions (Levis and Bigeleisen 1943) that phosphors consisting of dyes of closely similar 
structure exhibit greatly different properties in the pre-excited state. Hence it may be 
concluded that the quantities m and n (and m’ and n’) in eqs. (7) are related to varia- 
tions in the structure of the trypaflavine. molecule rather than to variations in the 
properties of the gelatine. 

~ Tt is seen from Fig. 3 that in the case of phosphor No. 1 the decay curves can be 
decomposed into three exponential functions only. This is the kind of decomposition 
to be expected in the present experimental conditions if one of the decay constants in 
one of exponential functions is, at least, of the order of 10° eos. i : 

In eqs. (5) and (6) fluctuations of the trypaflavine concentration in neither gelatine 
media were taken into account. Analysis of the decay curves shows that, for trypafia- 
vine in gelatine, the effects produced by interaction between excited and non-excited 
luminescent molecules (Jabłoński 1957 a, Frackowiak 1957 a) are yery small as compared 
to the influence of media on the trypaflavine phosphorescence in. two of differences 
phases of the gelatine. 'These effects can be accounted for by considering the use 
ties IT, B; and r, (i = 1, 2, 3, 4) as averaged over all the appropriate groups of lumi- 
nescent centres. In this case, the I’; and rj (i1, 2, 3, 4) should be expected to 
vary somewhat with the dye concentration (independently of the pH variations); 
however, the general character of eqs. (5) and (6) should not change. 
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Summary 


The authors investigated the absorption spectra of aqueous and gelatine solutions 
of trypaflavine. Monomers of the dye were proved to be practically the only forms pre- 
sent in both cases. The decay curves of phosphorescence of trypaflavine in gelatine 
coloured with dye in solution of different pH values were investigated. Three assump- 
tions, of the most simple kind, yielding a consistent account of all decay curves 
observed and of their composition, are discussed. The hypothesis by which trypaflavine 
molecule possesses one metastable level trypaflavine phosphoresces in both media 
(phases) of the gelatine, and exhibits pre-excitation effects, is proved to be the most 
plausible. All decay curves observed are proved to represent special cases within the 
framework of the foregoing hypothesis. 'The dependence of the decay curves on the 
pH values of the solution used for impregnation of the gelatine with the dye are account- 
ed for by changes in the properties of the gelatine. In the cases of solutions having 
pH > 10 or pH < 4, the shapes of decay curves are also influenced by persistent 
changes in structure of trypaflavine molecules produced by media of foregoing pH. 
The hypothesis proposed can be checked by investigating the pre-excited state of 
trypaflavine in gelatine and by studying the spectre of luminescence at various mo- 
ments following cessation of excitation. 

The authors wish to express their indebtedness to Professor Dr A. Jabtonski for 
his valuable discussions and for his interest throughout the present investigation. 
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The evolution in time of the hysteresis loop of (NH, CH,COOH),H,SO, was in- 
vestigated. Symples investigated immediately after preparation were found to present 
a double hysteresis loop. 

Only old samples can present double hysteresis loops. The term “old” is used for 
denoting samples that had not been subjected to investigation of any kind whatsoever, 
such as heating, an AC field, etc., for a long time. Young samples always exhibit a normal 
hysteresis loop. The process of converting an old sample into a young one is termed 
rejuvenation. It can be carried out by heating the sample above the temperature of the 
Curie point, or by placing it for a time in an AC electric field. 

The dielectric properties of triglycine sulfate doped with CuSO, were also investi- 
gated. In this case, rejuvenation in an AC field proceeds much more slowly. CuSO, doped 
triglycine sulfate samples present values of the polarisation P, and dielectric permitti- 
vity ¢ that are lower than the respective values for a sample containing no admixture. 

For doped samples, the coercive field depends on the amount of the admixture and 
rises with its concentration. 


Introduction 


Investigation of the ferroelectric properties of triglycine sulfate is due to Mat- 
thias, Miller and Remeika (1956). Results of investigations on the dielectric proper- 
ties of this compound have been published by Hoshino, Mitsui, Jona and Pepinsky 


Ger: : 
(1957), who found a spontaneous polarization of P, — 2,8 ES at 20°C. The Curie 


point lies at 48°C. A. C. Chynoweth (1959) found that defects from irradiation produce 
a decrease in the polarisation. The technology of obtaining large monocrystals of 
triglycine sulfate, together with results of investigations on their piezoelectric prop- 
erties (Konstantinova, Sylvestrova Aleksandrov 1959) and pictures of the domain 
structure (Konstantinova, Sylvestrova, Yureen 1959), are due to Konstantinova, 
Sylvestrova, Yureen and Aleksandrov, who also proved triglycine sulfate to exhibit 
memory of the stronger field. f 

Memory of the stronger A C field as investigated in BaTiO, by Rzhanov (1949) 
results from rejuvenation ‘with an A C electric field. Hitherto, ageing in ferroelectric 


(217) 
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materials was a process that had been known to take place in systems based on barium 
titanate. 

According to a hypothesis proposed by Pajak and Stankowski (1958), in an 
ageing ferroelectrics, domains of compensated spontaneous polarisation arise, and 
fully aged material is probably characterized by antiparallel configuration of the 
domains. This latter state has been termed domain antiferroelectric. Old samples 
exhibit a double hysteresis loop — a proof of the presence of the antiferroelectric 
state. A young sample always presents a normal hysteresis loop, and construction 
appears in the course of time. 

Yureen and Zheludev (1959) proved effects characteristic of ageing to take place 
in CuSO, doped Rochelle salt. A young sample of Rochelle salt containing an admix- 
ture of CuSO, exhibits a normal hysteresis loop which changes to a double loop in 
the course of 200 hours. The same authors found temperature hysteresis resulting 
from rejuvenation. Temperature hysteresis had been first observed by Piekara and 
Pajak (1952) in a system of the BaTiO, type. 


Method of measurements 


Triglycine sulfate crystals were grown from the saturated solution by gradually 
lowering the temperature (from 36° to 24°C). Triglycine sulfate with an admixture 
of CuSO, was obtained by the method of water evaporation at constant temperature 
(+8°C). The plates were cut perpendicularly to the ferroelectric axis (b-axis). The 
hysteresis loops were investigated with the Sawyer-Tower scheme. The dielectric 
permittivity was measured by the resonance method at a frequency of 1 MHz. 


Experimental results 


Immediately following preparation, triglycine sulfate exhibits double and 
asymmetric hysteresis loops (Fig. 1 a). If an AC electric field is applied to the sample 
for some minutes, the hysteresis loop is seen to assume the normal shape (Fig. 1b). 


Fig. 1. Oscillograms of hysteresis loops of pure triglycine sulfate: (a) sample immediately after produc- 
tion, (b) sample after some minutes of applying AC field 
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A more distinct rejuvenating effect of the AC field is shown in Fig. 2. Fig. 2a shows 
hysteresis loops in triglycine sulfate at a field of E = 600 V/cm. To C the 
sample, and AC field of E = 1,600 V/cm’ was applied for 5 minutes (Fig. 2b). It is 
seen that, on reducing the field strength to E = 600 V/cm (Fig. 2c), the polarisation 


Fig. 2. Oscillograms showing process of rejuvenation: (a) hysteresis loop of old sample, E = 600 V/cm, 
(b) after 5 minutes of applying a field of E = 1600 V/cm, (c) sample rejuvenated with AC field 
E = 600 V/cm j 


and losses for re-polarizing increase considerably. The sample is thus seen to have 
undergone rejuvenation by the AC field. 

By heating the sample above the Curie temperature (e. g. to 80°C) and subse- 
quent cooling, a young sample is obtained. The hysteresis loops for a thus rejuvenated 


C 
sample are shown in Fig. 3a, yielding a polarisation of P, — 2,74 E zand a coercive 
? cm 


field of E, = 400 V/cm. In the course of time, the loops become constricted, a symptom 
that ageing has set in (Fig. 3, 6 and c). 

Almost all samples exhibited asymmetry of the hysteresis loop. This asymmetry 
is related to the region within the crystal whence the sample had been cut (Fig. 4). 
The asymmetry of the hysteresis loop does not disappear in the process of rejuvena- 
tion. A similar effect has been observed by Holden, Merz, Remeika and Matthias in 
GASH crystals (1956). 
| If an admixture of CuSO, is added to triglycine sulfate, a material is obtained 
that is much less easily rejuvenated in an AC field. On applying a weak field to an old 
sample, the P(E) dependence was found to have the shape of a straight line (Fig. 5a). 
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Fig. 3. Oscillograms of the shape of the hysteresis loop as it changes in time, following rejuvenation 
at 80°C: (a) immediately after rejuvenation, (b) 260 minutes later, (c) 1330 minutes later 


SEE 
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Fig. 4. The asymmetry is related to the region whence the sample is cut. In both cases, the field was 
applied in the same direction 
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As the field intensity increases, a double loop and, subsequently, a constricted loop 
is observed (Fig. 5, b and c). The hysteresis loops obtained at the same field intensi- 
ties in a sample rejuvenated by heating above the Curie temperature are shown in 


Fig. 5, d, e and f. 


Fig. 6a. g (T) characteristics for sample A 


The e (T) dependence will be given here for three samples: A) one made of pure 
triglycine sulfate, B) one of triglycine sulfate with a small admixture of CuSO,, 
grown from a solution containing 0.05% CuSO,, C) one prepared from triglycine 
sulfate with a large admixture of CuSO,, grown a solution containing 0.25% 


CuSO, . 
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The e (T) dependence for samples A, B, C is shown in Fig. 6, a, b, c. Curve I 
corresponds to the old sample, whereas curve 2 is that of the young sample. The follow- 
ing table brings the values of the dielectric properties of sample A, B, C: 


TEC € peak 
Sample young sample (old) 
" 49.15 3010 
B 


75 ——1'C 
Fig. 6 b. e (T) characteristics for sample B. 
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Fig. 6c. e (T) characteristics sample C. 


The measurements of P, and E, were carried out at 20°C for young samples. 


Discussion 

The results of the present investigation can be summed up as follows: 

1. Triglycine sulphate undergoes ageing, since young material exhibits a normal 
hysteresis loop that becomes constricted in the course of time. 

. 2. Admixtures have the effect of fixing the domain antiferroelectric state, as in the 
cases of BaTiO, (Pajak, Stankowski 1958) and Rochelle salt (Yureen, Zheludev 1959). 

3. A larger admixture lowers. the Curie temperature T, of and old sample; the 
peak value of £ decreases, and so does that of P , Whereas the coercive field E, increases. 
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4. The asymmetry and electric properties depend on the region in the crystal at 
which the sample had been cut out. 

5. Temperature hysteresis in triglycine sulfate is somewhat different in shape 
from that occurring in doped BaTiO, (Piekara, Pajak 1952), and Rochelle salt 
(Yureen, Zheludev 1959), where a young sample exhibits a lower T, value than an 
old one (Yureen, Zheludev 1959), (Pajak 1957). 

6. In CuSO, doped triglycine sulfate, T, in the young sample is displaced 
towards higher temperatures, inversely as in BaTiO, and Rochelle salt. The value of 
e at the Curie point exhibited by a young sample exceeds the respective value for an 
old one. A similar effect has been observed in doped BaTiO, (Piekara, Pajak 1952) 
and Rochelle salt (Yureen, Zheludev 1959). 

It seems that all old ferroelectrics present the domain antiferroelectric state. 
Admixtures tend to fix this state. This is best seen in Rochelle salt. The latter presents 
antiparallel domain configuration; however, because of the low activation energy for 
reorientation, no double hysteresis loops are observed. Admixtures raise the activa- 
tion energy. On the other hand, an admixture homogeneously distributed throughout 
the material produces an increase in the coercive field, and does not lead to the ap- 
pearance of a double hysteresis loop. A young doped ferroelectric presents the usual 
hysteresis loop. It is only in the process of ageing that admixtures have the 
effect of fixing the domain antiferroelectric state. Some kind of specific ordering of the 
defects occurs here. It would seem that diffusion of the defects towards the walls 
separating the domains increases the activation energy for reorientation. This, 
probably, constitutes the microscopic mechanism of the process of ageing in doped 
crystals. In pure crystals, however, the ageing process, though not so pronounced, 
seems to be due to minor amounts of impurities or defects diffusing towards the domain 
walls. 

The authors wish to thank Professor Dr. A. Piekara for directing the present 
investigation and for his numerous valuable discussions. 
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ENERGY SPECTRA OF CHARGED PIONS PRODUCED 
IN pd-COLLISIONS AT 660 MeV* 


By H. HELFER, A. S. Kuznetsov, M. G. MESCHERYAKOV, W. SWIATKOWSKI, 
V. G. VOVCHENKO 


Joint Institute for Nuclear Research, Dubna 
( Received July 21, 1959) 


The energy spectra of charged pions produced in pd-and pp-collisions has been 
measured using a magnetic spectrometer. The differential cross sections for thep + d — nt, 
p + dm, and p+ pn" processes at 90? in the center of mass system of two 
colliding nucleons were found to be (5.0 + 0.6) x 10738, (0.57 + 0.08) x 107?8 and 
(6.7 + 0.7) x 107?5 cm?/sterad, respectively. The ratio of the efficiency of positive pion 
productions on protons in deuterons and on free protons is equal to 0.79 + 0.08. 


I. Introduction 


From the experimental data obtained at 460 and 660 MeV???) and at 970 MeV?) 
it follows that the proton-deuteron collisions without pion production are, mainly, 
a quasi-elastic scattering of incident protons by single nucleons in deuterons similar 
in general to the scattering of protons by free nucleons. Only in rare cases the incident 
protons interact simultaneously with both nucleons inside a deuteron. Such triple 
collisions involve also the elastic pd-scattering whose differential éross section at the 
energies 460 and 660 MeV and in the region of large scattering angles is equal to 
~ 10-29 cm?/sterad. Due to a proton-nucleon character of the high energy pd-colli- 
sions the possible reactions of a single pion procudtion are 


ptpt ira + pt n + (m), 0) 
p+ p+ () > m° + p + p + (n), (2) 
p + n + (p) > zz: +n + n + (p), (3) 
ptnrt(p>n + p + n + (B), fies) 
pnd(p-a ppt. (5) 


* Reported at the 6th Session of the Scientific Council of the Joint Institute for Nuclear Research, 
June 1, 1959. 
(221) 
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where (n) and (p) denote nucleons not colliding directly with incident protons, and 
therefore emitted with low energy. According to Batson et al.?) the probable energy 
of non-colliding protons in pd-collisions at 970 MeV is-about 5 MeV, their angular 
distribution being weakly correlated with the direction of the incident beam. In 
pd-collisions the production of pions is possible also in reactions in which two or 
three nucleons in final state are bound. These reactions can be as follows: 


p + p + () > mt + d + (n), (6) 
p + n + (p) > z° + d + (p). (7) 
pd mnt ?H (8) 
p+d>n?° + ?He (9) 


In the experiments described the production of charged pions in pd-collisions. 
was studied by measuring the energy spectra of positive pions from the reactions (1), 
(3), (6) and of negative pions from the reaction (5). The peak was also searched which 
would correspond to positive pions from the reaction (8). Besides, under conditions 
strictly identical to those of the observation of pions from pd-collisions, the energy 
spectrum of positive pions produced in free pp-collisions was measured. Such compara- 
tive study of high-energy pd-and pp-collisions can give informations on the effect 
of nuclear binding of nucleons in a deuteron on the process of pion production. 


2. Experimental Methods 


The present experiments were carried out with the six-meter synchrocyclotron 
of the Joint Institute for Nuclear Research. Use was made of the external unpolarized 
proton beam focussed by quadrupole lenses and passed through a deflecting magnet. 
Near the target, a thin-wall cylindrical container 4 cm in diameter and 12 cm high, 
filled with either liquid hydrogen or liquid deuterium, the proton density in the beam 
2 x 3 cm cross section was — 2.108 protons/cm? sec. The container was placed into 
a vacuum chamber cooled with liquid hydrogen. The external casing of the chamber 
was cooled with liquid nitrogen. It was supposed that the liquid deuterium and liquid 
hydrogen in the container were kept at the temperature of boiling hydrogen; their 
densities were equal to 0.169 and 0.0708 g/cm?, respectively). The effect of empty 
container was smaller than 3% of the particle emission in presence of liquid hydrogen 
or liquid deuterium. The proton energy in the center of the target was equal to 
(654 4- 5) MeV. 

The energy spectra of pions were measured by means of a magnetic spectrometer 
previously described by Azhgirei et al.°) Fig. 1 gives the schematic diagram of the 
spectrometer and its location with respect to the proton beam. In our experiments as 
well as in 5) two thin scintillation counters C, and C, placed at the input of the spectro- 
meter separated a particle beam emitted at the angle of 56? to the direction of the 
primary proton beam. This angle of emission of positive pions from the reaction 
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p + p — at + d corresponds to the angle of 90° in the center of mass system (c. m. s.) 
of the two colliding protons. The analysed beam on the horizontal plane was about 
+ 2°. The particles passed through channel II of the spectrometer were detected by 
a telescope constisting of four scintillation counters. The relative width Ap/p of the 
momentum interval selected by the scintillators was about + 3% for all measured 


M 
7 


" C] 
x: iron 


5 ///| lead 
Cz j 
C, EH aluminium 
ad uae Sd XP Cs B ^ 
Ce 
Fig. 1. Experimental arrangement. 1. Monitor. 2. Proton beam. 3. Target. 4. Magnet. 5. Filter. C,,.:.Cg 
scintillators. 


momienta, the subtended solid angle being about 0.5 >¢ 10-4 sterad. The threshold 
of pion detection was at 35 MeV. The dependence of a mean momentum of particles 
passing through the spectrometer from the magnet current was determined by cur- 
rent-carrying wire techniques. 

In the runs of the experiments the six-fold coincidences of pulses from the 
counters at the input and output of the spectrometer were detected simultaneously 
with four-fold coincidences from the counters of only one output telescope. The 
background of accidental coincidences was found by switching on the T lines 
shifting the pulses form separate pairs of counters C, Cs; Cs a C, Ce. In c pee 
the delay value the fine time structure of the synchrocyclotron brat: Rew was ta en 
into account. The detection efficiency of scintillations was close to 100% in all investi- 
gated energy regions. he small contamination of secondary protons was separated 
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from pions by means of a filter placed in front of the next to last scintillator of the 
output telescope. The relative contamination of j,-mesons and electrons in the particle 
beam passed through the spectrometer was determined by means of measuring the 
range curves in copper. To the results of measurements corrections were introduced 
which took into account decay in flight of pions as well as their slowing down and ab- 
sorption in the target, scintillators and in the filter separating the secondary protons 
from mesons. 

The absolute values of differential cross sections in spectra, d? o|d Q - d E, were 
estimated with the accuracy of about 7%, by comparing yields of positive and negative 
pions with that of recoil protons from the elastic pp-scattering. For 655 MeV at 56° 
in the lab. system (120° in c. m. s.), the differential cross section of the elastic pp-scat- 
tering was taken equal to (6.7 + 0.26) x 10-27 cm?/sterad (in c. m. s. (3.41 + 0.13) x 
x 10-2? cm?/sterad)9). 


3. Results of Measurements and Discussion 


The energy spectra of charged pions obtained from pd-collisions are given in 
Fig. 2 and 3 together with the statistical errors of measurements. For comparison 
in fig. 2 the spectrum of positive pions from free pp-collisions is plotted with a solid 
line without indicating the experimental points. This spectrum is measured with the 
same accuracy as the spectrum of positive pions from pd-collisions. All the rest of the 
errors connected with the determination of absolute values d? ojdQ - dE and with 


sterad. MeV 


700 200 500 E, MeV 


Fig. 2. Energy spectra of positive pions in the lab. system. The upper curve corresponds to free pp-col- 
lisions; the curve plotted according ta experimental points corresponds to pd-collisions; the dashed 
line — to pp-collisions in deuterons. 
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the estimation of corrections introduced were taken into account in calculating the 
differential cross sections d o/dQ of the observed reactions. Below 70 MeV the validity 
of experimental data on the form of spectra is not high, due to great difficulties with 
which the determination of u-meson contamination in this energy region is connected. 
'The usual transformation of the spectra to the c. m. s. showed that in this system the 
angles of emission of a great part of the pions lie in the range of 909—915. 
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Fig. 3. Energy spectrum of negative pions from pn-collisions in deuterons. 


When double pion production in one nucleon-nucleon interaction is rare, it is 
sufficient to consider that nagative pion production by protons on deuterons as 
occuring only in pn-collisions, while positive mesons are produced both in pp-and 
pn-collisions. Due to charge symmetry, the differential cross sections for the reactions, 
p+n>nat+n+nadprn>" +p +p are connected by the relation 


d š 
e pnm n + n) = ToU tam + p + p)180°-% 


where à is the angle of the pion emission in the c. m. s. At # = 90° the differential 
cross sections of the above reactions must be equal to each other, and the spectra of 
must be of the same form. Taking into account this 
inference from charge symmetry, one can expect that the result of subtraction of the 
spectrum of negative pions. produced in the reaction (5) from the observed total 
spectrum of positive pions of the reactions (1), (3) and (6) is the spectrum of positive 
-collisions in deuterons. The spectrum of positive pions thus ob- 
h a dashed line spreads towards high energies farther 


positive and negative mesons 


pions only from pp 
tained and plotted in Fig. 2 wit 
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than the spectrum of positive pions from free pp-collisions. This, apparently, is due 
to the motion of protons in deuterons. The absence of a peak from the reac- 
tion p + p— mt + d in the spectrum of positive pions from pp-collisions in 
deuterons, distinctly observed at 196 MeV in the case of free pp-collisions can be ex- 
plained in the same manner. Nevertheless, if one normalizes the spectra of positive 
pions produced on protons in deuterons and in free protons to equal areas, and then 
subtracts the contribution of the reaction p + p — zt + d from the spectrum of 
positive pions produced in free pp-collisions, it turns out that in the region of the peak 
corresponding to the reaction p + p > at + d the values d? ofdQ d E (ptp > at+ 
--..)p exceed the values d?o/dQdE (p + p > a+ + n + p)g. This fact 
indicates that in pd-collisions the reaction (6) is still contributing to the yield of 
positive pions. 

The absence of a peak corresponding to the reaction (8) in the spectrum of 
positive pions from pd-collisions is not a surprise as already at 340 MeV? the total 
cross section of this reaction is about 10-?? cm?, From the positive pion yield near the 
upper limit of the spectrum (E, = 310 MeV) where the peak of positive pions from 
the reaction (8) would be observed, one can conclude that if at 660 MeV the reaction 
(8) occurs its differential cross section for positive pion emission in the direction of 
56° to the proton beam does not exceed ~ 5.10-?? cm?/sterad. 

Comparing the spectra of positive and negative pions produced in pp-and pn-col- 
lisions in deuterons can see that to the right from the maximum the spectrum of 
negative pions falls more rapidly than that of positive pions. This difference, exceeding 
the limits of experimental errors, can be explained by the fact that in the reactions 
p+tp>n++n+pandp+n>n+p-+ p the effect of nucleon-nucleon inter- 
action in a final state is not the same. Indeed, in the first reaction the transition of 
T = 1— T = 0 of the initial and final isotopic spin states of the two-nucleon system 
leaves the neutron and the proton in 3S-state. Due to a strong interaction between 
‘nucleons in this state, the nucleons are emitted predominantly with low energies of 
relative motion and thus the energy taken by the pion increases. On the contrary, 
in the reaction p + n > m~ + p + p the transition T = 1 — T = | takes a consider- 
able part which leaves the protons in ?P state when the interaction between them is 
weak. 
The maxima of all the spectra (for pp-collisions only continuous spectrum from 
the reaction p + p > at + n + p is taken into account) are found practically at one 
energy, as is seen from the table in which the values of the mean meson energy in the 
spectra,are also given. As it should be expected, the obtained meson spectra differ in 
form from that calculated with taking into account of only one statistical weight of the 
final states. It is essential that the maxima in the spectra transformed to the c. m. s. 
are detected in the energy region closely corresponding to the energy of the maximum 
of z* + p-scattering cross section in the a+-p c.m. system. This fact indicates 
that in the processes of pion production in pd-collisions the resonance interaction be- 
tween a pion and a nucleon in the state with J — T — 3/2 plays an important role. 
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Energy of the | Energy of the Mean energy Mean energy 


Reaction maximum (lab. s.) | maximum (c.m.s.)| of meson (lab. s.) | of meson (c.m. s.) 


MeV MeV MeV MeV 


ptpoat+ptn 142 + 5 106 + 5 
p + d — mt + nud. 145 + 8 109+ 8 
p + d — z + null. 143 + 5 107+ 5 


Integrating the spectra in the c. m. s. over the energy, the following values of the 
differential cross sections of the observed reactions were obtained: 


do 
TET] (p + p mt + ...) = (6.7 + 0.7) x 10-8 em?/sterad, 


d 
de (p dor at + ...) = (69 + 06) x 1079 em?jsterad, 


d 
de (p+ d— ac) = (057 + 0.08) x 10 om°|sterad 


The differential cross section for positive pion production in free pp-collisions at 
660 MeV at 90? in the c. m. s. is close to the value (6,8 + 1.5) x 10-?8cm?/sterad 
obtained by Neganov and Savchenko® but distinctly smaller than the value (9.3 + 
+ 0.7) x 10-38 cm?/sterad reported by Meshkovsky at al”. 


The cross section z (p+ pz -c-Jp zi m je deut qe 


— E pt+a>m+..)=(53+ 0.7) x 10-28 cm?/sterad corresponds to the in- 


tegral yield of positive pions from the reactions (1) and (6). The ratio of probabilities 
of positive meson production in collisions of protons with free and bound protons is 


equal to 
do 
Jo Tp mt.) 
To = 0.79 + 0.08 
36 (P +p> "+ --)H 


As the spectra of pions from pd-and pp-collisions are measured in strictly identical 
conditions, the accuracy of determination of the value of this ratio is due only to statis- 
tical errors of measurements of pion yields and u-meson contamination. 

The fact that in case of pion production the protons in deuterons are less effective 
than the free protons is due to: (7) intranuclear motion of nucleons and their mutual 
“eclipses” in deuterons**; (ii) higher reabsorption of mesons by recoil nucleons in 


** According to Prokoshkin’s estimation (private communication) at the energy 660 MeV as 
a result of intranuclear motion the cross sections for meson production in nucleon-nucleon collisions - 
in deuterons must be about 4% smaller as compared to those on free nucleons. 


^ 
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case of pd-collisions; (iii) the effect of exclusion principle for recoil nucleons in pd-col- 


lisions. The ratio 
do ES "ES ( s: ) 
ig (P d pec =)D dQ De PR =P ose lee 


can in general differ from unity if only in processes (p + p> + ...)p and 
(p + p xt + ...),, the interference between amplitudes of possible meson-nucleon 
states is different. From the data obtained it follows that at 90° in the c. m. s. the 
relation of the number of positive pions to that of negative ones from deuterium is 
equal to 10.3 + 1.3. If in nucleon-nucleon collisions pions were produced only 
through the resonant state (J = T = 3/2), the ratio of total cross sections of positive 
and negative pion-production for pd-collisions would be equal to 11(19). 
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RANDWERTPROBLEME IM GEBIETE DER DURCH ELEMEN- 
TARGESETZE DEFINIERTEN LINEAREN SKALARFELDER IN 
ISOTROPEN UND HOMOGENEN RAUMEN 


Von A. RUBINOWICZ 


Physikalisches Institut der Universität Warschau 
Physikalisches Institut der Polnischen Akademie der Wissenschaften 


(Eingegangen am 9. September 1959) 


In der vorliegenden Note wird als das lineare Skalarfeld, das durch das Elementar- 
gesetz p = q (r) definiert ist, das Feld 


eco 
i5 Ny 9! o (r) y W 
ony Spitters eae iunt. 

a,b,c=0 

angesehen, wobei die Koeffizienten 4,5; beliebige komplexe Zahlen sein können. Es 
wird für solche Felder die Lösbarkeit von Randwertaufgaben untersucht, die durch die 


Vorgabe der Werte für 
OF JF QnaiF 

"an ONT UC Nn 
an der Begrenzungsflüche des Raumes definiert sind, in dem das Randwertproblem 
vorgegeben ist. Es wird bewiesen, daß ein solches Randwertproblem für endliche 
n-Werte nur in dem Falle gelóst werden kann, falls das Elementargesetz 9 — 9 (7) 
und daher auch das Feld F (x, y, z) Lösungen einer partiellen Differentialgleichung 
von der Gestalt 


a; Ai F (x, y, z) = 0 (a, Z 0) 


Ms: 


lI 


i=0 


sind. 


§ 1. Definition von linearen Skalarfeldern, die durch Elementargesetze erzeugt 
werden 


In den Zeiten als die Physiker noch nicht davon träumten, die Gesetze der Feld- 
physik durch partielle Differentialgleichungen darzustellen, da ja doch diese damals 
noch unbekannt waren, bediente man sich der Elementargesetze zur Darstellung 


; (235) 
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des Naturgeschehens. Es sei etwa an das Newtonsche Elementargesetz der Gravitation 
erinnert. 

Die Gesetze der Felder, denen wir in der modernen Physik begegnen, werden 
hingegen fast ausschließlich mit Hilfe von partiellen Differentialgleichungen festge- 
legt. Im Falle daß diese Differentialgleichungen linear sind, ergeben sich aus ihnen, 
je nach ihrer Ordnung, ein oder mehrere Elementargesetze, die gewöhnlich als die 
Grund- oder Fundamentallösungen bezeichnet werden. Auf diese Weise ist es stets 
möglich die Lösungen der linearen Feldgleichungen mit Hilfe von Elementarge- 
setzen darzustellen. Diese Möglichkeit besteht nicht nur im Falle der Felder, die 
durch partielle Differentialgleichungen vom elliptischen Typus beherrscht werden, 
sondern auch im Falle von hyperbolischen Differentialgleichungen. Auf diese Weise 
kehren wir im Falle linearer Felder eigentlich zum Standpunkte der Elementar- 
gesetze zurück. 

Wie aus einer früheren Arbeit des Verfassers (Rubinowiez 1951, wird im fol- 
genden als A zitiert) und den nachfolgenden Überlegungen hervorgeht, sind jedoch 
die mit Hilfe der Elementargesetze definierten linearen Felder allgemeiner als die 
durch partielle Differentialgleichungen darstellbaren. Angesichts der Schwierig- 
keiten, denen man heute in der Feldphysik begegnet, ist es deshalb vielleicht von 
Interesse die Frage zu stellen, wodurch die linearen Felder, die durch partielle Differen- 
tialgleichungen definiert werden können, sich gegenüber den nur durch Elementar- 
gesetze angebbaren auszeichnen. 

Wie man auf Grund von Elementargesetzen lineare Felder in einem isotropen 
und homogenen Raume festlegen kann, wurde in A für den Fall von skalaren Feldern 
gezeigt. Als das zum Elementargesetz @ = g(r) gehörige! in dem Quellpunkte Q 
(je nach dem Elementargesetz) reguläre oder hier auch Singularitäten aufweisende 


F eld haben wir das durch 


F (12,952): > Aabe an Ban (| — a + b + c) (1.1) 


a,b , c= 0 


gegebene Feld anzusehen. Die Koeffizienten Æ, sind hier komplexe Zahlen deren 
Wahl nur durch die Forderung eingeschränkt wird, daß (1.1) konvergiert. Die ein- 
zelnen in der Summe (1.1) auftretenden Ableitungen wollen wir nun als Felder von 
Gegenüberstellungsmultipolen? bezeichnen, weil man sie im Sinne der Maxwellschen 
Konstruktion von Multipolfeldern interpretieren kann. Da man das Quellsystem eines 
durch eine Ableitung /-ter Ordnung gegebenen Feldes als aus 2’ einfachen Quell- 
punkten bestehend ansehen kann, soll es das Quellsystem eines Gegenüberstellungs- 
multipols 2'-ter Ordnung heißen. 


1 Das Elementargesetz p = g (r) kann auch von der Zeit abhängen, also z. B. eine Kugelwelle 
f(r-et)/r darstellen, die einem beliebigen Zeitgesetze entspricht (vgl. A, § 7). 
2 In A wurden sie Felder mehrfacher Quellen genannt. 
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Außerdem haben wir in A den Begriff von Feldern eingeführt, die wir nun als 
Felder von Kugelfunktionenmultipolen® bezeichnen wollen Sie sind als die Felder 
(1.1) definiert, die sich in der Gestalt 


Yim (9, 9) R” (r) (1.2) 


darstellen lassen. Ein Feld (1.2) das eine Kugelfunktion Y,,, (2, y) I-ter Ordnung 
enthält, soll ein Feld eines Kugelfunktionenmultipols 2"-ter Ordnung heißen. 
Um ein solches Feld aus Feldern von Gegenüberstellungsmultipolen aufzubauen, 


à 9 9 9 
hat man sich eines Operators Y,,, pe Dy? = zu bedienen, der in der nach- 
stehenden Weise hergestellt werden kann: Bezeichnen wir mit Y,,, (z, y, z) die in 
_kartesischen Koordinaten dargestellte räumliche Kugelfunktion, die der flächen- 


haften Kugelfunktion Y,, (2, y) entspricht. Es ist also 
Y, (% y, 2) = r Ys, (0, 9). (1.3) 


Bekanntlich ist (1.3) ein in x, y, z homogener Ausdruck /-ten Grades, der der Poten- 
tialgleichung A Y, (x, y, z) = 0 genügt. Der gesuchte Operator wird dann erhalten, 
wenn man in Y,, (x, y, z) die Koordinaten x, y, z durch die Differentiationsoperatoren 
9/9x, 9[Oy, 2/92 ersetzt. Wie man etwa mit Hilfe gruppentheoretischer Methoden be- 
weisen kann, (vgl. A), besteht dann die Beziehung 


9 I 2 14-2s 7 f 
ee mall A (n) s R , wo Kram 
T € ° dy > 2) 4 9 (r) 14 (0, p) l (r) N ( ) 


; 1783. dec] M^ operto 
ni (rer (s 2j A:g(r =r = 2) [a c es g (r). (1.5) 


In‘ der von uns benutzten Terminologie stellt die rechte Seite der Beziehung (1.4) 
das Feld eines Kugelfunktionenmultipols 2!ter Ordnung dar, dessen Richtungs- 
abhüngigkeit durch die Kugelfunktion Y;, (0, g) bestimmt wird. Gemäß der linken 
Seite von (1.4) kann dieses Feld durch eine Summe von Feldern von Gegenüberstel- 
lungsmultipolen von der Ordnung 21 *5 dargestellt werden. Bei beliebig vorgegebenem 
Elementargesetz 9 = g (r) gibt es somit verschiedene Felder von Kugelfunktionen- 
multipolen 2"-ter Ordnung mit einer gegebenen Kugelfunktion Y,, (9, 9) je nachdem 
sie durch Felder von Gegeniiberstellungsmultipolen der gleichen sowie der beliebig 
höheren Ordnungen 2/*?, 2’+4, 2**, ... dargestellt werden. Sie unterscheiden ‘sich 
selbstverständlich durch die Radialfunktionen Rit (r) (1.5). 

In A wurden nun die Konsequenzen der Forderung untersucht, daß die Anzahl 
der linear unabhängigen Gegenüberstellungsmultipolfelder, die zu einem gegebenen 
Elementargesetz 9 = gf) und einer bestimmten Kugelfunktion Her Ordnung ` 
Yim (9, 9) gehören, nicht größer sein solle als n. Es wurde gezeigt, daß in diesem Falle 
ee INS NEED E ee ee SB: 

3 In A hießen sie einfach Multipolfelder. 


= 
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das Elementargesetz 9 = g(r) und daher auch das mit seiner Hilfe aufgebaute Feld 
(1.1) Lösungen einer partiellen Diffrentialgleichung von der Gestalt 


Da ATRE 92 0 (1.6) 
i=0 
sein müssen. 

Wir sehen somit: Stellt man an ein lineares, in einem isotropen und homogenen 
Raume definiertes Skalarfeld (1.1) die Forderung, daß die Anzahl der Felder von 
Kugelfunktionenmultipolen mit einer bestimmten Kugelfunktion Y,, (9, p) die Zahl n. 
nicht überschreitet, so müssen diese Felder der partiellen Differentialgleichung (1.6) 
genügen. Die Gleichung (1.6) ist somit die allgemeinste partielle Differentialgleichung 
in einem isotropen und homogenen Raume, der ein lineares Skalarfeld genügen muß, 
damit es nur eine endliche Anzahl, nämlich höchstens n Felder von der Gestalt (1.2) 
gibt, die eine bestimmte Kugelfunktion Y,, (0, v) enthalten und einem bestimmten 
Elementargesetz entsprechen. 

Wie in A gezeigt wurde, kann das Feld F (x, y, z) (1.1) statt durch Gegeniiber- 
stellungsmultipole auch durch Kugelfunktionenmultipole dargestellt werden. Es wird 
dann durch 


Fra) => > Vig (Or 9) X cin RI () (1.7) 


gegeben. 

Wird das Feld F (x, y, z) durch eine Lösung einer Differentialgleichung vom 
Typus (1.6) gegeben, so kann es, wie in A ausgeführt und bereits oben bemerkt 
wurde, je nach der Ordnung dieser Differentialgleichung im allgemeinen nur mit 
Hilfe mehrerer Elementargesetze beschrieben werden. Schon im Falle der Wellen- 
gleichung Au +k? u = 0 muß man etwa die beiden Elementargesetze e”"/r und 
sin kr/r verwenden, um nicht nur im Quellpunkte singuläre, sondern auch hier 
reguläre Lösungen zu erhalten. Im Falle der Differentialgleichung (1.6) von der Gestalt 
(A + k?) (A — x?) u = 0 ist man überdies gezwungen auch uen etwa die beiden 
Elementargesetze e+”/r und e~—”’/r zu benützen. 

Im Falle von Feldern F (x, y, z), die durch Differentialgleichungen (1.6) gegeben 
werden, treten, wie in A gezeigt wurde, in manchen Fallen auch Felder von Kugel- 
. funktionenmultipolen von der Gestalt Y; (0, p) S;*?* (r) auf, die im Quellpunkte Q 
regular sind, die jedoch nicht mit Hilfe der Operatoren (1.4) aus einem Elementar- 
gesetz p = g(r) herstellbar sind. Als Beispiel diene die Potentialgleichung Au = 0, 
die Lósungen von der Gestalt 


1 : 
Yim (9, 9) ps: ; (1.8a) 


und 


Yim (0, g) r (1.8b) 
aufweist. Während die Felder (1.82) mit Hilfe des Operators (1.4) aus dem Elemen- 
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targesetz g(r) = l/r erhältlich sind, muß man zur Gewinnung von (1.8b) ein anderes" 
in A im § 5 angegebenes Verfahren benutzen. 

Um eine Verallgemeinerung der Klasse der Felder zu erhalten, die durch Differen- 
tialgleichungen von Typus (1.6) dargestellt werden, wäre es daher folgerichtig den 
Begriff der zu Elementargesetzen gehörigen linearen Felder weiter zu fassen: 

1. Man müßte auch solche lineare Felder in Betracht ziehen, die durch Summen 
von Feldern (1.1) mit verschiedenen Elementargesetzen g = g(r) festgelegt werden. 

2. Man müßte neben den Funktionen Y,,,(0, y) Rj*?' (r) auch die nicht mit Hilfe 
der Operatoren (1.4) herstellbaren Funktionen Y,,, (9, p) Srt2s (r) in Betracht ziehen. 

Wir wollen jedoch hier die durch Berücksichtigung der beiden oben angeführten 
Umstände entstehenden Komplikationen vermeiden und definitionsgemäß das Feld 
F (x, y, z) (1.1) oder F (r, 0, p) (1.7) als das durch das Elementargesetz 9 = (7) 
gegebene Feld ansehen, Eine Änderung des Endresultates unserer Überlegungen 
dürfte dadurch wohl nicht eintreten. 


§ 2. Randwertaufgaben im Gebiete linearer Felder, die durch Elementargesetze 
definiert werden 


In der vorliegenden Note wollen wir zeigen, daß in der Klasse der linearen Skalar- 
felder, die in einem isotropen und homogenen Raume durch Elementargesetze bestimmt 
sind, sich die Felder die einer partiellen Differentialgleichung (1.6) genügen, auch noch 
durch eine andere Eigenschaft auszeichnen. Um sie aufzufinden stellen wir uns die 
bereits in A formulierte, dort aber nicht beantwortete Frage: Welche Bedingungen 
muß ein lineares, in einem isotropen und homogenen Raume durch ein Elementar- 
gesetz p = g(r) definiertes Skalarfeld erfüllen, damit man für dieses Feld ,vernünf- 
tige“ . Randwertaufgaben stellen kann, d. h. solche, deren Lósungen durch eine 
endliche Anzahl von Randbedingungen eindeutig festgelegt werden. Als solche 
wollen wir etwa die Vorgabe der Werte von 
əF QF op-1F (2.1) 


an den Begrenzungsflächen des Raumes wählen, in dem das Randwertproblem vor- 
gegeben ist. Mit N bezeichnen wir dabei die äußere Normale an diese Flächen. 
Wir nehmen an, daß der Raum, für den wir das Randwertproblem zunächst 
lösen wollen, durch den Außen- oder Innenraum einer Kugel gegeben wird, in deren 
Mittelpunkt sich der Quellpunkt Q befindet. Es könnte vielleicht scheinen, daß man 
aus einem so speziellen Beispiel keine weitreichenden Schlüsse ziehen könnte. In 
Wirklichkeit können wir jedoch schon von vornherein feststellen: Sollte es sich zeigen, 
daß Felder F (x, y, z) (1.1) oder F (r, 9, p) (1.7) existieren für die wir keine „vernünf- 
tigen“ Randwertprobleme für den Außen- oder Innenraum einer Kugel lösen können, 
so dürfen wir selbstverständlich nicht behaupten, dab solche Probleme für Bereiche 


von heliebiger Gestalt lésbar sind. 


F, 
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Die vorgegebenen Randwerte können wir auf der Kugeloberfläche nach Kugel- 


funktionen Y, (0, p) entwickeln. Es sei etwa 


= > Yun (2.2) 


=0 m=-I 


Entscheiden wir uns für den Innenraum einer Kugel so ist 9/9) N = 0/dr und aus (1.7) 
und (2.2) folgt dann daß 


AS A 
din = Dn qa R” | (2.3) 
wo rg den Radius der Kugel bedeutet. Mit Rücksicht auf (1.5) gilt nun: 
RI (r) = a (r) sa aod ie (r), 
d? š l = 1 z 21 -- 1 5 : 
Sa Rit? tpe — ) gi" ma T 1+1+2 (en Rite+2 (n. 
SRE EP ON pats He Per PREND obesse (2.4) 
Aus (2.3) und (2.4) folgt jedoch das Bestehen der Beziehungen: 
-X Cm RIS Gr 
T Lota dot t 1+1+2s 
dim Do To dim Tr ^ Clm Ria (ro). 
Au PIE L(+ 2) = Lan 
dim 9 To dim zd 2 din = Di Clm Rit (r 0); 
dm +... =F £, pede (ro). (2.5) 
s=0 ' 


Falls nur p Randbedingungen vorgegeben sind, so müßten bei festen l- und 
m-Werten die p Konstanten dí, (c = 0, 1, 2, ..., p — 1) eine unendliche Anzahl von 
Konstanten cj,, festlegen, vorausgesetzt daß 20 Funktionen Rj,* * (ro) ( zLhid-L 
l+ 2, ...) mit einem bestimmten /'-Wert fü? alle s-Werte von einander linear un- 
abhängig sind. Die Lösung des Randwertproblems würde also in diesem Falle nicht 
eindeutig bestimmt sein. 

Insbesondere müssen die linearen Gleichungen (2.5) auch im Falle / = m = 0 
erfüllt sein. Gemäß (1.5) ist jedoch in diesem Falle 


R$ (r) = 4g (r) (2.6) 
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Damit also insbesondere in der ersten Gleichung (2.5) im Falle 1 = m = 0 nur eine 
endliche Anzahl linear unabhängiger R7 (ry) (2.6) auftritt, muß zwischen ihnen eine 
lineare Beziehung bestehen, die wir stets in der Gestalt 


Dardi or) =0 (a, 0) (2.7) 


i=0 
voraussetzen können. 
Wir wollen uns zunächst überzeugen, daß beim Erfülltsein von (2.7) in der 
ersten Gleichung (2.5) im Falle / = m = 0 nur n linear unabhängige 


RF (re) = 4e 1,2. 670 1,2, sn 1) o 


auftreten. Mit Rücksicht auf (2.8) geht ja die erste Gleichung in (2.5) in unserem 
Falle über in 


d = 21d 4 90 | r= r (2.9) 
4 Mit Hilfe von (2.7) kónnen wir hier A” o(r)|, = ,, durch die 4: gr), - „mit i = Oma, 
nl ausdrücken. Wendet man auf (2.7) den Laplaceschen Operator 4 an, so kann man 
A" *! g(r) =n durch die A’ g(r) | rar, mit i = 1, 2, 3,..., n darstellen u. s. w. Auf 
diese Weise erkennt man, daß beim Bestehen der Beziehung (2.7) in der Summe 
(2.9) für dj, nur die ersten n Ausdrücke A’ g(r),..,, (s = 0, 1, 2, ..., n — 1) auftreten 
und daher diese Summe die Gestalt 


n—1 a 
doo T 2 Coo A’ 9(r)l,.- (2.10) 


annimmt. Die Koeffizienten Cy werden durch gewisse unendliche Summen gegeben, 
in denen die Koeffizienten cj, und a, auftreten. 

Besteht nun eine Beziehung von der Gestalt (2.7) zurecht, so reduzieren sich im 
Falle 1 — m = 0 alle Summen in (2.5) in der Weise, daß sie zu endlichen Summen 
werden, die jedoch alle die gleichen, nämlich die in (2.10) auftretenden Koeffizienten 
Ci, enthalten. Die lineare Beziehung (2.7) zwischen den A’ g(r) = Re (r) besteht 
nämlich unabhängig vom /-Wert auch zwischen den RÍ*? (r). Man braucht ja bloß 


den Operator Yin x zur 2 | auf die Differentialgleichung (2.7) anzuwenden, um 
Ox’ dy’ Oz] ^ f x 
mit Rücksicht auf (1.4) die Beziehung 
DR) (a, 2 0) (2.11) 
i=0 J 


für beliebige l mit den in (2.7) auftretenden Koeffizienten a; zu erhalten. 
Aus der Relation (2.11) kann man aber auch folgern, daß für jedes beliebige 


l, m-Wertepaar die Summen in (2.5) die Gestalt 


n-1 TE à 2 
TERR roe Haee bei oho s x2 (2.12). 


s=0 


040 A. Rubinowicz 


annehmen. Es treten somit bei gegebenem J, m in jeder Summe in (2.5) die gleichen 
Koeffizienten Cj, auf. Das lineare Gleichungssystem (2.5) geht daher in ein Gleichungs- 
system für die n Koeffizienten C7, über. Ist somit die Anzahl p der Randbedingungen 
(2.1) gleich n, so ist das aus den n Gleichungen (2.5) bestehende lineare Gleichungs- 
system im allgemeinen auflösbar. Es können somit p = n Randbedingungen (2.1) 
befriedigt werden, so daß das gestellte Randwertproblem für den Innenraum der 
Kugel lößbar ist. ; 

Eine Nichtauflósbarkeit des Gleichungssystems (2.5) ist nur für solche Werte ro 
des Kugelradius zu erwarten, für die die (2.12) entsprechenden Determinante Rit AN] 
verschwindet. Dies kann selbstverständlich nur für besondere Elementargesetze 
p = ọ (r) eintreten. Ein bekanntes Beispiel ist das Elementargesetz g(r) = sin #r|r 
der Wellengleichung Au + k? u = 0, wo n = l ist. 

Mit Rücksicht darauf, daß das Feld F (x, y, z) in der Gestalt (1.1) darstellbar 
ist, folgt aber aus (2.7) ebenso wie in A, daß es eine Lösung der Differentialglei- 
chung (1.6) ist. Unter dieser Voraussetzung ist das Randwertproblem (2.1) selbst- 
verstándlich nicht nur für die Kugel, sondern auch für einen Raum von beliebiger 
Gestalt lósbar. 

Die Tatsache, daB wir hier ebenso wie in A auf die Differentialgleichung (2.7) für 
das Elementargesetz p = g(r) stoßen, bedeutet daß die in A gestellte Forderung 
mit der hier formulierten äquivalent ist. Sowohl die Forderung, daß die Anzahl der 
linear unabhüngigen Kugelfunktionenmultipole mit einer bestimmten Kugelfunktion 
Yim (0, g) eine endliche ist, als auch die Forderung der Lösbarkeit der Randwert- 
aufgabe (2.1) führen zu dem Ergebnis, daß das Feld F (x, y, z) (1.1) oder das ihm 
äquivalente Feld (1.7) der partiellen Differentialgleichung (1.6) genügen muß. 


Interessant ist die Frage, warum wir hier sowie in A auf eine partielle Differen- 
tialgleichung stoBen, die gerade die Gestalt (1.6) hat. Diese Frage beantwortet sich 
dahin, daß dies bekanntlich die einzige lineare Differentialgleichung für: ein Skalar- 
feld im dreidimensionalen Raume ist, die invariant ist gegenüber der Gruppe der 
Verschiebungen und der Rotationen. Statt F (x, y, z) können in (1.6) bei den ein- 
zelnen Gliedern Summen von beliebigen zeitlichen Ableitungen der Funktion F (x,y,z) 
auftreten. 

Bezüglich der Feldtheorie der Elementarteilchen sei hier nachstehendes bemerkt: 
Die in der vorliegenden Note dargestellten Überlegungen sind im allgemeinen auf sie 
nicht anwendbar. Es hat ja doch in dieser Theorie keinen Sinn Randwertaufgaben 
zu stellen. Von diesem Gesichtspunkte aus waren somit hier Felder die nur durch 
Elementargesetze definiert werden können, zulässig. Man könnte dann jedoch eine 
solche Theorie nicht auf einer Langrangeschen Funktion aufbauen, da ja aus dem 


Lagrangeschen Formalismus zwangläufig als Feldgleichungen sich partielle Differen- 
tialgleichungen ergeben. 


Es ist ferner auch nicht abzusehen, wie man ein, nur auf einem Elementargesetz 
aufgebautes Feld quantisieren könnte. Daß skalare Felder von der Gestalt (1.2) 


M 
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gemeinsame Eigenlösungen der Operatoren des Quadrates des gesamten Impuls- 
momentes und seiner z-Komponente sind, bedeutet ja doch noch nicht eine vollstän- 


dige Quantisierung des Feldes. Es fehlt hier doch noch die Quantisierung der 
= Radialfunktion. 


L Krynica, im August 1959. 
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RELATIVISTIC HYDRODYNAMICS OF ROTATING FLUID 
MASSES MOVING WITH THE VELOCITY OF LIGHT 


Bv Prerre HILLION, TAKEHIKO TAKABAYASI, AND JEAN-PIERRE VIGIER 
Institut Henri Poincaré (Paris) 
(Received 13, October, 1959) 


Wir verallgemeinern eine urspriinglich von Weyssenhoff vorgeschlagene Methode 
und beobachten das Verhalten von relativistischen flüssigen Tropfen bei Geschwindig- 
keiten, die derjenigen des Lichtes gleich sind (oder nahe kommen). Es folgt daraus, 
dass für v = c neue qualitative Eigentümlichkeiten auftauchen und dass Tensor- Grössen 
durch Pseudotensoren, im Sinne Weyssenhoffs ersetzt werden. Als Beispiel geben wir 
die vollkommene, aus Trépfchen bestehende Flüssigkeit mit einer nicht nennenswerten 
Restmasse und zeigen dass, wenn wir sie mit einer einfachen Verbindung zusammen- 
bringen, die daraus entstehende Flüssigkeit genau der vor kurzem von Lee und Yang 
vorgeschlagenen, zweikomponenten Neutrinotheorie entspricht. à 


Introduction 


In two recent papers the general theory of relativistic fluid masses was treated 
from two different points of view. In the first, Bohm? and one of us (J.-P. V.) have 
established directly the physical signification of the well-known Weyssenhoff equa- 
tions of motion for such bodies. With the help of the new concept of center of matter 
density (which corresponds to the classical notion of geometrical body center in non 
relativistic mechanics) it can be shown that these equations represent the relative 
behaviour of this new center with respect to the usual center of mass. In the second 
paper? we gave (in collaboration with F. Halbwachs) a Lagrangian formulation of these 
equations with the help of a set of kinematical variables describing the rotational 
orientation of the element of the fluid which we call Einstein-Kramers variables. The 
utilization of these variables is valuable for various reasons which we have exposed 


‚elsewhere.? 


1 D. Bohm and J.-P. Vigier, Phys. Rev., 109, 1882, (1958). = 
2 P, Hillion, F. Halbwachs and J.-P. Vigier, Lagrangian formalism in relativistic hydrodynamics 
° of rotating fluid masses published in Nuovo Cimento, 10, 817 (1958). ne 

3 See ref. 2 and also F. Halbwachs, Thesis (Paris) to be Published at Gautier-Villars, (1959) under 
the title ,, Recherches sur la dynamique du corpuscule tournant relativiste et l'hydrodynamique des 
fluides à spin“. ; 
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In this paper, however, we want to study a special case of the latter formalism, 
namely the physical behaviour of rotating fluid masses for velocities approaching the 
velocity of light. 

This problem is interesting for three main reasons. In the first place, the study of 
the classical behaviour of extended bodies at such velocities might give useful indica- 
tions on the behaviour of matter in the high energy domain: in particular for temper- 
atures corresponding to such velocities. This should eventually be applied to the 
theory of collective excitation in high energy plasmas, and to the so-called cosmic 
“jets”. 

Secondly, it might. pave the way for a quantization of classical extended. bodies 
when v = c. ! 

Lastly, as we shall see, a perfect continuous fluid of a special type of relativistic _ 
fluid droplets with vanishing rest mass connected by simple tensions satisfies exactly 
the laws of motion of a two component neutrino field. 


We wish to stress, however, that independently of these applications the question 
of the behaviour of particles extended in space for velocities approaching the velocity 
of light is interesting as such. So far as we know it has never been treated except in 
a qualitative way.* This is only natural as long as one remains attached to the basis 
idea that particles are point-like and devoid of internal structure. 


Evidently in that case all that is needed classically is the description of the particle 
world-lines. On the contrary if the particles can be described by an energy momentum 
density T „and current density j, enclosed within a time-like tube it is clear that they will 
be modified by external motion. For example, the particles shape flattens perpendic- 
ularly to the velocity as a consequence of the Lorentz contraction and, at the limit 
v = c, is reduced to a plane section. For high acceleration its various points contract 


at different rates thus modifying the internal distribution of matter and energy. 


In this paper, we shall show that the kinematical variables introduced by Ein- 
stein and Kramers (which we have used to describe the internal states of the relati- 
vistic fluid masses) provide a convenient way to treat this problem mathematically. 
In section 1 we shall give a brief summary of the Lagrangian formalism in terms of the 
E. K. variables. In section 2 we shall see what happens to these variables when the 
velocity is increased up to the velocity of light. In section 3 finally, according to the 
preceding discussion, we shall study a particular case of relativistic fluid droplets 
which when accelerated up to the velocity of light and coupled in a simple way, be- 


have like neutrinos satisfying the two component spinor theory recently proposed by 
Lee and Yang.’ 


4 With the exception of an attempt by Weyssenhoff, Acta phys. Polon., 9, 7 (1947). 
5 T. D. Lee and C. N. Yang, Phys. Rev., 105, 1671, (1957). 

A. Salam, Nuovo Cimento, 5, 299 (1957). 

L. Landau, Nuclear Physics, 3, 127 (1957). 


Relativistic Hydrodynamics of Rotating Fluid Masses 247 


STE 


As is well known the study of extended relativistic fluid droplets rests on the 
idea that one can introduce integral quantities which correspond to average or global 
properties of the body. In this way the study of the internal complex motions can be 
left aside and replaced in the first approximation, by the laws of motion which govern 
these global quantities. 

We thus start? from the basic assumption that the particle is comparable to 
a fluid mass with conserved energy momentum density 7,, and current density J, 
which vanish outside of a time-like tube defining the physical limit of the droplet. 
One then defines immediately six physically important global quantities, namely: 
— The total energy momentum vector 


| 
Gu = ETC 
> 


(where du denotes a three dimensional volume element of a plane section 2) 
which determines the inertial Lorentz frame IT, (in which G, = 0). 
— The centers X and Y of mass and matter density defined in IT, by the relations 


XG (f Tha dug} = f 38 7% ae; 
> > 
and 
y) (f Ban) = f tas 
! > Le 
the latter moving along a world-line with the unitary tangent four velocity v, . The 


index „ denotes the fact that the quantities are evaluated in the frame IT. 
— The angular momentum M,, with respect to the center of matter density 


4 : 
My, — = [(x» x u) T,, ae (x E ») T,.] dv 
uc 
D 


Now, if we denote by A the derivation of a quantity A along the world-line fol- 
lowed by the center of matter density (A = v,9, A) we can establish the fundamental 
laws of motion 


G,=0 
M, = Guy — G, u, (1) 


6 See references (1). From now on we shall use the same notations. Greek indices u vary from 
one to four and denote the usual tensor indices — latin indices vary from one to three and denote space — 
like components — indices repeated twice are summed over all possible values — We work in Minkow- 


| gki's space time. ; 
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which combined with a Weyssenhoff-like condition (such as M,, v, = 0) determine 
completely the behaviour of the global quantities attached to the droplet. If instead 
of a single droplet one considers a perfect macroscopic fluid" made of such droplets 
one can define, instead of the preceding tensors, density relations obtained by multi- 
plying (1) by the scalar density o and adding to them the relation 


0 T3 9, (ov,) =0 


As was pointed out elsewhere, this procedure has the advantage of introducing 
into the formalism the mathematical methods of field theory and involves no disad- 
vantage since the droplets are assumed not to interact. 


The passage to Lagrangian formalism rests on the introduction of the Einstein- 


(5 
u 


(The index £ varies also from one to four though it does not mean a vector index) 
one of which a? is identical to v. This frame is attached to the center of matter density 
of each droplet which constitutes the macroscopic fluid, and its rotation represents the 


Kramers variables. They consist of an orthogonal set of four unitary vectors a 


instantaneous rotation of the droplet through the relations 
Oy, = $ (ay? aP — ai ay”) 


which determine the angular four velocity. 

We are now in a position to demonstrate the following general theorem®. If 
in the first approximation M,, is a general function of the a(? and not of their deriv- 
atives, the equations of motion (1) are just the conservation equations associated 
with the Lagrangian 


L = ko Map 0,5 + eM c? + ico a, Oat A, (a? ae uv) 


The physical signification of each term is clear. the first term $@M,,@,, COT- 
responds to the rotation energy, the second mz, c? to the rest mass energy, the third 
to the conservation of current and the last to the orthogonolity and unitary condi- 
tions of the a® variables. 

The demonstration of this theorem is simple. It has been established by Belin- 
fante and Rosenfeld that the canonical energy momentum tensor deduced from 


a Lagrangian L of variables q, namely: t,, = € 94 q, — Sue L 
e(d, qa) 
is conserved as a consequence of the field equations, namely 


I, t, | = 0 (2) 


If we can write o M,, as m,, where ma (a) is the average angular momentum 
of the droplets which constitute the macroscopic fluid at a given point (so that the 


7 In hydrodynamics this means a fluid without sources or sinks made of such elements without 
interaction; the center of matter density of each droplet following a given line of flow. 
8 Which generalizes the results of reference (2). 
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total Lagrangian is just the Lagrangian of the individual droplets 3 M,,@,,+ moc? + 
(b aG) ‘oli he I | 

e ed Gs ag! — 6,,) multiplied by o plus the Lagrange conservation condition 

icpa(? 9, S we get by variation with respect to o and using MU cec 


DIN 


As a consequence the energy momentum tensor /,, becomes: 


bay = Eu v, 


with g,=0(,S+%m, ae) 9, 2) = o G, and the corresponding conservation 
relation ð, ¢,, = 0 can be written: 


g, = 0 (3) 


In the case of density quantities, the dot implies, as it is well-known, the derivative 
along the local current line, that is: g, = 9, (v, &,) 
` The next step is to write the second Rosenfeld-Belinfante conservation relation: 


2, uvi Fri i (t,, Tw tyu) (4) 
where 


a War ee! 
fina pm G uy] 3 (8; Gr) q 


$€ representing the infinitesimal Lorentz transform on the variables a(? that is 


a =} (Sau 0p, = Cer Op) 


This gives immediately as a consequence of our assumptions 


ic 
uei = $9 Mu af 
and relation (4) becomes 
0 M = m mU RU RET Sy Up (5) 
(with v, ES ica®, ho = |M,gl and o Mep = Map). This demonstrates our theorem 
since relations (2) and (3) are evidently identical with equations (l). 
We notice further that any given form of m,; will then correspond to a particular 


type of internal motion. For example one could start from a Lagrangian satisfying 


Weyssenhoff’s condition 
L — o rige + ico a(98, S Fico hy aP aP ai + eA, (aa — du) (6) 
which we shall utilize later. In this Lagrangian: 


PM (at 
2 = Vaa aE Zep hô Euvap Vv Map (7) 


2. (3) G4) a (0) 
al? = eua aa; a ae 
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§ 2 

According to our program the next step is to see what happens to this formalism 
for velocities approaching the velocity of light. This could be attempted in various 
ways. The first and most straight forward is just to accelerate the particle in a given 
direction and apply the usual corresponding Lorentz formulas to all tensor quantities. 
We have attempted that but it leads to very complex calculations which finally give 
the same results as the more simple and powerful method which we shall now de- 
velop; so we shall not discuss it in this paper. 

This other method rests on a very simple remark made by Weyssenhoff. If one 
follows the world-line of any given particle (in our case the motion of the center of 
matter density), the most natural parameter is the propertime 7 on that line (—c?dr? = 
= dx, dx,). Now, if that line becomes tangent to the light cone dr = 0 and r loses 
all physical significance. This shows that relations (1) which depend on v must be 
transformed for v — c and that v is not a suitable parameter for such velocities. 
Accordingly Weyssenhoff proposed to drop r as parameter and to replace it by another 
parameter p on the world-line, such that p will still flow even when it becomes tangent 
to the light cone. The only restriction on p is to impose upon it the restriction 


vn dr. 
T manc 


so that the condition t’ — 0 implies dr = 0 which means the particle moves with 
the velocity c. One can then define the world-line followed by the center of matter 
density not by the unit vector: 


ER CD dxa 
ica, = — 
š: dr’ 
but by another vector: 
das 
aec 
dp 
so that: 
— n n 
W icT a. 
and 
We We = — CT? 
If 7' > 0 we then get 
Wa Ww, = 0 


and the world-line becomes tangent to the light cone. In order to avoid confusion it is 
necessary to stress the fact that w, is not an ordinary four vector but a four vector 
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depending on the parametrization; its four components transform like components 
of a four vector when the coordinates are transformed without change of parametri- 
zation, but they are all multiplied by a common factor dp/dp when the parametrization 
is changed from p to p. 

As Weyssenhoff remarked it is even possible in a given frame to take as para- 
meter the ordinary time t = x,/c. In this case we change the parametrization for every 
change of coordinate system in such a way that p remains equal to ¢ in the coordinate 
system of interest. Then w, becomes u, = x, [dt with the components? (u, ic). It is 
then evident that the four quantities u, do not form a four vector neither in the ordi- 
nary sense, nor a four vector depending on the parametrization (this would only be the 
case if p was a scalar parameter) but form instead a new geometrical object which we 
can call *pseudo-vector". The components of such a pseudo-vector transform like the 
components of vector but with an additional multiplication factor dt/dt = dx*[dx, 
when we pass from the frame 2 to a frame 2. The magnitudes of a pseudo-vector 
(u„u,)”, or any scalar product such as (Grua) behave like *pseudo-scalars" and are 
also multiplied by dt/dt in the case of a Lorentz transformation. As a consequence, 
the pseudo-vector u, has, in all possible frames, the ordinary velocity v as space com- 
ponents and c as time-component. In the limit v — c one gets evidently: 


Uy Ug = 0. 


Weyssenhoff's method for the passage to velocities approaching the velocity of 
light is now clear. In any equation one replaces derivatives with respect to the proper 
time t by the expression: 


d dt d 

dr dt dt 
and pass to the limit by pushing 7’ to zero. More generally one can write laws of 
motion in terms of any parameter p, satisfying the condition 7’ = dt/dp > 0 and 
pass to the limit t’ = 0. 


To illustrate this general method, let us treat briefly the problem of the limit of 
Weyssenhoff’s equations for a particle moving with the velocity of light. 


š Š dii mm 
Let us choose as parameter the ordinary time ¢ (r = =) and denote by f’the 


time derivative of a given quantity f so that: 


fa € -).f + was TEN 


If we now take for example the relation 


® The u designs a three vector. 


252 P. Hillion, T. Takabayasi and J. P. Vigier 


i b itt E 
ıt can eG WITLLOIL — 
dt dt 


When tv’ > 0 this implies 


Applying the same reasoning to the other equations (1) one obtains the follow- 


ing laws of motion: 


c =0 
M == G, ug — G; Ug 
Meg Up = 0 


1 . : 
If we start instead from G, =m, v, + a M,, v and multiply both members by 
T A we get: 
12 , 1 4 
G,T?-—mgu,t + a Mw w. 
The left hand side goes to zero with t’ and multiplying by the matter density ọ 
we find for v>c 


+ 1 + 
Ho Uu + — Mur us = Ü 
with 


, , 
Mo — 73 gu Un = Q Tig Ee 


This implies that we introduce a density u, of proper mass which goes to zero: 
with 7’ since it would take an infinite energy to push a non zero rest mass particle 
to the velocity of light. 


As a consequence the equations of motion for a fluid of Weyssenhoff particles 
can be written: 


@=0  g,-—0 mi, = Bi Uy BU 
Myv Up =O. g,u,—0 m, u, = 0. 
Of course, these relations are not independent and reduce to the set: 
„=0 (8a) M, =G,u,—G,u, (8d) 
G, u, = 0 (8b) M, u, = 0 (8e) 


py ^p 


2 O ERNO) (8) 
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However, this procedure of Weyssenhoff is rather cumbersome since it implies 
going to the limit of every equation of motion separately and in more complex cases, 
generalizing his simple model!9, it proves often very difficult. 

In order to avoid any trouble we therefore propose to generalize Weyssenhoff's 
idea in the following way: 

If we utilize a Lagrangian formalism in terms of the E. K. variables the behaviour 
of particles for v = c is clearly determined by the corresponding Euler equations. 
As a consequence if we can determine the limit of these variables we can determine 
the corresponding limit of L for v — c and deduce therefrom the corresponding 
laws of motion, which are just the limit of Euler's equations. 

This reduces our problem to two steps: 

a. Determine the limit of the E. K. variables for v — c; 
b. Write the corresponding Lagrangian. 

'The determination of the limit of the E. K. variables can be worked out with the 
help of Weyssenhoff's pseudo-vectors when p — t. If we recall the definitions (7) 
we can write instead of a the vector a — u,[ic that is 


(4) 


Dr 
a, —T a, 


so that cae a = 0 for an isotropic world-line (e = (by 
This relation can be written: 
a) a9 — af a + a) af? =0  (k—1,2,3) | (9) 
Let us now put p = t and work out the limit in that case. 


We have u, = {u, ic} a(? — 1 and we can replace a by LO =, 1} 
with the relation 6|® pP = — 1 


In the same way, instead of CSE 262605 tuvap Uv Maß 
we can write: 
a(O — ER PB t u,m 
E 20 cho pra aa gen 
that is: 
Bo a with 22 af = 0. (10) 
As we have: 
US a) = () 
we can also write: 
a(9 al? — 0. an 


Ber: 72 a 2 ranp E a E 


10 Considered for example. in ref. 2 and 3. 
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The three relations (9) (10) (11) evidently imply that 
a — à of! (12) 


where A is a constant. 
This results immediately from the fact that 


ial? VEF, io? o VE GP: 
i j 


and 


Cx ge. V > (a)? y » (a9)? — a a® >0 
: J 


which results immediately from Schwartz’ inequality. As one knows the equality sign 
is only possible when relation (12) is satisfied. 

As we shall see later this relation implies that the spin is either parallel or 
antiparallel to the velocity, when v— c. 

Let us now consider the limit of ae and ot These limits, which we shall call 
bP and b® satisfy evidently the relations 


pO p = 0 bD = 1 (13) 
pO p = 0 pO po — 1 (14) 
BED —=90  b@)b@ =O. (15) 


Let us denote by f7, (r = 1, 2, 4) (k = 1, 2, 3) the space components of the 
quantities z). The relations (13), (14), (15) become for v > c 


peo po po p = = bw) 
B® pO = bP bY = y» 
1) 402 1) p2 
BE pp = DP oP (16) 
1) ed 1 
BP ap = 1 — WP) 
2) A(2) _ 2 
pe pe = 1 — egy 
where the first member is the three dimensional scalar product. These five relations 
are equivalent to four independent relations. We still need four relations to determine 
b and b? completely. To get them let us pass to the limit for v — c of the definition 
of a(? , We have: 
3 4) (1) 4G 
a suu OOD oP. 
The time component of a becomes: 
a = ezip a a aW. 
Multiplying both membere by 7’ and taking the limit we get: 
dB Me 8415 9) b» io 


pe E po š (go A B2) 
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where the symbol A denotes the usual three dimensional vector product!!. As 
b® = A bY with A = constant. 
Then: 
A = BO (B A B®) (17) 
The space components of a“ become 
Gy a) (4) p) 7 (2) 4) 2 (1) p2 
B; SS Ej4jk boue DIS ne bj Lb. Oise Eijk4 de og ip 


= b (fi^ A Be) + bg S + OP BAB), 


Multiplying these relations by t’ and going to the limit we get the vector relation 
(where b? are numbers): 


BO = at (^ A Be) + BP BO AB) + BY (BAB) 
replacing go by Ago and noting ß® = 1, we find: 
AB = po ABO + BO A (09 po — b» ga». 
If we calculate now the product A (BO A B®) we find 
2 (BO AB) = BO ABO A B® + B® A [POA (09 BO — bW gay 
taking into account the well-known vector relation 
Pse RER RG 
the preceding relation becomes 
à (BO A po) = (B® . B®) gio — (BO - BO) go + 
= [oo (p? : p?) E bY (p? : p? peo = 
— (B® . B®) (bP fo — wp go. 
Taking Pato account relations (16) this equation can be written 
1 BO ABH) = ug up go — BO — egy BO + BP + (OP)? bi” — 


— (by? bP] BO — bP (bP go — p» £e A B® A B®) = — Bo + yp po 
so that = Fl » 
Mp po = Bm — 1(89 A 89) (18) 


If we eliminate 6“ between (17) and (15) we get 


Bp a = [BO — A GO A B) (9 AB) 


as 
po : (Bo A po) 20 
11 The dot denotes the three dimensional scalar product. 
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then 
ID = — (Bo ^ p) . (B Ag). 
Utilizing: 
A-(B RO) = CAN B) 
we get 


bP = — B® (BO A B® A gà) = 0 


and a similar calculation gives bP = 0. 
This essential result shows that for v — c the vectors go to system of 3 pseudo- 
vectors in the sense of Weyssenhoff. The relations 


Go gine Nin eT 3 

ge ae = 6% My = 0, 1, 2,3 
valid for v Æ c become for v — c%: 

pp bf? = 0, mo == ll, 24 4 

USD T a Lp 
and the system a(? tends towards the pseudo-vector system: 


ay > bp = (B®, 1) 


ay ADD (19) 
ay? > bf) = (B®, 0) 
ap. > de = (895, 0). 


The next step before we determine the limit of L is to replace the fluid density g 
which is equal to the ordinary density in the local rest frame of the fluid at a given 
point by a new quantity for the velocity of light. The new quantity must evidently 
also satisfy the limiting form of the conservation relations 


3 e (a) = 0 (20) 
as 
a® = 1 a^ 


this becomes 
0 4 
9, [£ a) = 0 


12 These relations are similar to the equations which define a unitary orthogonal system of ordinary 
vectors in three dimensions; but here they apply to pseudo-vectors at the velocity of light which we can 
now treat as ordinary vector variables in our Lagrangian formalism. 
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and if t’ — 0 we can assume that o/7' > S, Æ 0 so that the limit form of (20) is just 
ST or eco 20 (21) 


As S, depends of t = x,/ic the relation (21) must be a consequence of the equa- 
tions of motion and cannot be assumed a priori as in the usual case with v Æ c. As 
we shall see this follows from our treatment where there will be no Lagrange supple- 
mentary condition implying the conservation of Sy. 

As an example of step b. we are now in a position to pass to the limit of Weyssen- 
hoff’s Lagrangian (6). We define this limit as the sum of the limit of various terms 
of the initial Lagrangian and shall justify this procedure by showing that we get cor- 
jd the relations (8) which describe Weyssenhoff particles moving with the velocity 
of light. 


Starting from (6) 
L = my c? + ic p a 9, S + ico hy at? al? 9, aP + Au (a? af — ô) 
for v — c, we have (with m, = 0): 

Q = So T 
I) (4) 

Tid, 220, 
1 1 

aD => 6) (22) 
(2) (2) 

a =D, 


€) (D) (r) ur) 
a, Qy — Og b; by "E Ó; 


and, as here, according to a preceding remark concerning Sọ we do not need the 
Lagrange parameter 9, S the equation (6) becomes at the limit v — c. 


L = ic Sp hy b 019, b2 + A, (b b: — 55). (23) 


The deduction of relations (8) from the Lagrange function (23) is very simple. 
The canonical energy momentum tensor associated to L is just: 


Iw = 9 (9, b) 9, bi 


= ic hg S, b BM 9, bf? 


S ME 


= Eu u 
with: 
gu = So (ic ho bj? 9, bj). 
The conservation equation 9, t,, = 0 is then just (5) 


a 


g, = 0 
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Remarking now that Euler's equation with respect to the variable S, is 9/95 = 0 
we find 

Eu Uu = 0 

that is precisely (5). 

The relation Sọ = 0 results immediately from a simple combination of Euler’s 
equations, derived from the preceding Lagrangian. 

Variation of (23) with respect to b{® gives: 

Variation of (23) with respect to bP gives: 

h c bP 9, (So b4) + hoc So bY” — 2A, p = 0 (b) 
So if we multiply (a) by 5%, (b) by b® and add the result we find: 
2A, (bj b, — b? bj) + hac So (B B2 + BM b(1?) + hy c bf BM 8, (Sy BY) = 0 


that is: 
9, (So bP) == (e 


since the first and second term vanish because bí? b” = b® b = 0 and A,; is 
a symetric tensor. 

The last two relations m,, b* = Oand Myy = g, V, — g, u, will be obtained from 
the Belinfante-Rosenfeld tensor /,,,,. 

Indeed utilizing definition (4) we find: 


Im te hg Sa (D Ve = pr Ge ee d 
with 
mus = the So (bo? DOO — DOS 6); 
so that relations 
9, des > buy = Loy 
give 
m, Ze g, 5; 8 bi. 
The last condition m,, b; = 0 results immediately from the orthogonality rela- 
tions f; B, = 6, deduced from the variation of the 4;, s. i 
To finish this section we would like to make two remarks. 
The first is that as a result of the Lorentz contraction Weyssenhoff particles 
moving with the velocity of light flatten into disc shaped forms in the plane orthogonal 
to the velocity. The spin is then either parallel or antiparallel to the velocity. Taking 


into account the fact that 61 - 6? = 0 one sees immediately that the Lagrangian takes 
the simple form 


hos db db r 
L=— = = ES ae obi ai LÀ (o? bf? Es) 
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If we remark that the first term in the Lagrangian b b{® 9, bf? is just the scalar 
product of the angular velocity c on the four velocity b? we see that the equations 


of motion imply that this term vanishes so that the particle is reduced to a flat disc 
which does not rotate around its axis. 


§ 3. 


According to our program, the last section will include the comparison of the 
preceding perfect fluid constituted with Weyssenhoff particles with the hydrodynam- 
ical representation of the two component spinor theory of the neutrino® 13 (recently 
proposed by Lee and Yang). As we shall see, this comparison shows that the two are 
identical when we introduce a simple type of tensions between those particles. We can 
thus demonstrate the following theorem: 


Theorem 


A perfect continuous fluid of Weyssenhoff particles with negligible rest mass 
which move with the velocity of light satisfies with suitable tensions the two component 
neutrino equation: the two component spinor density defining all tensor densities 
needed to characterize the fluids behaviour. To prove this result let us recall briefly 
certain mathematical elements! concerning two component spinors. 

The first step is to establish a representation of 2 component spinors by tensor 
quantities. If one starts from a right handed coordinate frame a right handed neutrino 
field is represented by a right handed 2 component spinor g^ (in all that follows 
spinors are considered as c numbers). Such a spinor transforms under a proper 


Lorentz transformation 
Xu = yy X, Ds 
by 
p > ^ = A gf FD 
where A? is a 2 x 2 matrix defined by the relation | 
| (49)* e, At = a, c, (3); 
where 
o, = {0p 0) | (4)s 
18 Cf, Takabayasi, Comptes-Rendus, 246, 1010 (1958). i °: 
; 14 The transformation properties of the 2 component spinor and its associated bilinear quantities 
. were developed in: Takabayasi, Nucl. Physics, 7. 231 (1958) The general mathematical techniques 
to represent a two component spinor by tensor quantities is included in Takabayasi, Prog. Th. Physics, 


14, 283 (1955) and Takabayasi and Vigier, Prog. Th. Physics, 18, 573 (1957) (especially in the 
Appendix). 4 
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c, are just the Pauli matrices and o, = io, is the so-called time Pauli matrix with 


a of aed’ 
Pol VOMI 


indices * + T have here the usual signification and denote complex conjugation, 
hermitian conjugation and transposition. 
In the case of the neutrino we know that 9 obeys the wave equation: 


9, 9^ + 0,9, 9^ = 0 (5)s 


0,9, gy? = 0. 
The complex conjugate of (5) is: 
9, M p, = 0. : (6) 


Similarly, a left-handed neutrino field gf is represented by a different 2 compo- 
nent spinor (left-handed spinor) gf which transforms, under the proper Lorentz 
transformation (1), like 


gro of = AF gt Ms 
where 
Att c, A! = ap o, 
with 
o, = (— op 04) = ae 
4! is related with A? by: 
Af = (Ar) (8); 
The wave equation satisfied by @% is: 
9, 9* — 0,9, 9? = 0 (9); 
i. e., 
0,9, 9, = 0. (10); 
If we consider a transformation 
gi of = Dh (11), 
where J” is a 2-by-2 matrix defined by 
Ti = Dau Tt i ee hear Ae (12), 


(note that (+ = 97 r» 


then we can show that g% is a left-handed spinor. So we call the transformation (11) 
a chirality-conjugation. 
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Now, the only possible non-vanishing bilinear covariants formed with g^ are 
the following (in so far as we do not use derivatives): 


Si = gt a, g^ (13), 
and 


d ë, = ("J+ a, p = gfT I+ o, g? (14); 
S? is a real 4-vector, while £, describe 3 independent components of a self-dual tensor. 
In other words, writing real and imaginary parts of £, as: 


1 x 
P= (ED 


1 : 
Ag >= DE E = P) (15); 


1 Š Š 
we see that AU and A constitute the space and time components of a real anti- 
symmetric tensor &,,; i. e., 


Ap = Š; 
iAP = éna (Bb) ~ Q, 2, 3) (16); 


Note that each of S2, 2 and Al” transforms as a 3-dimensional vector for 
pure space rotation. 


We can now prove that eR and £, satisfy the following identical relations: 


Sd si —0 (17); 
£, Ë, = 0 (18); 
E, St = 0. (19); 


It can immediately: be shown? that the spinor g^ is equivalently represented by 
the set of quantities : 


Se | 
restricted by the conditions (17) — (19). Relations (17) — (19) are proved by using 
the formula concerning the components of Pauli matrices: 
(Op) a8 (oup = 20, a's 5 Sab Sag 


(see ref. (4)), and also the properties of J” matrice (12). 
In terms of A? and A, in place of &,, equations (18) and (19) mean that AQ, 
AO and S2 are mutally orthogonal. Relation (17) can be rewritten as 


(S3 = (S)? (S= Si 0) 
BEAD S E aa 
16 See reference (14). £ 
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Furthermore, we can prove that: 


a? = APY? = 9 


with 
AD aD ap 
AD aP 49 | = (89)? (20), 
S 5 5 
If we therefore introduce unit vectors b(? by 
JO up E (21), 


bt) (r= 1, 2, 3) span a 3-dimensional right-handed orthogonal axes; thus: 
b, b, = 4,, b; bk = Oy. (22); 


This completes the first step of our demonstration: the original spinor g? is now 
equivalently represented by the set of variables: 


AP} G=1,2,3) (23), 
i e, by ig 
bye (24), 
and it is now clear that any expression of ø? must be representable explicitly in terms 
of (24). For instance bilinear quantities including 1-st order derivatives are represented 
as follows: 
gpt 8,9 —8,9^* g^ = iS; bi 9, b; (26), 
and: 
kis 


pt 0,9, 94 — 9, 9^ * o, g^ — — 5 


S; Erim Of Oy bs. (27); 


The second step is to show how the equation of motion (5), is representable in 
terms of our variables (24),. First, we multiply (5), by @°* from the left and multi- 
ply (6); by g^ from the right, then add or substract the resulting equations, then 
we obtain: 


Io (Pt 9^) + 9, (ptt o, 9?) = 0 (30); 
and 
(°t By 9* —9,97* 9) + (9^* 9, ge —9,9** o, 9) —0 (FD; 
Equation (30) is nothing but 
9, S£ = 0 


15:6, 


Jo S; + 9, (S; bb) = 0. (32); 
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On the other hand, equation (31) is expressed by means of the reduction formulae 
(26) and (27), as 
BD 9, BP + $ eu Hf? 3, BO = 0 (33), 
We then proceed in the following way: first we multiply (5) by o; to obtain 
0,9% g +9; g^ + i (G, 9; g* Os 9, g^) —0 (34), 
Then we multiply (34) by 9?^* from the left, obtaining 
gt 0:99 94 + g^* 9,9! + i (Gt 04,9; 9 — 97 09,9) —0 — (34), 
then add or substract (34’) with its complex conjugate. The results are: 
9, (pt o; g^) + 9, (ptt 99) + i ((g^* 049; 97 — 9; 9% oy 9) — jl) = 0 9) 
and 
(pto; 2o 94— Io g^* o; 9?) + (g? 9; — 8, g? 9?) + i te) — (9) — 0 
(36), 


where the symbol (j, k) denotes the preceding term after permutation of the indices. 
By means of the reduction formulae (26), (27),, these are re-written as 


8, S1 + 9, (S$ bF) + S3 bP 9; Up = 0 37)3 
and 
ine 


The 6 real equations (37); and (38); are not all independent of (32), and (33). The 
relations, which are contained in (37); and (38), which are independent of (32), and 
. (33), are only the following: 


8, (Sd b2) = Sd bk Do bi (39), 
and | : 
93V (S45]) = 51 05.99 0p. (40); 


For instance, relation (37); results from (39), and (40), by the procedure 
(39), x DP + (40) x bP. 
In this way we have obtained 4 independent real equations: 
| (82)... (83), CMa (40)4) (a 


which are equivalent to the original spinor equation (5); (which also contains 4 real 
equations). 
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The last part of our demonstration is just to remark that our equations of motion 
(41), (together with the constraints (22),) result from the following Lagrangian density: 


he 


1 
5 So (bk Oo bk + — Sus bi 010p) + As (bk di — Ors) (42) 


L= 2 


From this L, we can derive the energy-momentum tensor and intrinsic angular 
momentum tensor of the field by the normal Belinfante-Rosenfeld procedure. 
In the spinor formalism, the Lagrangian is: 


hc 
L = — >s; (9^ ou Ip p? — 9, gt a, q^). (43); 


Using the reduction formulas (26), and (27), , (43), is rewritten, then we get just 
the first term in the right-hand side of (42). Adding 9,, (b; b — 6,,), we could thus at 
once reach the form (42)5. From this form the equations of motion (41), can be im- 
mediately derived by variation. 

The same procedure can be applied, of course, to the gë field!$. 

The physical interpretation of the Lagrangian (42), is clear if we take into account 
the vector identity: 


bi DRO, b] + Erim bk 9, dm = 3 Erim Or 91 OF, 
which results immediately from the orthogonality conditions on the br. This relation 
implies that if we denote by Ly, and Ly the Lagrangians (10—2) and (42—3), we have: 
Ly = Ly + ho Sg Erim bz 9; bs, 
= Ly Oy 


If we denote by ©,, the tensor ©,, = ho Sg Eru bR 9, bj that is the contribution of 
the ©,, term to the canonical energy momentum tensor. 

We see immediately that such a term can be interpreted as a tension term binding 
together the Weyssenhoff droplets (moving with v = c), constituting the neutrino 
fluid. Indeed ©,, satisfies the well-known property of ordinary tensions namely 


N 


The physical meaning of these tensions can be understood if we remark that the 
divergence 0, of the term 
UNE 
JT 


(namely the Belinfante-Rosenfeld angular momentum density) is zero. 


fi UA — Gi 


16 We-want to remark here that it is not at all astonishing to find that our particles are described 
by a two component spinor with a definite chirality. As we have seen only two possibilities arise for the : 
velocity of light: the spin is either parallel or antiparallel to the velocity corresponding to opposite chira- 
lities. This was also pointed out by Wigner, Rev. Mod. Phys., (1957). 
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Calling f,,, the term corresponding to Ly we have then: 


al == ak 3 
2 So m, m 9, pa 2 So 7: I b3 des 3 (t, = t5.) (44), 
which is a known consequence of the neutrino equation. 

Physically this means that the tension angular momentum is zero so that the 


orbital angular momentum (x, £, — x, t,,) is just equal to the variation of the partic- 
, . Cs 3 . B 
les’ angular momentum. Noting that t,, = g, b; + O,, this can also be written: 


m,, = g, b, — By ba + $ (Ow — On) 


which means that the variation of the total angular momentum (orbital plus spin) is 
just equal to the torque exerted on any given droplet by the tensions of the neigh- 
bouring fluid. We do not propose to discuss here any specific physical model which 
would lead to such an interaction but remark, however, that it is certainly one of the 
simplest possible types. Any other type resulting from a general ©,, (function of the 
9, b*) would in general contribute to the Fun, and 9, fua introducing supplementary 
angular momentum. 

As we said before, this result can be interpreted in the following way: Consider 
a continuous distribution of a two component spinor satisfying the neutrino's equa- 
tion (5). With the help of this spinor you can determine variables by the relations (13); 
to (15) which satisfy (42), and therefore characterize completely a continuous fluid 
made of Weyssenhoff particles with negligible rest mass moving with the velocity of 
light. The lines of flow tangent to 5% are followed by the center of matter density of 
the corresponding droplet and the three frame br which varies along that path de- 
termines completely the evolution of all physical properties (angular momentum, etc.), 
which characterize Weyssenhoff’s motion for such velocities. 

In other words, in this last section we have used a new set of parameters (9°) 
instead of the E. K. variables to describe the droplet’s behaviour. A similar idea had 
already been proposed by two of us (T. T. and J.-P. V.) to describe the motion of 
Weyssenhoff’s particles at ordinary velocities!”, Only in that case one uses 4 component 
normalized Dirac spinors v and Dirac's matrices y, The correspondence with the E. K. 
variables aš is given by a, = yt y, p; ds = V* y, Ys V and al + ia? = y* y, y. This 
has also been proposed by Gürsey!*. 


Conclusion 


As we propose to discuss in detail in another paper some of the physical conse- 
quences of the preceding results we will limit here the discussion to a few remarks. 
As we said before the first two solutions give a general method to treat the be- 
haviour of any continuous energy-momentum and current distribution enclosed 


1? T. Takabayasi (unpublished). Unal and J. P. Vigier. Comptes-Rendus. 245, 1787, 


1891 (1957). 
18 Gürsey. Nuovo Cimento 5, 784 (1957). 
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within a time like tube by using E. K. variables. This can be applied to many physical 
problems and an attempt to utilize this in the case of the beams of high energy acceler- 
ators is under progress at the Institute Henri Poincaré. 

The second remark is that we have obtained in sections 2 and 3 a concrete com- 
pletely clear relativistic model of the hydrodynamical representation of the neutrirto 
equation. It includes a simple “quantum coupling" corresponding to “quantum poten- 
tial" terms in the hydrodynamical representations of the other quantum equations. 
It shows that the neutrino equation can be interpreted 1) as governing the evolution 
of a probability distribution of neutrinos (with density Sp) as in the usual interpreta- 
tion of quantum mechanics or 2) as describing a real perfect continuous fluid of Weys- 
senhoff with zero rest mass moving with the velocity of light connected by the ten- 
sions @,, 3) as giving a real physical model of the behaviour of the real physical 
wave u associated with the neutrino in the causal interpretation of quantum mechanics 
proposed by de Broglie, Bohm and one of us (J.-P. Vigier)!?. 

The third and last remark is that such an hydrodynamical representation of the 
neutrino equation paves the way to many interesting developments. For example in 
a subsequent paper we shall demonstrate that a mixture of two such perfect fluids 
with a spin-spin interaction leads to hydrodynamical equations which when expressed 
in spinor variables correspond exactly to the non-linear equation recently proposed 
by Heisenberg and Iwanenko”®. 

We wish to express our gratitude to Professors L. de Broglie, D. Bohm and to 
Dr. F. Halbwachs for the interest shown in this work and for various valuable sug- 
gestions. Professor Bohm in particular has greatly contributed to the clarification of 
our model. 


Appendix 


To demonstrate that at the velocity of light the set of E. K. variables a; transform 
into the three orthogonal pseudo-vectors bj, we can proceed in the following way: 
we have the following relations: 


(à — Un DIE (4) „aß (iy (2) (3) „(4 
a, = i ay = Sch, Enygp 55 Uy,” == Eh Um) ay (1) 
and i 
4) 26. 4 a 1 
Derg Go ps OF does OU (u, y, ë, 7 — 1, 2, 3, 4) (2) 


we defined the parameter b’ 


ach Mya 
Pera fed oss (3) 


19 See for example: De Broglie „Une Interpretation Causale et non-linéaire de la Mécanique". 
David Bohm, Phys. Rev. David Bohm and J.-P. Vigier, Phys. Rev. The theory is also summarized 
. in Halbwachs, see ref. 3. 

2 That is y,9, P — y {PtP —(Wt y My} P = 0. 
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if we take: 
= bum a and a) become b’ a”, b' a® 
so: 
(E a) (b^ af?) = (^ aff) (6 a) = v W 
now if we define: 
= imn bap, p d b' aco (5) 
and we have: 
BD AD oa (6) 
» BO (7) 
it results immediately from Schwarz inequality that: 
QUIT 8 


where A is a real scalar. 
Let us choose: 


BP =1 P= (123 0) 
b is like a two dimensional unitary vector and 
4) p4 
BD D = (10) 
Let us prove that A = + 1. For this, let us start with the relation: 
2) (8 SE DE. 4) „A 
a a — af? aP = any a af? 
which multiplied by b' and going to the limit (taking into account the relation (8)) 
-3 (DRA (2) Aa — 4) „(A 
iA (aA — aP BP) = eus BO af? 
the space components taken out of this expression are: 


5 1 
iA (a? Bf? — af? BP) = ga BP ah? 


and, as we have: a = ¿a pf? = — il? = 
à (E AB) = a BO — a (11) 


where the symbol — denotes a three dimensional vector. 
From (11), we have: 


a2 (42 A BM) (a9 A 5M) = (ap pO — 40) (a B® — a) 


2? [(a® - a) (o. p — (a0 - 52] = a al p . p + GM . aO — 2a aX . BH 
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taking into account relation (2): 


40). aO = 1 + a aW G® . 30 = 1 + a a® 


a £ po T al aO» š po == a® 


and relation (10): 
it comes for the preceding expression: 4? = 1 
and finally: 


3 2 4 
P = + ipu? 


(1) 


Now to obtain the limits for a{ and a{? let us suppose b' sufficiently weak to 


enable us to write: 


b'a® = ip + B26 + O, (0) (12) 
bra® = BO + ib’? + O, (b”2) (13) 


where 6 and 6™ are space-like vectors, O,(b’?2) and O,(b’?) are negligible quantities. 
For b > 0, b'a®, and b'a(? are not isotropic vectors and their norm must be 
positive; of course, they are always orthogonal so we have 


105200 520 ” (14) 
TO S (Ü (15) 
iB (Oe id) 5.0 (16) 


Moreover 69 =0 ~as pP — 1 and: 


y al = 


I 
=: 


T 
oer 
we conclude: b’? oo (D - or 09 —( 


let us apply (12) and (13) to the relation: 


BA 1 4 
b a = © ep ay) ap. (ba) 


we obtain: 
TUE MI Act CTI a ape OD C) 

or: 

|o _ D, 

ip = Ewap a” ag? e (17) 

3 B 

5 — ie, s aP aP (9 (18) 
The relation (18) gives: 

— à = 39 . (qO) AFD) (19) 


of? = — o (a A a2), — a (5® Aa), — a) (a A 5). (20) 
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Now, let us start with the relation: 
Ha? — eus P a af? 
The same calculation brings to: 
— à = 6 . (a A a®) (21) 
of? ee e (a A a), = am (09 A a), HET a) (a? A 80), (22) 
As 6S = 0 the equation (21) has two solutions: 
&9 — Aa + j a (23) 
69 = 0 (24) 


Let us consider the first solution. Let us multiply (20), by a® and by a? taking into 
account relation (19); we have: 


69 - O = a a 9. aD = a (9 (25) 
Because the orthogonality conditions, when we put 23) into (25) we obtain the 
8 Y P 
system: 
2+ à (aP)? + naQ? a = af? i? 
haf? of? + p + pa aP = of? ap 6 
The solution is: 
a D, p=aP, PIH Ca (py en 
Now, with 68 it is possible to calcutate 5. It’s sufficient to multiply successively 
(22) by a“ and a®, this gives: 
jo. go = 59.29 — 0 
or f 
3 = v (a? Aa) (28) 
If we bring (28) into (19): E 
ger. o E (a A a?) : (a? N a?) = [(a® : a) (a? 4 a) = (a a? = 
= v [1 + (a)? + a] = 980? 


This implies: % = — 1 
and: 


SS ea A 42) (29) 
But the relation (14) writes: 
iA. 6% — 0 (30) 


y 
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and the relation (17) gives: 
pi. (a9 Aa) = 1 (31) 
We see that (29), (30), and (31) are not compatible. 


Now we have the second solution 6 = 0 


From (22): 
6) = — 69 (a( A a®); (32) 
which implies (16), moreover (14), (15), (31) are verified if e? < 0. But bringing (12) 


and (13) into the expression: 
(a a® — a® a) 2 = e, (UP) (b'a) 
a direct calculation gives: 
a Aa — 39 + à — (89 + 6%) (9 = — ID (aD Aa + p) 

let us multiply this result by be, taking into account the relation (31) and b.b — 1 
we have: š 

s = — 26@ or op = — 1 
that is what we intend. Besides putting (32) into (19), we find: 

(a A a?) : (a A a enl 
which implies 

(a)? + (a)? = 0 
or 
a = a? —0 


this is the resulte which we wanted to prove. 


ERRATA 


«Nuclear Charge Effect on the Free-Electron Model Energy States” S. Olszewski 
(Acta phys. Polon. 18, 107 (1959)) : 
Formula (13): At the beginning of each of the three last lines of the formula: 


°” 


instead of the “+” sign it should be “—” sign. 


Table III & IIIa, fifth column: instead of Z, = a should be Z, = = sixth 


1 1 
column: instead of Z,, = 3 it should be Z= m 


“On the Coulomb and Exchange Operators in the Free-Electron Mode = S, Oll 


szewski (Acta phys. Polon. 18, 121 (1959)) 
Formula (1b) should have the form: 


Xp) = e (fee un ab) 


Formula (æ), footnote 2, should have the form: 


C= lim GERNE cos = cos az —— E — dr dr = | l 
C TEF ir-ri s 


= 


"Ez rg + 25 + Siln— m)} + 5— in + m) Siti) 


A Formula (not numbered), coming after formula (11), should be 

EE os az é dz = T S (al) Š B $ up 

e: din. of (id it a be: a- 
— A wee: 


= 


E Pe 


x 
i 
y 
¿ te. < 
` "ue 
- = ` I 
1 f d wy. wn * 
d. 2 3 . i E pm y DE ee E sont 
: ' iš r # n uui 
' - fo ~a 
4 


- = ipt di i ose") kipa. x m F K. E 
ç 14 ie — imn 45, 4 (Cuy SQ DANS h TANG -— yal a i 5 hi ah F?” LA NC [s 
2 9e. i d $ Ë L AS x ° ^ - 


i Be ç 


« "Per 


= 


"m 


A af Ft M y 


. 


E ` — AS me a” 
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